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Abstract 



A geometric evolution equation is a partial differential equation that evolves some kind of 
geometric object in time. The protoype of all parabolic evolution equations is the familiar heat 
equation. For this reason parabolic geometric evolution equations are also called geometric heat 
flows or just geometric flows. The heat equation models the physical phenomenon whereby heat 
diffuses from regions of high temperature to regions of cooler temperature. A defining characteris- 
tic of this physical process, as one readily observes from our surrounds, is that it occurs smoothly: 
A hot cup of coffee left to stand will over a period of minutes smoothly equilibrate to the ambient 
temperature. In the case of a geometric flow, it is some kind of geometric object that diffuses 
smoothly down a driving gradient. The most natural extrinsically defined geometric heat flow is 
the mean curvature flow. This flow evolves regions of curves and surfaces with high curvature to 
regions of smaller curvature. For example, an ellipse with highly curved, pointed ends evolves 
to a circle, thus minimising the disttibution of curvature. It is precisely this smoothing, energy- 
minimising characteristic that makes geometric flows powerful mathematical tools. From a pure 
mathematical perspective, this is a useful property because unknown and complicated objects can 
be smoothly deformed into well-known and easily understood objects. From an applications point 
of view, it is an observed natural law that physical systems will move towards a state that min- 
imises some notion of energy. As an example, crystal grains will try to arrange themselves so as 
to minimise the curvature of the interface between them. 

The study of the mean curvature flow from the perspective of partial differential equations 
began with Gerhard Huisken's pioneering work in 1984. Since that time, the mean curvature flow 
of hypersurfaces has been a lively area of study. Although Huisken's seminal paper is now just over 
twenty-five years old, the study of the mean curvature flow of submanifolds of higher codimension 
has only recently started to receive attention. The mean curvature flow of submanifolds is the main 
object of investigation in this thesis, and indeed, the central results we obtain can be considered as 
high codimension analogues of some early hypersurface theorems. The result of Huisken's 1984 
paper roughly says that convex hypersurfaces evolve under the mean curvature flow to round points 
in finite time. Here we obtain the result that if the ratio of the length of the second fundamental 
form to the length of the mean curvature vector is bounded (by some explicit constant depending 
on dimension but not codimension), then the submanifold will evolve under the mean curvature 
flow to a round point in finite time. We investigate evolutions in flat and curved backgrounds, and 
explore the singular behaviour of the flows as the first singular time is approached. 
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CHAPTER 1 



Introduction 

The mean curvature flow is a well-known geometric evolution equation. The study of the mean 
curvature flow from the perspective of partial differential equations commenced with Huisken's 
seminal paper HHull on the flow of convex hypersurfaces. Since the appearance of that paper the 
mean curvature flow of hypersurfaces has been a lively area of study, and indeed continues to be 
so. Although this seminal paper is now just over twenty-five years old, the study of the mean 
curvature flow of submanifolds of higher codimension has only very recently started to receive 
attention. This thesis is concerned with the mean curvature flow of submanifolds of arbitrary 
codimension, and the main results we obtain can be considered high codimension analogues of 
some early hypersurface results due to Huisken. To give these high codimension results some 
context, we first briefly survey the relevant hypersurface theory. 

Let F : ^ iV"+'= be a smooth immersion of a closed manifold S, and H{p, t) be the mean 
curvature vector of := F(S(p), t). The mean curvature flow of an initial immersion Fq is 

given by a time-dependent family of immersions F : S x [0, T) — )■ N"'~^^ that satisfy 

{§iF{p,t) = H{p,t), pei:,t>o 
\f{;0) = Fo. 

The mean curvature flow equation determines a weakly parabolic quasilinear system of second 
order. We refer to the initial-boundary value problem (11.1) as 'MCF.' We advise the reader that 
we shall sometimes refer to MCF as an equation, and at other times, as a system. We also point 
out that by hypersurface or submanifold, we mean an object that has dimension greater than or 
equal to two. For the entirety of this thesis the reader is to assume that n > 2. Flows of space 
curves have been studied before, however the techniques are not the same (the Codazzi equation 
is vacuous for a curve). The main theorem of HHuU asserts that the mean curvature flow evolves a 
convex hypersurface of Euclidean space to a round point in finite time. Huisken's approach to this 
problem was inspired by Richard Hamilton's seminal work on the Ricci flow |Ha2], which had 
appeared two years earlier. Because the normal bundle of a hypersurface is one-dimensional, both 
the second fundamental form and the mean curvature can be very profitably viewed as essentially 
scalar-valued objects. The second fundamental form can be treated as a scalar-valued symmetric 
(1, 1) -tensor, similar to the Ricci tensor, and many of the techniques developed by Hamilton in 
his study of the Ricci flow can be used. The first crucial step in BHulll is to show that convexity 
of the surface in preserved by the mean curvature flow, and this is achieved by Hamilton's tensor 
maximum principle. After tackling the problem of hypersurfaces of Euclidean space, Huisken next 
went on to investigate the flow of hypersurfaces in a general Riemannian manifold, and slightly 
later, of hypersurfaces of the sphere. The Riemannian case showed that negative curvature of the 
background hindered the flow, whilst positive curvature helped. Although in this thesis we do not 
investigate the case of arbitrary Riemannian backgrounds, we mention that in [Hu2 | convergence 
results similar HHuU are still true provided the initial hypersurface is sufficiently positively curved 
to overcome the negative curvature of the ambient space. On the other hand, Huisken's results in 
IIHu33 are particularly relevant to some of the work in this thesis. Since the sphere has positive 
curvature this helps the flow, and in this case Huisken was able to attain convergence results when 
the initial hypersurface satisfies a non-convex pinching condition. The pinching condition we 
work with for submanifolds is very similar that of HHuBI . 

A feature of [ Hul l was that at the finite maximal time of existence, the entire hypersurface dis- 
appeared into a point at the same time. The 'roundness' of the point is made precise by magnifying 
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the hypersurface as the singular time is approached. This distinguishing feature is a manifestation 
of the convexity of the initial hypersurface. If this condition is relaxed and the initial hypersurface 
is only assumed to have positive mean curvature, then in general more highly curved regions will 
shrink faster than less curved regions, and a singularity will develop at some point before the en- 
tire hypersurface disappears. This naturally leads one to ask what are the possible limiting shapes 
of an evolving hypersurface as the (first) singular time is approached. It is customary to break 
up the kinds of singularities that can form into two categories depending on the rate at which the 
singularity forms. For the present discussion is suffices just to say these are called type 1 and 
type II singularities. It turns out that type I singularities are much easier to analyse than type II 
singularities, and in the type I case Huisken was able to obtain a complete classification. This 
was carried out in two papers, [Hu4 1, where compact blow-up limits were classified, and in [Hu5 |, 
which treated the more general complete case. A key element of this singularity analysis was the 
monotonicity formula introduced in IIHu4L 

Having briefly sketched the first developments in the study of the flow of hypersurfaces, we 
now turn to the study of the mean curvature flow of submanifolds, what is known and the results 
contained in this thesis. Much of the previous work on high codimension mean curvature flow has 
used assumptions on the Gauss image, focussing on graphical IICLT1ILL21IW21 IW41 . symplectic or 
Lagrangian submanifolds [ S WIICL21IW 1 [|Sm2l lN I . Another line of approach has been to make use 
of the fact that convex subsets of the Grassmannian are preserved [TW,W3,W6|. In this thesis 
we work with conditions on the extrinsic curvature (second fundamental form), which have the 
advantage of being invariant under rigid motions. Several difficulties arise in carrying out this 
program: First, in high codimension the second fundamental form has a much more complicated 
structure than in the hypersurface case. In particular, under MCF the second fundamental form 
evolves according to a reaction-diffusion system in which the reaction terms are rather compli- 
cated, whereas in the hypersurface case they are quite easily understood. Thus it can be extremely 
difficult to determine whether the reaction terms are favourable for preserving a given curvature 
condition. Second, there do not seem to be any useful invariant conditions on the extrinsic curva- 
ture which define convex subsets of the space of second fundamental forms. This lack of convexity 
is forced by the necessity for invariance under rotation of the normal bundle. This means that the 
vector bundle maximum principle formulated by Hamilton in HHaBII . which states that the reaction- 
diffusion system will preserve an invariant convex set if the reaction terms are favourable, cannot 
be applied. The latter maximum principle has been extremely effective in the Ricci flow in high 
dimensions BBWL BS, B2I where the algebraic complexity of the curvature tensor has presented 
similar difficulties. For arbitrary reaction-diffusion systems, the convexity condition is necessary 
for a maximum principle to apply. However, in our setting the Codazzi equation adds a constraint 
on the first derivatives of solutions that allows some non-convex sets to be preserved. As we 
have already mentioned, a similar situation arose in IIHu31 . where a non-convex condition was 
preserved. 

The content of this thesis is as follows. In the first chapter we summarise some standard 
facts on the geometry of submanifolds in high codimension from a 'modem' perspective. A key 
aspect of this is the machinery of connections on vector bundles, which we employ extensively in 
deriving the evolution equations for geometric quantities. In particular we introduce connections 
on tangent and normal bundles defined over both space and time, which prove very useful in 
deriving evolution equations and allowing simple commutation of time and space derivatives. This 
connection also provides a natural inteipretation of the 'Uhlenbeck trick' introduced in HHaBII to 
take into account the change in length of spatial tangent vectors under the flow. 

The second chapter fills in some details in the proof of short time existence for fully nonlinear 
parabolic systems of even order and applies this to the mean curvature flow. This a 'standard' 
result that is frequently quoted in the literature, yet a complete proof, especially in the setting of 
equations defined on a manifold, continues to remain elusive. In this regard, we draw attention 
to Lamm's Diploma Thesis, where he comprehensively proves local existence for fully nonlinear 
parabolic systems of even order in Euclidean space. We reconstruct some details of the following 
theorem: 
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Main Theorem 1. Let E x (0, w) be a vector bundle over M x (0, u), where M is a smooth closed 
manifold, and let U be a section T[E x (0, w). Consider the following initial value problem: 

{P{U) := dtU - F{x,t,U,VU, . . . ,V^'^U) = in E x (0,w) 
\u{M,0) = Uo, 

with Uq € C^'"'^'"(£'tj)- The linearised operator of P at Uq in the direction V is then given by 
dP[Uo]V = dtV + i-ir Yl A\x,t,Uo,VUo,...,V^"'Uo)ViV. 

\I\<2m 

Suppose that the following conditions are satsified: 



1) The leading coefficient A^^-'^' 



satisfies the symmetry condition ^f-*!-?!' ■ 



2) The leading coefficient satisfies the Legendre-Hadamard condition with constant A 

3) There exists a uniform constant A < oo such that X]|/|<2ml^^l";£^iJ — ^ 

4) F is a continuous function of all its arguments 

Then there exists a unique solution U € C'^"^'^'^ (E^), where /3 < a, for some short time > 
to the above initial value problem. Furthermore, if Uq and all the coefficients of the linearised 
operator are smooth, this solution is smooth. 

Chapter 3 contains what is the main result of this thesis, which is a high codimension analogue 
of Huisken's original theorem on the flow of convex hypersurfaces: 

Main Theorem 2. Suppose Sq = Fq{T,"') is a closed submanifold smoothly immersed in R"+'^. 
IfTiQ satisfies \H\,„i„ > and < c|-ffp, where 

if2<n<A 
ifn > 4, 



c < 



then MCF has a unique smooth solution F : S x [0, T) — )• R""*"*^ on a finite maximal time interval, 
and the submanifolds converge uniformly to a point q G M"^"'"*' as t ^ T. A suitably normalised 
fiow exists for all time, and the normalised submanifolds converge smoothly as t ^ oo to a 
n-sphere in some {n + k)-subspace ofW^^^. 

As the following simple example shows, the pinching ratio in Main Theorem |6] is optimal in 



n-li 



dimensions greater than or equal to four. Consider the submanifolds S 
where e is a small positive number. The second fundamental form is given by 

\ /O \ 



xS^(l) C 



h\ 



{ex,y) 



(x,0) + 



0/ 

_£{n-2) 
(n-l)2+e2 







V 



(o,y) 



V 



These submanifolds collapse to under 



and so they satisfy \h\ 

the mean curvature flow and do not contract to points. In dimensions two and three the size of 
the gradient and reaction terms of equation (14. 5t prevents the optimal result from being achieved. 
This is similar to the situation in [iHu3 J . where in dimension two the difficulty in controlling the 
gradient terms prevents the optimal result from being obtained. We remark that contrary to the 
situation in IIHu3ll . one cannot expect to obtain such a result with c = l/(n — 1) = lin the case 
n = 2 in arbitrary codimension as the Veronese surface provides a counter-example: This is a 
surface in that satisfies = but which contracts without changing shape under the 

mean curvature flow. We are not aware of any such counter-examples in dimension three (there 
are none among minimal submanifolds of spheres HCOl ). 

Curvature pinching conditions similar to those in our theorem have appeared previously in 
a number of results for special classes of submanifolds: In [01] Okumura shows that if a sub- 
manifold of Euclidean space with parallel mean curvature vector and flat normal bundle satisfies 
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|/ip<l/(n — then the submanifold is a sphere. The equivalent result for hypersurfaces 

of the sphere with < ^j^l-ffp + 2 (where the flat normal bundle condition is vacuous) was 
proved by Okumura in f02]. Chen and Okumura [C0 | later removed the assumption of flat nor- 
mal bundle and so proved that if a submanifold of Euclidean space with parallel mean curvature 
vector satisfies < l/(n — then the submanifold is a sphere (or, in the case n = 2, a 

minimal surface with positive intrinsic curvature in a sphere, such as the Veronese surface). The 
broad structure of the proof of Main Theorem 0closely follows IIHulll . which in turn, draws upon 
Hamilton's seminal paper on Ricci flow IIHa2i . 

After presenting the case of a Euclidean background we progress to discuss the situation where 
the ambient space is a sphere of contant curvature K. We obtain the following theorem, which can 
likewise be considered a high codimension analogue of [iHu3l : 

Main Theorem 3. Suppose Sq = Fq{Ti) is a closed submanifold smoothly immersed in S"^'^. If 
So satisfies 

\h\^<^\H\^ + '-^K, n = 2,3 
\h\^ < ^IHI"^ + 2K, n>4, 
then either 

1) MCF has a unique, smooth solution on a finite, maximal time interval < t < T < oo 
and the submanifolds St contract to a point as t ^ T; or 

2) MCF lias a unique, smooth solution for all time < t < oo and the submanifolds S^ 
converge to a totally geodesic submanifold Sqo- 

The assumptions of Main Theorem 2 required that |/?|mm > 0. In Main Theorem [3] no as- 
sumption on the size of the mean curvature is made, so the initial submanifold could, for example, 
be minimal. In this case the positive curvature of the background sphere still allows us to obtain 
convergence results. For similar reasons to the Euclidean case, the second main theorem is also 
optimal in dimensions greater than and equal to four. 

In the final chapter we follow Huisken's work in rHu41 and fHu5l to give a partial classi- 
fication of type I singularities of the mean curvature flow in high codimension. We pursue a 
slightly different blow-up argument than that used in IIHu41 and IIHu51l ; in particular, we consider 
a sequence of parabolically rescaled flows rather than a continuous rescaling. We also provide 
an alternate proof of the Breuning-Langer compactness theorem for immersed submanifolds of 
arbitrary codimension using the well-known Cheeger-Gromov compactness theorem. 

Main Theorem 4. Suppose F^ : SJ^ x (— oo, 0) — ?■ M"+^ arises as the blow-up limit of the mean 
curvature flow F : S" x [0, T) — )• M"+^ about a special singular point. IfT,o satisfies \H\,„in > 
and < 4/(3n)|-frp, then at time s = —1/2, FooC^oo) must be a sphere §"^(m) or one of the 
cylinders S*" (m) x M""™, where 1 < m < n — I. 

We close out the last chapter by showing how a simple blow-up argument can be used instead 
of the convergence arguments of Section 14.71 of Chapter |4] to determine the limiting spherical 
shape. 

The results of Chapters 2 and 4 appear in the paper 'Mean curvature flow of pinched submani- 
folds to spheres', which is coauthored with the author's PhD supervisor, Ben Andrews. This paper 
has been accepted to appear in the Journal of Differential Geometry. 
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Submanifold geometry in high codimension 

In order to work with the normal bundle we first discuss vector bundles, including puUback 
bundles and sub-bundles. The machinery we develop is useful and new even in the codimension 
one case, as we work with the tangent and normal bundles as vector bundles over the space-time 
domain, and introduce natural metrics and connections on these. In particular, the connection we 
introduce on the 'spatial' tangent bundle (as a bundle over spacetime) contains more information 
than the Levi-Civita connections of the metrics at each time, and this proves particularly useful in 
computing evolution equations for geometric quantities. 

2.1. Connections on vector bundles 

2.1.1. Vector bundles. We denote the space of smooth sections of a vector bundle E by 
r{E). If £^ is a vector bundle over A^, the dual bundle E* is the bundle whose fibres are the dual 
spaces of the fibres of E. If Ei and E2 are vector bundles over N, the tensor product Ei (g) E2 is 
the vector bundle whose fibres are the tensor products (£'i)p ® (-£'2)p- 

2. 1. 1. 1. Metrics. A metric 5 on a vector bundle is a section of E* (g) E* which is an iimer 
product on Ep for each pin N. A metric on E defines a bundle isomorphism =f^g from E to E*, 
defined by 

for all ^,r] e Ep. If g is a. metric on E, then there is a unique metric on E* (also denoted g) such 
that the identification #g is a bundle isometry: For all ^, 77 G Ep, 

If 5i is a metric on K„ i = l,2,thenc/ = gi(g>g2 G T({E'l(g)E'^)(g){E^(g)E^)) ~ r{{Ei(g) E2)* (S> 
{El (g) E2)*) is the unique metric on Ei (g) E2 such that g{^i (g) ^2 <^ m) = 6)5'2(??i, 

2.1.1.2. Connections. A coimection V on a vector bundle E over is a map V : r{TN) x 
r(^^) — > T{E) which is C°°(Ar)-Unear in the first argument and R-linear in the second, and 
satisfies 

vt/(/o = /vc/^ + (f//)e 

for any U G r(rA^), ^ G r(£;), and / G C~(iV). Here the notation Uf means the derivative of 
/ in direction U. Given a coimection V on E, there is a unique coimection on E* (also denoted 
V) such that for all ^ G r{E), oj G r{E*), and X G r{TN), 

(2.1) X{coiO) = i^xcom + oj{VxO- 

If V* is a connection on Ei for i = 1, 2, then there is a unique connection V on £^1 (g E2 such that 

(2.2) Vx(6 ® 6) = (v3f6) <8 6 + 6 <8 (v|6) 

for all X G r(riV), ^ier{Ei). in particular, for S e T{El ^ E2) (aa ^^a-valued tensor acting 

on £"1), VS G r(T*iV (g) El (g) E2) is given by 

(2.3) (Vx5)(0 = "^iHSm - S{Vl^i). 

A connection V on £J is compatible with a metric g if for any ^, r/ G T{E) and X G V{TN), 

If V is compatible with a metric g on E, then the induced connection on E* is compatible with 
the induced metric on E*. Similarly, if Vj is a coimection on Ei compatible with a metric gi for 
i = 1,2, then the metric g\ (g 52 is compatible with the coimection on Ei x E2 defined above. 
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Another important property of connections is that they are locally determined. 

Proposition 2.1. Let E be a vector bundle over N and p a point in N. If^i and ^2 cif^ two section 
of E such that ^1 = ^2 on an open neighbourhood U ofp, then 

for all X £ T{TN). 

Proof. It is obvious from the definition of a connection that the covariant derivative only 
depends on X at the point p. To show that it depends locally on ^, let p be a smooth cut-off 
function with support in U. Then p^i = p^2 on U and hence Vxp£,i{p) = ^ xP^2- Futhermore, 
the Leibniz property of a connection gives 

Vxp6(p) = {Xp){p)ii{p)+ p{p)Vxii{p) = Vx^p). 
The a same holds for VxpC2{p) too, thus Vx^iip) = ^x£,2{p) as stated. □ 

2.1.1.3. Curvature. Let E' be a vector bundle over N. If V is a connection on E, then the 
curvature of V is the section e T {T* N T* N E* E) defined by 

Rv{X,Y)C = Vy(VxO - Vx(Vye) - ^[Y,X]C 
The curvature of the connection on E* given by Equation (12.11 ) is characterized by the formula 

= {R{X,Y)u){()+uj{R{X,Y)O 

for all X,Y e r{TN), u E r{E*) and ^ e r{E). 

The curvature on a tensor product bundle (with connection defined by equation (12.21 )) can be 
computed in terms of the curvatures of the factors by the fonnula 

i?v(^,n(ei®6) = {Rv^{X,Y)^i)^^2 + ^i^{Rv<X,Y)^2)- 
In particular, the curvature on E^ (g) E2 (£'2 -valued tensors acting on £'1) is given by 

(2.4) {R{X,Y)Sm = R^2{X,Y){S{0) - S{R^,{X,Y)0. 

2.1.2. Fullback bundles. Let M and N be smooth manifolds, and let £ be a vector bundle 
over N and / a smooth map from M to A^. Then f*E is the pullback bundle of E over M, which 
is a vector bundle with fibre {f*E)x = Ef(^xy If ^ e r(£'), then we denote by the section of 
f*E defined by ^f{x) = Cifi^)) for each x G M (called the restriction of ^ to /). 

The pull-back operation on vector bundles commutes with taking duals and tensor products, 
so the tensor bundles constructed from a vector bundle E pull back to give the tensor bundles of 
the pull-back bundle f*E. In particular, if g( is a metric on E, then gis a. section of E* (>i) E*, and 
the restriction gj £ r{f*{E* (g) E*)) ~ T{{f*E)* (g) (/*£)*) defines a metric on f*E. 

Proposition 2.2. IfV is a connection on E, then there is a unique connection -^V on f*E, called 
the pullback connection which satisfies •^V„(Xj) = V f^uX for any u G TM and X G r(£'). 

Proof. Suppose that ^ is an arbitrary section ^ G T{f*E) and p £ M. Let Zi be a local 
frame for E about f{p). The sections {Zij} thus form a local frame for f*E about p so we can 
write ^ = ^^{{Zi)f. The properties of a connection and the pullback then give 

= v{e){z,)f+e^v,{z,)f 
= v{e){z,)f+evf,,z,. 

A further computation shows that this is independent of the local frame used, and because connec- 
tions are locally defined by Propostion 12.11 the pullback connection is well-defined. □ 
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At first one might think that for the puUback connection to be well-defined, it would be nec- 
essary to extend the section f^v to a neighbourhood of f{p) in order to operate on it locally. The 
above proposition shows that this is not necessary, although in order to define the pull-back con- 
nection, we had to define it terms of a local frame. Often in submanifold geometry the induced 
connection is defined in terms of a projection in N onto the image of the tangent space of M. This 
definition is frame independent, however it is necessary to extend the vector fields in order for the 
definition to make sense. One can then show afterwards that the definition is independent of the 
extension used. 

Proposition 2.3. If g is a metric on E and V is a connection on E compatible with g, then is 
compatible with the restriction metric gj. 

Proof. V is compatible with g if and only if Vg = 0. We must therefore show that ^Vgf = 
if V5 = 0. But this is immediate, since ^Vv{9f) = ^ f^^vd = 0. □ 

Proposition 2.4. The curvature of the pull-back connection is the pull-back of the curvature of the 
original connection. Here Ry G T{T*N T*N (g) (g) E), so that 

f*{Rv) e r(r*M (8) T*M (g) f*{E* (g) E)) = T{T*M ® T*M (g (/*^)* (g f*E). 

Proof. Since curvature is tensorial, it is enough to check the formula for a basis. Choose a 
local frame {Zp}^^^ for E. Then {{Zp)j} is a local frame for f*E. Choose local coordinates 
{y"} for N near f{p) and {x*} for M near p, and write f"' = y"'o f- Then 

Rf^{di,dj)iZp)f = fV9,i^V9AZp)f) - (i o j) 

= ^Vj{VudAZp))-{i^3) 
= N .{diPV aZp) - {i ^ 3) 

= {djdif'')VaZp + dif''Nj{{VaZp)f) - (i ^ j) 
= dif''Vf^9^{VaZp)-{i^j) 

= dirdjf\Vb{VaZp) - (a o 6)) 

= dif''djf''Ry{da,d,,)Zp 

= Ry{f*di, f*dj)Zp. 

□ 

In the case of pulling back a tangent bundle, there is another important property: 

Proposition 2.5. If V is a symmetric connection on TN, then the pull-back connection •'^ V on 
f*TN is symmetric, in the sense that for any U,V E T{TM), 

^Vf/(/,y) - fVv{f*U) = /.([[/, y]). 
Proof. Choose local coordinates for M near p, and j/" for near f{p), and write U = 

^Vu{f*V) -{U^V) = ^Vuiv^d^rda) -{U^V) 

= {u'diV^ - v%W)djf''da + u'v^didjr - djdif'')da 



□ 
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2.1.3. Subbundles. A subbundle K of sl vector bundle E over M is a vector bundle K over 
M with an injective vector bundle homomorphism lk K ^ E covering the identity map on M. 
We consider complementary sub-bundles K and L, so that E^ = lk{Kx) ® ll{Lx), and denote 
by ttk and ttl the corresponding projections onto K and L (so ttk ° lk = Id^-, 7r/,oi/^ = Id/,, 
T^K ° t'L = T^L ° I'K = 0, and i-i^- o + ^L^T^L = Id^;). If V is a connection on E, we define a 
connection V on and a tensor G r(T*M ® K* ® L) (the second fundamental form of K) 
by 

(2.5) V„C = vri^(V„(t/^e)); h^{u,0 = ^L{yu{iKi))] 
so that 

(2.6) V„(ii^O = iKi^uO + iLih^'iu, e)) 



K 



for any G TM and ,f G r(Er) C r(£'). The curvature of V is related to the second 
fundamental form and the curvature of V via the Gauss equation: 

(2.7) R'^'iu, v)^ = ttk{Rw(.u, v){tKO) + h''{^^, h^'iv, 0) - h^v, h^'iu, 0) 

for all u,v G TxM and ^ G ^{K). The other important identity relating the second fundamental 

form to the curvature is the Codazzi identity, which states: 

(2.8) 

T,L{Rv{v,u)iKi) = Vu{h''{v,0) - V,(/i^(^x,0) - h''{u,V,0 + h''iv,VuO - h''{[u,v],0- 

If we are supplied with an arbitrary symmetric connection on TM, then we can make sense of the 
covariant derivative Vh^ of the second fundamental form , and the Codazzi identity becomes 

(2.9) Vuh'^iv, e) - V.h^'iu, e) = 7:l{Rv{v, u){iKO)- 

An important case is where K and L are orthogonal with respect to a metric g on E compatible 

K 

with V. Then V is compatible with the induced metric , and and are related by 

(2.10) /(/i^(n,O,r?)+5'^(e,/i''(n,?7)) = 
foralU G r{K) and G r(L). 



2.2. The tangent and normal bundles of a time-dependent immersion 

The machinery introduced above is familiar in the following setting: If F : — )• A^"+'= 
is an immersion, then F^, : TM — )- F*TN defines the tangent sub-bundle of F*TN, and its 
orthogonal complement is the normal bundle NM = F^{TM)^. If ^ is a metric on TN with 
Levi-Civita connection V, then the metric g'^^^ is the induced metric on M, and V^*^ is its 
Levi-Civita connection, while h^^^ G V{T*M ® T*M (gi NM) is the second fundamental form, 
and /i^^^ is minus the Weingarten map. The Gauss identities (12.71 ) for TM are the usual Gauss 
equations for a submanifold, while those for NM are usually called the Ricci identities. The 
Codazzi identities for the two are equivalent to each other. 

Denote by vr the orthogonal projection from F*TN onto TM, by vr the orthogonal projection 
onto NM, and by l the inclusion of NM in F*TM. For u,v £ TM, equation (1231 ) is exacdy the 
usual Gauss relation: 

^VuF*v = F*(V„f) + ih{u,v). 
whilst for ^ G NM we recover the usual Weingarten relation: 

^VuF.v = L{VuC)-F,iW{u,0)- 

At the moment, the right hand side of both of these relations is really just notation expressing the 
fact that we have the decomposition F*TN = F^TM © iNM into orthogonal sub-bundles. We 
want to show, as the notation suggests, that the tangential component is the induced Levi-Civita 
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connection on M, and that /i is a symmetric bilinear form. Let a and j3 be functions on A^, then 
since is a connection and a and (3 restrict smoothly to functions on M, 



F 



VuPF^v = iuP)F,v + P^VuF^v. 



Therefore, 



and 



^VauF^:V = F^.{y auF^:v) + Lh{au,v) 
= aF^,(yuF^v) + aLh{au,v) 



^Vu^F^v = {uP)F^v + ^{F^VuF^v) + ih{u,v)) 
= F^{Vu(3F^v) + Lh{u,/3v). 
After projecting these equations onto TM and NM we get 

VauF.v = aVuF.v, Vu^F, = {u^)F,v + ^VuF.v 
h{au,v) = ah{u,v), h{u, f3v) = /3h{u,v). 

This shows that V := vro^VoF* is indeed a connection on M, and that h is bilinear. Furthermore, 
since is torsion-free and using Proposition I2.5l we have 

= ^VuF.v + ^VuF.v - [F^u, F^v] 
= F^{VuF^v) — F^,{VyF^,u) — F^{[u, v]) + i/i(n, v) — ih{v, u), 

which shows V is also torsion-free and h is symmetric. Finally, since is metric-compatible, 

for u,v,w e T{TM), 

Vu,{g{u,v)) = ^Vu,{g{F^u,F^v)) 

= g{^VwF^u,F^v) + {u v) 

= g{F^(yi,u),F^v) + {u^v) 

= g{VwU, v) + g{u, V^t>), 

thus by uniqueness of the Levi-Civita connection, the induced connection V is the Levi-Civita 
connection on M. Similar calculations show that the Weingarten map is bilinear in both its ar- 

guments, and that V := vro^Voiisa metric compatible connection on the normal bundle. 
Differentiating g{F^:U, i^) = shows that g{h{u, v),^) = g{{u, ^),v). In local coordinates {x*} 
for M near p and {y'^} for N near f{p) the Gauss-Weingarten relations become 

d'^F" dF" ^„ OF'' dF^ 



dx'-dxi dx^ '^^ dxi dx^ 

where rf^ are the Christoffel symbols of the submanifold, F"^ the Christoffel symbols of the 

ambient space, and C^^ the normal connection foms. The Christoffel symbols of the ambient 
space are obviously zero if the background is flat, and the normal connection forms are zero if M 
is a hypersurface. 

In this thesis we want to apply the same machinery in a setting adapted to time-dependent 
immersions. If / is a real interval, then the tangent space T(S x /) splits into a direct product 
© Rdt, where J{ = {u G r(S x /) : dt{u) = 0} is the 'spatial' tangent bundle. We consider 
a smooth map F : S" x / — A^"+'= which is a time-dependent immersion, i.e. for each t G I, 
F{.,t) : S — is an immersion. Then F*TN is a vector bundle over Ex/, which we can equip 
with the restriction metric gp and puUback connection coming from a Riemannian metric g on 
N and its Levi-Civita connection V. The map : J{ ^ F*TN defines a sub-bundle of F*TN 
of rank n. The orthogonal complement of F^,(;K) in F*TN is a vector bundle of rank k which 
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we denote by N and refer to as the (spacetime) normal bundle. We denote by vr the orthogonal 
projection from F*TN onto IK, and by tt the orthogonal projection onto 'N, and by l the inclusion 
of N in F*TN. The restrictions of these bundles to each time t are the usual tangent and normal 
bundles of the immersion Ft. 

The construction of the previous section gives a metric g{u, v) = g{F^u, F^,v) and a connec- 
tion V := vr o o on the bundle IK over S x /, which agrees with the Levi-Civita connection 
of g for each fixed t. We denote by g the metiic induced on given by ^(,^, r]) = g{L^, irf). The 

construction also gives a connection V := tt o o on !N". We denote by /i G r(IK* ® IK* ® IN") 
the restriction of /i^ = vr o ^ V o to IK in the first argument. Proposition 12.5 [ implies that /i is a 
symmetric bilinear form on IK with values in INf. The remaining components of /i^ are given by 

h^{dt,v) = ^{^VtF,v) 

= TT{^y,F,dt + F,{[duv]) 

(2.11) = V^{TtF^dt) + h{v,TrF^dt) 

where we used Proposition 12.51 Henceforward we restrict to normal variations (with irF^dt = 0), 
since this is the situation for the mean curvature flow. We also define W G r(IK* N (g) IK) by 
W(tt,0 = -h^{u,i) = --iri^VuiO for any u e r(IK) and ^ e r(N) (we refer to this as the 
Weingarten map). The Weingarten relation (I2.10l i gives two identities: 

(2.12) UHu,v),0 =9{v,W{u,0y, 

(2.13) g{h^{dt,^),v) = -^{V^.dt,^) 

where the latter identity used (12.111 ). The Gauss and Codazzi identities for IK and IN give the 
following identities for the second fundamental form: First, if u and v are in IK, then the Gauss 
equation (12.71) for IK amounts to the usual Gauss equation at the fixed time, i.e. 

(2.14a) R{u,v)w = W{v,h{u,w)) - W{u,h{v,w)) + Tr{R{F^u, F^v)F^w) 

(2.14b) R{u, V, w, z) = g{h{u, w), h{v, z)) — g{h{v, w), h{u, z)) + F*R{u, v, w, z). 

If u = dt but V ^ Ji, then we find: 

(2.15) R{dt,v,w,z) = g{Vu,7TF^dt,h{v,z)) - g{V zirF^dt, h{v,w)) + F* R{dt,v,w, z). 

X 

The Gauss equation for the curvature RofJ4 also splits into two parts: If u and v are spatial 
these are simply the Ricci identities for the submanifold at a fixed time: 

(2.16) R{u, v)C = h{v, W{u, 0) - h{u, W(t;, i)) + tt{R{F,u, F,v)iiO); 
while if u = dt and f G IK, then we have the identity 

X _ x^x x^x 

(2.17) R{dt,v,C,r]) = R{F^dt, F^v, l^, ltj) - c/(V-w(,,,o7rF*5t, ??) + 5(Vw(t.,,,)vrF*9i, ^). 
Finally, the Codazzi identities resolve into the tangential Codazzi identities, given by 

(2.18) Vuh{v,w) -Vyh{u,w) = Tt{R{F^v, F^u)F^w) 

for all u,v,w G r(IK), and the 'timelike' part, where u = dt and v,w ^ r(IK): 

(2.19) 7r{R{F^v,F^dt)F^w) = VaM^^w) - V^V^T^F^dt) - h{w,W{v,7rF^dt)). 

X 

Note that here Vh G r(r*(S x /) IK* (g) IK* ig) N) is defined using the connections V and V as 

in Equation (12.31 ). that is Va^/i(u, v) = Vg^{h{u, v)) — h{Vd^u, v) — h{u, Vg^v). 

We remark that by construction we have Vg = and = 0. In contrast to the situation 
in other work on evolving hypersurfaces, we have Vg^g = 0. That is, the connections we have 
constructed automatically build in the so-called 'Uhlenbeck trick' IIHa31 Section 2]. 
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Proposition 2.6. The tensors G r(J{* (g) F*TN), i G r(3^* F*TN), tt eT(F*TN ^ 'K) 
and Tt G T{F*TN (g> J{) satisfy 

(2.20) {VuF,){V) = ih{U,V) 

(2.21) {Vuim = -F*^{U,0 

(2.22) (Vc/vr)(X) = W{U, ttX) 

(2.23) {VuTt){X) = -h{U,TTX) 

for all U,V £ T{J{), ^ £ r(N) and X £ T{F*TN). 

Proof. These follow from our construction and Equation (I2.3l i: For the first we have (since 
F^, is a F*TA^-valued tensor acting on J{) 

(yuF,){V) = ^Vu{F,V) - F.iVuV) = F*{VuV) + ihiU,V) - F.iVuV) = ihiU,V), 

where we used the definitions of h and V. The second identity is similar. For the third we have: 

{Vu7T)iX) = Vu{7^X) - TrfVuX) 

= Vu{7tX) - tt{^Vu{F^ttX + lttX)) 
= Vt/(vrX) - Vc/(vrX) + W{U, 7tX) 
= W{U,7rX). 

The fourth identity is similar to the third. □ 

We illustrate the application of the above identities in the proof of Simons' identity, which 
amounts to the statement that the second derivatives of the second fundamental form are totally 
symmetric, up to corrections involving second fundamental form and the curvature of A^: 

Proposition 2.7. 

zh{u,v) — V yV vh{w , z) = h{v,W{u,h{w, z))) — h{z,W{w,h{u,v))) — h{u,'W{w,h{v, z)) 

+ h{w, W{u, h{v, z))) + h{z, W{u, h{w, v))) — h{v, W{w, h{u, i 

— h{u, t:R{F^v, F^,w)F^,z) — h{w, TrR{F^,u, F^,z)F^,v) 

— h{z, t:R{F^u, F^w)F^v) — h{v, TrR{F^,u, F^,w)F^,z) 
+ TiR{ih{u, v), F^w)FtfZ — 7rR{Lh{w, z), F^u)F^v 

+ 7rR{F^,u, F^w)ih{v, z) + 7rR{F^,v, F^z)ih{u^ w) 
+ TtR{Ft,v, F^:w)ih{u, z) + TtR{Ft,u, F^:z)ih{v, w) 
+ TtVptuRiF^v, F^w)F^.z — ttV F,wR{F*z, F^u)F^.v. 

Proof. Since the equation is tensorial, it suffices to work with u,v,w, z G r(J{) for which 
Vn = 0, etc, at a given point. Computing at that point we find 

zh{u, v) = Vwiyuh{z, v) + 7rR{F^,u, F^,z)F^,v) 

= ^vS^wh{v, z) + {R{u, w)h){v, z) + V^(7r^(F*n, F^,z)F^,v) 

= V^yh{w,z) +7tR{F^v,F^w)F^z) + {R{u,w)h){v, z) +V J^R{F^u,F^z)F^v) 

= V^^h{w,z) + {R{u,w)h){v,z)+V^R{F^u,F^z)F^v)+V.^R{F^v,F^w)F^z) 

where we used the Codazzi identity in the first and third lines, and the definition of curvature in 
the second. Since /i is a N-valued tensor with arguments in 'K, the second term may be computed 
using the identity (12.41 ) to give 

± 

{R{u, w)h){v, z) = R{u, w){h{v, z)) — h{R{u, w)v, z) — h{v, R{u, 'w)z). 
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This in turn can be expanded using the Gauss identity (I2.14al ) for R and the Ricci identity (12.161) 
for R. In the third term (and similarly the fourth) we apply the identity (12.31 ) to vr: 

V^{ttR{F^u,F^z)F^v) = V^Tt{R{F^u,F^z)F^v) + 7r{^V^{R{F^u,F^z)F^v)). 

In the first term here we apply the identity (I2.23l l. In the second we can expand further as follows: 

^V^{RiF^u,F^z)F^v) = {^V^,R){F,u,F,z)F^v + R{{Vn,F^)u,F,z)F,v 

+ R{F^u, V^F,{z))F^v + R{F^u, F,z)(V„F,(t>)). 

In the terms involving VF* we apply (12.201 ). and we also observe that ^V^R = VF^.wR^y the 
definition of the connection ^V. Substituting these identities gives the required result. □ 

In subsequent computations we often work in a local orthonormal frame {ej} for the spatial 
tangent bundle Ji, and a local orthonormal frame {i^a} for the normal bundle We use greek 
indices for the normal bundle, and latin ones for the tangent bundle. When working in such 
orthonormal frames we sum over repeated indices whether raised or lowered. For example the 
mean curvature vector H G r([N') may be written in the various forms 

H tVg h Q ^ hij hi 9 ^ ^ij hiify U(y . 

Similarly, we write = g''^ g^^ g^^h.ij'^hki^ = hijahija- The Weingarten relation (I2.12l i be- 
comes 

while the Gauss equation (12. Hal ) becomes 

where we denote Rijki = R-iF^^ei, F^Cj, F^^ek, F^ei). The Ricci equations (12.161) give 

X 

^ijajB — hipahjpp hjpahipp -\- Rijaj3) 

where Rijap = R{F^,ei,F^ej,iVa, i^i^p), and the Codazzi identity (12.181) gives 

^ ihjk ^ jhik — Rjika^a- 

In this notation the identity from Proposition I2.7l takes the following form: 

^ ihij — ^ jhi^i + hi^io^hipcthjp hijQ.hiip(yhip 

~l~ hjiQ.hipQi^hkp + hji^^yhipahip hii(yh}^pQ^hjp hjiahj^pahip 

~l~ hj^iQ,Ri(yj jjV hijo^Rf^Q^if^iyf^ + Rkjiphip + RkUphjp Rnjphkp Rikjphip 

+ '^iRjkip^p — ^kRlij(5^(5- 

Particularly useful is the equation obtained by taking a trace of the above identity over k and I: 

^hij — \/ i\/ j F[ -\- F[ ' hiphpj hij • hpqhpq -\- Ihjq • hiphpq hiq ' hqphpj ^jq ' hqphpi 
+ HfyRioijpVj^ hijaRkakfi^fi ~\~ Rkjkphpi ~\~ Rkikphpj IR^pjqhpq 
(2.24) + 2hjpaRipal3l^[S + 2hipaRjpal3^li + ^iRjkkfB^P — ^ kRkij fi^^ P ■ 

Here the dots represent inner products in !N. 
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2.3. The method of moving frames 

The Gauss, Ricci and Codazzi equations can also be quite nicely derived using Cartan's 
method of moving frames. We shall only need to use this machinery once, and then only briefly, 
in Chapter 6, but it is nonetheless instructive to see how this can be done. We work in the setting 
of a fixed immersion F : M" — )■ A^"+'^, and let {ca : 1 < a < n + fc} be an adapted local frame 
for N, so that {cj : < i < n} are tangent to M, and {cq, :n + l<a<n + fc}are normal to 
M. In the Cartan formalism, the Levi-Civita connection on N is given by the structure equations 

duj" = -u:l A 

The first structure equation determines a torsion-free connection, and the second guarantees the 
connection is metric-compatible. The curvature of the connection is given by the structure equation 

(2.25) dw^ = -w;? + 

We now restrict the indices to M, so = 0. From the above structure equations we obtain 

(2.26) du:' = -co] A J 

(2.27) +ojj = 

(2.28) doj} = -uji A wj^ + n) 

which determines the Levi-Civita connection on M. Furthermore, since t^;' = in the normal 
bundle, we also obtain 

dw" = -o;^ A oo^ 
+ = 

which determines a metric compatible connection on the normal bundle. Let us show how to 
derive only the Codazzi equation. Since = duo'^ = —uf A lo^ , by Cartan's Lemma the uif can 
be expressed as a linear combination of the lo^: cof = hij°^io^ , and also hij" = hji"'. In equation 
(12.251) we restrict a to a and & to i to get 

(2.29) dwf + uj'f A oj{ + A t^f = Clf . 

Exterior differentiation of ujf = /ijj^u;' gives dhu"' Aco'' — hij°'coj A w', and putting this together 
with ( 12.291 ) we obtain 

We define 

hijk^ui^ = dhij°' — hij'^ojl — hii°'u}j + hij^co'^, 
which is just the first covariant derivative of h, and assuming a flat background, we recover the 
usual Codazzi equation: /ijjfc" = 



CHAPTER 3 



Short-time existence theory 

The mean curvature flow equation determines a weakly parabolic quasilinear system of sec- 
ond order It is now well-known that many geometrically-defined partial differential equations 
possess zeroes in their principal symbol because of some kind of geometric invariance displayed 
by the equations. In order to assert short-time existence to the mean curvature flow we use the 
well-known 'DeTurck trick' to first solve a related strongly parabolic equation, and we then re- 
cover a solution to the mean curvature flow from this related solution. The DeTurck trick was first 
invented to solve the Ricci flow, however the method applies to many other geometric flows. In 
his lecture notes fHa41, Hamilton shows how the DeTurck trick can be applied to the Ricci, mean 
curvature and Yang-Mills flows. We have also seized this opportunity to fill in a few details in 
the proof of short-time existence for fully nonlinear parabolic systems of even order defined on a 
manifold. In IIHall Hamilton gives a proof of local existence for the harmonic map heat flow (a 
strongly parabolic quasilinear system) using Sobolev spaces and the inverse function theorem, and 
we were inspired to adapt his proof to fully nonlinear systems in the Holder space setting. Towards 
the end of this task Tobias Lamm pointed out to us the he proved the short-time existence of fully 
nonlinear operators in Euclidean space using Schauder estimates in his Diploma Thesis [L]. Given 
that we started reconstructing this theory on our own, we have still decided to include this chapter 
Since we are not claiming anything essentially new in this chapter, we have freely borrowed from 
Lamm's thesis to improve our own exposition. In particular, we now use Simon's method of scal- 
ing to derive the Schauder estimates for parabolic systems, as opposed to Trudinger's method of 
mollification which we had originally used. We still show how Trudinger's method can be com- 
bined with the mean value property of subsolutions to the heat equation to provide a remarkably 
simple proof of the Schauder estimates in the case of single equations. We emphasise the global 
aspects of solving the problem more than Lamm, and we work in the setting of parabolic systems 
defined in sections of a vector bundle over a closed manifold. 

The strategy for proving such an existence theorem is well-known: one begins with a solution 
to the heat equation and then uses the method of continuity and the Schauder estimates to prove 
existence for general linear operators. The short-time nonlinear existence result then follows by 
linearising the nonlinear operator and applying the inverse function theorem. In reconstructing 
the and linear theory, our main reference has been the Chinese text [Gu|. As we have men- 
tioned above we use Trudinger's method of mollification to derive the interior Schauder estimate 
for second order parabolic equations, and we simply cite Lamm's thesis for the derivation of the 
Schauder estimates for even order parabolic systems in Euclidean space. The application of the in- 
verse function theorem to yield the nonlinear existence result was inspired by HHali Our method 
is different to Lamm's, however in showing the solution is an appropriate Holder space with ex- 
ponent /?, for /3 < a where a the Holder exponent of the initial data, we have benefited from 
D . The application of the DeTurck trick to the mean curvature flow first appears in [Ha4|, and 
we have simply expanded on these notes of Hamilton's, adding in a few calculations. Combining 
the harmonic map heat flow with the mean curvature flow to give a simple proof of uniqueness of 
the mean curvature flow is also due to Hamilton; the equivalent result for the Ricci flow first ap- 
peared in pHa6|. We adapt the Ricci flow result to the mean curvature flow, following the detailed 
expositions given in l,HA.iCLNi . 



15 



16 



3. SHORT-TIME EXISTENCE THEORY 



3.1. Short-time existence for fully nonlinear parabolic systems of even order 

In this section we give a proof of short-time existence for fully nonlinear parabolic systems of 
even order. The nonlinear existence result that is our ultimate goal is attained by an apphcation of 
the classical inverse function theorem in Banach spaces, and is in fact quite short once we have 
all the linear theory in place. The linear theory plays an essential role in the nonlinear theory and 
most of the following is devoted to establishing the linear theory. 

The setting for our study of systems of partial differential equations defined on a manifold is 
slightly different to that of the more familiar Euclidean case. Here we are interested in differential 
operators that act on sections of a vector bundle over a manifold, and not simply functions defined 
on some domain of Euclidean space. So that the reader can accustom to this setting, let us first 
consider linear systems of second order Let E and F be two vector bundles over M. Let the 
indices a,b,c,... range from 1 to N, and the indices k,... range from 1 to n. Suppose {ca} 
is a local frame for the bundle E over a coordinate neighbourhood of M with local coordinates 
{x*}. It may be helpful to keep in mind the specific example of the mean curvature flow, in which 
case the section we are interested in the the position vector of the submanifold, N is dimension of 
the background space and n the dimension of the submanifold. A linear differential operator of 
second order is a map L : V{E) r(F) which in any coordinate chart is of the form 

L{U) := ^U- - Al'\x,t)d^djU' - Bf{x,t)dkU' - Ct{x,t)U\ 

where A £ T{Sym^{T*M) (g) E (S) E*), B G TiJ*M (g) E (g) E*) and C £ T{E x E*). The 
notion of parabolicity is defined in terms of the principal symbol of the differential operator The 
principal symbol a of the above linear operator L in direction ^ G T{TM) is the vector bundle 
homomorphism 

The system is said to be strongly (weakly) parabolic if the eigenvalues of the principal symbol 
are positive (non-negative). The principal symbol encodes algebraically the analytic properties of 
the leading term of the differential operator. A linear differential operator of order 2m is a map 
L : T{E) r(F) which in any local coordinate chart is of the form 

(3.1) ^u'^ + i-ir Yl ^'^i^' 

I<2m 

or in full 

^C/a + {-ir{Al'''--'-{x,t)di^ ■ ■ ■ di.^U'^ + ■■■ + Bt\x,t)dkU'' + C^{x,t)U'). 

The principal symbol is defined in a similar manner as before. Let us now move on to fully 
nonlinear differential operators. A fully nonlinear differential operator of order 2m is a map 
L : T{E) r(F) which in any coordinate chart is of the form 

(3.2) ttU"- F^ix, t, U'" , dU^ d^'^U''^'^), 

at 

or equivalently, is globally of the form 

(3.3) ^U" - F^ix, t, U''-' , V[/^^ , . . . , V^^f/''^'"). 

We will often drop the indices relating to the section, and simply write the above equation as 

(3.4) -U-F{x,t,U,VU,...,V^"'U). 

We have chosen to categorise the equations (linear, quasiUnear, etc) in terms of their form in a 
local chart. This is important for operators that are not fully nonlinear, because if the connection is 
varying in time, the initial global appearance of the equation can be betraying. As we shall see in 
the next section, this is indeed the case for the mean curvature flow. This difference is immaterial 
for fully nonlinear operators, and so nothing is lost by saying such an operator looks globally of 



3.1. SHORT-TIME EXISTENCE FOR FULLY NONLINEAR PARABOLIC SYSTEMS 



17 



the fomi (I3.4I ). We mention that as an alternative to passing to local coordinate descriptions, we 
could define the equations globally with respect to a fixed connection, which then reveals the true 
nature of the equations. 

The notion of parabolicity for a nonlinear operator is defined in terms of its linearised operator 
The linearisation of a nonlinear operator L about some fixed function Uq in the direction V is the 
linear operator given by 

dL[Uo\{V) = ^L([/o + sy)[^^ 

= d^v - F'''-''^^-{x, t, Uo,VUo, . . . , v2™;7o)Vi, • • • V^,„y 
F'^ix, t, Uo, VC/o, . . . , V^'^Uo)VkV - F{x, t, C/q, VC/q, . . . , V^^C/q)!/. 

A fully nonlinear operator is said to be parabolic if its linearisation evaluated at the initial time is 
parabolic. The principal symbol of the linearised operator evaluated at the initial time is 

a{dL[Uo]m = 0, U^' ,VU'q\ . . . ,V^^U','-)^,, • • • e..„e„ e*''- , 

and thus the nonlinear operator is strongly (weakly) parabolic if a{dF[Uo]){^) > (>)0. Closely 
related to the notion of parabolicity is the Legendre-Hadamard condition. The linear operator 
(13.321 ) is said to satisfy the Legendre-Hadamard condition if there exists a positive constant A > 
such that coefficient of the leading term satisfies 

for all r/ e r(£') and r/* G T{E*). We can now state the local existence theorem we wish to 
prove. 

Main Theorem 5. Let E x (0, uj) be a vector bundle over M x (0, a;), where M is a smooth closed 
manifold, and let U be a section T[E x (0, w). Consider the following initial value problem: 



(3.5) 



\P{U) := dtU - F{x,t,U,VU, . . . ,V^"'U) = in E X (0,w), 
\u{M,0) = Uo, 

with Uq G C^'"'^'"(£'aj)- The linearised operator of P at Uq in the direction V is then given by 
dP[Uo]V = dtV + i-ir Yl A\x,t,Uo,VUo,...,V^"'Uo)ViV. 



\I\<2m 

Suppose that the following conditions are satsified: 

tailjl---imjm 



1) The leading coefficient ' satisfies the symmetry condition ^f-*!-?! ' *™-?™ 



2) The leading coefficient satisfies the Legendre-Hadamard condition with constant A 

3) There exists a uniform constant A < oo such that ^\n<^2m 

4) F is a continuous function of all its arguments 

Then there exists a unique solution U G C"^"^'^'^ (E^^), where (3 < a, for some short time > 
to the above initial value problem. Furthermore, if Uq and all the coefficients of the linearised 
operator are smooth, this solution is smooth. 

3.1.1. Function spaces and preliminary results. We introduce some more of our notation 
and the necessary function spaces. We denote differentiation in space by di or dx, and differen- 
tiation in time by dt. As an example of this notation, the spatial gradient '}2!^=i X]f=il^«^"P 
l^i^up. The kih order derivative is denoted by a raised index: d^. We set 5'^'' equal to + d\\ 
note the lack of a subscript for this combined derivative. Usually this combined derivative will 
be used when referring to the 2mth order derivative in space and the first derivative in time: 
Q2m,i ._ Q2m ^ ^j^^^ need to work in both parabolic Holder and Sobolev spaces, ini- 

tially on R'^^^ and later on a closed manifold. We reserve Q. for an open domain contained in W^, 
and M for a closed manifold. We shall work with the parabolic domains P := Q. x (0,a;) and 



18 



3. SHORT-TIME EXISTENCE THEORY 



:= M X (0,w). For any two points X = {x,t), Y = {y,s) G P, the parabolic distance 
between them is given by 

d{X,Y) = max{|x — y\, \t — s\^}. 

We shall also work with the parabolic domain Pg := {X = {x,t) e P : dist{x,Cl) > 6} and 
the backwards parabolic cylinders Qr[Xq) := {X G x R : (i(X,Xo) < < to} = 
Br(xq) X [to - i?^™, to)}. Let u : P ^ . For a G (0, 1), the Holder semi-norm is given by 

[u]a;p:= sup ■ 

The Holder norms are given by 

2m 

k|2m,l;P := ^Id^ulo; P + \dtu\o;P 

k=0 

\u\2m,l,a;P ■= |«|2m,l;P + [9^"'"^«]a; P- 



The set of functions 

{ueC^'^P) : [u]2m,l,a;P<00} 

endowed with the norm |'u|2m,i,a; p is called a Holder space. Written out in full the norm is 

2m 

\u\2m,l,a;P ■= ^\dx'^\o; P + \dtu\o- P + [dl'^uja-P + [dtu]a-p. 
k=0 

These Holder spaces are Banach spaces. Next we define the analogous spaces on a closed mani- 
fold. Let E he a vector bundle over a closed manifold M, and let E and TM be equipped with 
metrics g and connections V. Let dg{X, Y) be the geodesic distance on M measured by g, and let 
ig be the injectivity radius of the manifold M. For any two points X = {x, t), Y = {y, s) G M^, 
the parabolic distance between them is given by 

d{X, Y) = max{dg{x, y),\t — s\^}. 

The definition of the Holder space of functions u on mimics that of functions on W^'^^. For 
a G (0, 1), we define the semi-norms 

2m 

\u\2m,l;M^ ■= ^\^xU\o;M^ + \dtu\o. 
k=0 

|^^|2m,l,a;Ma, := N2m,l;M,^ + SUp . 

Here we use the notation V^"*'^ := V^"* -I- df. To define Holder spaces of sections of E we need 
to be a little more careful: The points X and Y live in different vector spaces above M, and so 
parallel translation is needed to indentify the spaces in order to perform the subtraction. As we are 
working on a closed manifold, geodesies always exists between any two points, however beyond 
the injectivity radius the geodesies may not be unique. Let 'J'y,x denote the parallel translation 
along a geodesic from Y to X. For U G T{E^), we define the norms 

2m 

\U\2m,l;E^ := ^[Vl'^U {X)\o.,E^ + \dtU {X)\o; 
k 

irn |V2™'lC/(X) -TyxV2'"'l[7(F)| 
\U \2m,l,a; E^ ■■= \U \2m,l; E^ + SUp . 

dg{x,y)<ig 

The norm |V^'"'^C/(X) - 'J'y,x'^^'^'^U(Y)\ is measured by the bundle metric g, but for conve- 
nience we shall supress this dependence in our notation. We mention in passing that these Holder 



3.1. SHORT-TIME EXISTENCE FOR FULLY NONLINEAR PARABOLIC SYSTEMS 



19 



spaces are well-defined on closed manifolds, since on a closed manifold all metrics are equivalent 
and the injectivity radius is always positive. Both these conditions fail to be true on arbitrary com- 
plete manifolds. The remaining definitions are entirely analogous to the Euclidean case and we 
repeat them to avoid confusion. The Holder norms are given by 

2m 

\U\2n.,l;E^ ■-= Y.\^'.^\o; + \dtUW, 
k=0 

\U\2m,l,a;E^ := \U\2m,l-E^ + [^^""''UU; E^- 



The set of tensor fields 

{[/ G C72-'l(£;^) : [C/]2m,l,a;i?. <oo} 

endowed with the norm |f/|2m,i,a; E^j is again called a Holder space and it is easily verified that it 
too is a Banach space. Written out in full the norm is 

2m 

\U\2m,l,a;E^ := Y.\'^'a:U\o; E^ + \dtUW,E^ + [Vl^^U^E^ + mUE^- 
k=0 

Next we introduce the anisotropic Sobolev spaces we wish to work in. The set 

{u : did{u e L'^{P),i + 2mj < 2m} 

endowed with the norm 

1/2 




\did{u\'^dxdt 



^ i+2mj<2m 



is the Sobolev space denoted by W2^'^ (P). On a manifold the norm is given by 
Mw?-'He.)--= ill E I'^M^l'dV.dt' 

^^•^ i+2mj<2m 

These spaces are also Banach spaces. If we interchange the order of the covariant and time deriva- 
tives in our definition we obtain an equivalent norm. Note that each time derivative counts for 2m 
space derivatives, and that we have again suppressed the dependence on the bundle metric. We 
shall also need the following spaces, in which the highest order spacial derivatives is of order m: 

Mwr'\E^) ■■=([[ {Y.\'^S\' + \dtu\')dv,dt) ' . 

We also define the space 

V{E^) = {[/ G wI^\E^) : Vdt G h'iE^)], 
* 1 1 

and note that V{Ei^) is dense in W2 ' (E^). 

Let C'^{^t) be the set of all smooth functions that vanish near the spatial boundary {(x, t) : 

• 

X G dn,t G (0,0^)} of P, and let C°°{^aj) be the set of all smooth functions that vanish near 
the parabolic boundary {{x,t) : x G dO,,t G (0,T)} U {{x,t) : x G H, t = 0} of P. Denote 
by t^2™'\P) the closure of C'°°(IIt) in W^'^'^iP), and by W^'^'^iP) the closure of ^^(II) in 
W2"^'^{P). We similarly define the spaces W^'^{P) and W^'^{P). 

To close out this section we recall some important results that are used in the following. The 
first result is a well-known covering lemma that allows us to patch local Euclidean estimates 
together to give global estimate on the manifold. For a proof we refer the reader to WHall Corallary 
4.12] and IHel. 



20 



3. SHORT-TIME EXISTENCE THEORY 



Lemma 3.1. Let (M, g) be a Riemannian manifold, p £ M and tq G (0, ig (p) / 4). Suppose that 
for each q > there exist constants Aq such that |V^-Rm| < Aq in Bp{rQ). Then in normal 
coordinates {x*} on Bp[rQ) there exist constants Cq = Cq{n,ig, Aq, . . . Aq) such that for each q 
the estimates 



^^ij < Qij < 2(^ij and 



dxP 



<Cq 



hold in Bp{m\n{Ai / \/Aq, fo})- 



The second result is known as Garding's inequality. Garding's inequality on M is a well- 
known result. The inequality also holds on a closed manifold, where one uses patching arguments 
similar to those we shall use later on, to lift the Euclidean estimate on to the manifold. A proof of 
Garding's inequality in the Euclidean case can be found in many places; for example [Gi|. 

Lemma 3.2 (Garding's inequality). Suppose that A'^^-^ is a smooth section of a tensor bundle E 
over a smooth closed manifold M that satisfies the Legendre-Hadamard condition. Then there 
exist positive constants Aq and Ai such that the bilinear form defined by 

'B{U,V) := I Af-^V'pU^VjU^dVg 

JM 

satisfies the inequality 

'B{U,V)>Xo [ \V'^U\^ dVg - \i f \U\^dVg. 

JM JM 

We shall also require the Poincare inequality: 

Proposition 3.3 (Poincare inequality). Let M be a smooth, closed manifold. For any u G Wl 
there exists a positive constant C such that 

(3.6) f \u\^dVg<C f \V^u\^dVg. 

JM JM 

For a proof this proposition we refer the reader to [TTe? pg 40.]. Since the Poincare inequality 
holds at each tuneslice of we can integrate (13.61 ) in time to get 



l-up dVgdt<C / / iV^up dVg dt, 

Mu; J JMu; 

and then by the Kato inequality |V|Vn|| < |V^n| we also obtain 
(3.7) [[ \u\'^dVgdt<C [[ \V^u\'^dVgdt 

J J Mui J J Mu, 

for any q>l. 

3.1.2. Hilbert space theory. We commence our existence program by studying differential 
operators in divergence form. Consider the problem 



(3.8) 



UU- + (-1)1-^1 Eo<|/MJ|<„^ Vj(^^"(X)V,C/'') = F-{X), XeM^ 
\uiM,0) = Uo, 



If Uq is is sufficiently smooth then we can consider the problem for y := U — Uq, so without 
loss of generality we can assume Uq = 0. Ultimately we are interested in smooth solutions, so 
for us Uq will always be smooth and this transformation is always possible. Henceforth, we will 
usually assume Uq = 0. For simplicity, we assume that the connection does not depend on time, 
so we can commute time and space derivatives without introducing derivatives of the Christoffel 
symbols. We now want to introduce the notion of a weak solution to the above problem, and then 
recast the problem in terms of bilinear form on a Hilbert space. 
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Definition 3.4. A section U £ W^' {T{Ei^)) is called a weak solution of the initial value problem 



3.81 ) if for any ip G C°°{T{Ei^)), the equation 



f^tV+ Al^'^ViU^Vj^" + Y B^^-^ViU''Vj(pA dVgdt 

|/| = |,/|=m 0<|/|,|J|<m ^ 

/ / F V dVg dt 



^^■^ ^ |/| = |,/|=m 0<|/|,|J|_ 

(3.9) |/| + |J|<2m-l 



holds. 



Since (7°° (E^^ ) is dense in W^' {M^ ) , the test function can in fact be any function in W^' {E^ 
For ease of reading will again often drop the indices running over the section U and simply write 

Ut^+ Y A^-^VfUVjip+ B^'^ViUVjip\dVgdt= 1 1 FipdVgdt 

^ \I\ = \J\=m 0<|/|,|J|<m ^ JJM^ 

\I\ + \J\<2m-l 

We have the following two characterisations of weak solutions. 

Proposition 3.5. A section U G W^' {E^^) satisfies ( 13.91 ) if and only ifU satisfies 
(3.10) 



UtUt+ Y At'^^fU'Vjipl + Y B^'-'ViU'VjipAdVgdt 

^ \I\ = \J\=r~~ n^lrll7l^™ / 

F'^if^ dVg dt 



\I\ = \J\=m 0<|/|,|J|<m 

|/| + |J|<2m-l 



for any if e C°°[E^). 



Proof. Suppose that U G W^' {E^) satisfies ( |T9l ) for any (/? G C'^{E^). Because (/? is 
smooth, ipt is a valid test function, and so (13.101 ) holds. Conversely, suppose thatf/ G W^'\E^) 

o ^ 

satisfies (13.101) for any G C°°(Maj). Then since Lp is smooth, (/7(x,s)ds is a valid test 
function, and choosing the test function as such in (13.101 ) shows ( 13. 9b holds. □ 

Proposition 3.6. A section U G 14^2"*' (^i.^) satisfies (13.91) if and only ifU satisfies 
(3.11) 

// Y At''VTU'Vj^1 + Y Bt''ViU''Vj^i\e-''dVgdt 

■^■^^^'^ ^ \I\ = \J\=m 0<|/|,|J|<m ^ 



|/| = |J|=m 0<|7'|,|„ 

|i'| + |J|<2m-l 



[[ F-^te- 

J J M,., 



dVgdt 



for any ip G C°°{Ei^), where 6 is a positive constant. 

Proof. Suppose that U G W^''^{E^) satisfies ^3 for any ip G (7°°(£'^). Because 99 
is smooth, v^te^^* is a valid test function, and so (13.111) holds. Conversely, suppose that U G 
W^'\E^) satisfies (IXTTT) for any ip G C°°(Fa;)- Then since 99 is smooth, ip{x, t)e 6* /q 99(2;, s)e 
is a valid test function, and choosing the test function as such in (13.1 II) shows (13.101 ) holds, and 
thus □ 



We shall use the Lax-Milgram lemma to prove existence and uniqueness of a weak solution to 
problem (13.91 ). Our approach is similar to that of f.HPJ and llShl . where slightly different function 
spaces were used. 
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Theorem 3.7. Let H be a Hilbert space and V an inner product space continuously embedded in 
H. Let H : H xV ^M.be a bilinear form with the following properties: 

1) ForallU G H andW G V , there exists a constant C such that \'B{U,W)\ < C\\U\\h\\W\\v 

2) !B is coercive, namely, there exists a X > such that 23(VF, VF) > A||l^||y 

Then for any bounded, linear functional F{W) in H, there exists a U & H such that F(W) = 
"Bill, W)for each W £ V. Moreover, ifW is dense in H, then U is unique. 



For a careful proof of this result we recommend to the reader llShi pg 118.]. 

Theorem 3.8. If F, Uq £ L'^(E^), then the initial value problem (13.81 ) admits a weak unique 
solution U G W^'\e^). 

Proof. Let U G W^'^{E^), V G V{E^) and 6 be some constant greater than zero that will 
be fixed later on. Consider the bilinear form associated to the differential operator in problem (13. 8t 
(3.12) 

B{U,V):= ff (utVt+ ^fe^"^Vf[/VyF"+ B^^-WiU^VjvAe^'^'dVgdt. 

^^'^ \I\ = \J\=m 0<|/|,|J|<m ^ 

\I\ + \J\<2m-l 

We want to show that the bilinear form S satisfies the conditions of the Lax-Milgram Lemma. 
First, it's easy to see 

and so S is bounded. Next we show 23 is also coercive. For convenience, write S = Ii + I2, 
where Ii and I2 refer to the two summation terms of (I3.12l i. Focussing on /i, for V G V{Ei^) we 
have 



Ii= [[ V At^-^ViV''VjVt''e-^^dVgdt 



|/| = |J|=m 

>\jj ^^{A%''ViVVjV)e-'' - \\dtA\^ 1 1 ViVVjVe-''dVgdt 
^ \ H ^ {A%''ViVVjVe-'') + ^ [[ At ''ViVVjVe-'' dVg dt 



\\dtA\o jj^^ ViVVjVe-''dVgdt 



Upon integrating the first term on the right we find both terms are non-negative: the endpoint 
t = T from Garding's inequality and t = because V G V{E^), and we discard these terms. We 
are left with 

By choosing 6 sufficiently large the first term on the right can be made positive. Now we deal 
with 12- By using the Peter-Paul inequality on the terms of I2 they are either of the order |V"^yp 
multiplied by an e, or lower order terms divided by e. In the case of the former, they can again be 
absorbed by choosing 9 sufficiently large. In the case of all lower order terms, they can also be 
absorbed by using the Poincare inequality and then choosing 6 sufficiently large. After all such 
estimation we obtain 

s(y,y)><5||y||^p,i(^^) 

for some constant S > 0. This shows !B is coercive and we may now apply the Lax-Milgram 
Lemma. In Theorem l3.7[ choose W^'^{E^) as the space H, V{Ei^) as the space W, and F{V) = 
J f^j FVte^^^ dVg dt. By the Lax-Milgram Lemma, there exists a unique U G VF™'^r(£'tj) such 
that S(C/, V) = F{V) for all V G ViE^). Thus U is the unique weak solution to problem dM) 
by Propostion 13.61 □ 
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Next we discuss the regularity of the weak solution to problem (I3.8I ). We first need to first 
recall some basic facts about difference quotients. Difference quotient approximations to weak 
derivatives are a common tool in PDE and proofs of the following facts can be found in many 
texts, for example BGillGTl . Let u : Q ^ be function and ft' CC 0,. The difference quotient 
in direction Cfc is defined for all x £ Q' by 

u{x + hek,t) - u{x,t) 

Dk.h{x) := , 

n 

where < \h\ < dist{n', n) and k = 1, . . . , N. 

Proposition 3.9. 1) Suppose u G W2{ft) and \ < p < oo. Then for each ft' Gd ft the 
estimate 

(3.13) \\Dh,k\\L'\n') < C\\du\\L2(^n) 

holds for some constant C and allO < h < {l/2)dist{fl' , dfi). 
2) Suppose that u G LP{fl'), \ < p < oo, and that there exists a constant C such that 

\\Dh,ku\\LP{ni) < C 
holds for allO<h< (1/2) dist{n' , dU). Then du G L'^{VL') and 

\\du\\L2(jii) < C. 

Proposition 3.10 (interior regularity). Suppose that F G L'^{E^) and U G W^'^{E^) is a weak 
solution to problem (13.81) . Then U G VF^ '"'^(-E'a;) (^nd the estimate 



holds. 

Proof. We give the proof for the case m = I. The proof for systems of even order follows in 
a similar manner way, with small changes needed to incorporate the scaling of the system; in this 
regard see HHPl . Regularity is a local problem, so we derive the necessary regularity estimates on 
Euclidean space and then lift them to the manifold using patching argument. As our starting point 
we therefore work with the following definition of a weak solution 

(3.14) J J u^^p" + Al'^diu^dj^p'' + B^'^diU^djip''dxdt = J J f^p^dxdt, 

O ^ Q 

which holds for all ip G W^' (P). Rewrite this as 

u^if" + Al'^div^djif" = jj dxdt, 

where g"" = f"" — B'^^-' diu''djip°'. Choose if = —D_fi,k{'n'^ Dfi^ku"')X[o,s]- This is a valid choice as 
we have restricted h to be sufficiently small. With this choice of if equation (13.141 ) reads 

/ / u1{-D.f,,kirj^Df,^ku'')) + Al''diU%{-D_h,kiv^Dh,ku''))dxdt 

(3.15) -^-^ 

g''{-D_h^k{v^Dh,kun)dxdt. 
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We focus on the term involving the time derivative. Using the properties of difference quotients 
we have 

dt[ri'{Dh,ku''f) dxdt 



{rj\Dh,ku{x,trf) dx'~ 



in 

[ r]\Dh,kuix^sTfdx. 
In 

Now focus on the second term on the right of (I3.15I ). By Proposition 13.91 in order to achieve 
the desired spatial regularity it suffices to suitably bound the norm of Dh,kdu. Using various 
properties of difference quotients we estimate 

Al'^d,u''dj[-D.h^k{ri'Dh,ky')) dx dt 

Al'^diu'' - D^h.kdj{ifDh,kU^)dxdt 

p 

Dh,k{K''9iu')dj{7]^Dh,kU^) dx dt 

p 

= jj^ [aI'\x + hek,t)Dh,kdiu' + iDh,kAl''d^u'y,{ri^Dh,kun dxdt 

= jj^ {aI'\x + hek,t)Dh,kdiv!' + {Dh,kAl'' diu^) {r]^djDh,ku'' - 2rjdjr]Dh,ku'') dx dt 

= 11 r]'^Al'^{x + hek,t)diDh,kU^djDh,ku'' dxdt 

- lijdjT] j j Al'^x + hek, t)Dh^kdiu"-Dh,kU^ dx dt 

+ ll^{Dh,kAl'^)d:{v^djDh,ku'' - 2r]d,vDh,kun dx dt 

>rj^\o jj \Dh,kdu\^ dx dt - j j SI dxdt. 

In going to the last line we have used Garding's inequality and grouped the remaining terms into 
the term By using the properties of difference quotients and the Peter-Paul inequality, S'^ as 
well as the term involving g on the right hand side of equation (13.151 ) can both be estimated by the 
norm of du and /. Recombining this estimate on the spatial derivatives with the estimate on 
the time derivative gives 

sup f rj\Dh.ku''{x,s)fdx+[[ \Dh,kdu\'^dxdt<C{\\du\\L2^p) + \\f\\L2(P)). 
o<s<LuJn J J Ps 



From Proposition I3.9l it now follow that 

|a2n|2dxdt<C(||9n||i,(p) + ||/||i,(p)). 

Ps 

The estimate on the time derivative can be proved in a similar fashion using the difference operator 
in time. For the time derivative we obtain the estimate 



jj^ {dtufdxdt < C(||5u||i.(p) + ||/||i2(p)), 



'Ps 

and combining the space and time estimates complete the proof in the case m = 1. □ 
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The near-bottom boundary estimate can also derived in a similar fashion, and then the local 
interior and near-bottom estimates can be lifted to a closed manifold using a patching argument to 
give an estimate holding globally on E^. Higher regularity estimates can also be obtained using 
standard bootstrap arguments. As a simple consequence of the higher regularity estimates and 
the parabolic Sobolev inequality (on a manifold) we have the following existence theorem in the 
smooth category. 

Corollary 3.11. Suppose that U E W^'\E^) is a weak solution to problem (13.81 ). If F E 

C^{E^), then U £ C^{E^). 

3.1.3. Schauder theory. In this section we derive the interior Schauder estimates in Eu- 
clidean space, and then lift these local estimates to the vector bunlde E^^ to obtain a global 
Schauder estimate holding on the bundle. A number of methods can be used to derive the Schauder 
estimates; we shall present two of these. Trudinger's method of mollification offers a simple proof 
of the Schauder estimates for elliptic and parabolic equations of second order. This method ex- 
tends to systems of even order, and indeed we pursued this route in an early draft of this thesis. But 
perhaps an even easier and cleaner method of the deriving the Schauder estimates is Leon Simon's 
method of scaling HSimli We shall use Simon's method to derive the estimates for systems of 
even order. 

Trudinger's method was introduced in fV^ where he treated both equations and systems of 
elliptic type. For second order equations, the method is remarkably simple, and makes use of 
the solid mean value inequality. His method of mollification extends to systems of even order, 
where the application of the mean value inequality is replaced by an estimate and the Sobolev 
embedding theorem. Wang [Wa| has used Trudinger's method of mollification to derive Schauder 
estimates for second order parabolic equations, where the application of the solid mean value 
inequality was replaced by estimates coming differentiating the fundamental solution of the heat 
equation. Here we show how the mean value property of the heat equation can be used in exactly 
the same way as the solid mean value inequality to provide the desired estimates. Simon's method 
of scaling is remarkably simple, with the transition from second order equations to high order 
systems made by essentially only changing notation. Simon's method, first published in journal 
form in [ Siml] , can also be found in his book [Sim2| (which appeared some years earlier), and 
complete details can also be found in Simon's lecture notes on PDE IISim3ll . In MSimll Simon's 
indicates how his method adapts to encompass equations and systems of parabolic type, and this 
is pursued in Lamm's Diploma Thesis [L |. Before proceeding, we first recall the the Holder space 
interpolation inequality, which we shall use often in the derivation of the Schauder estimates. 

Proposition 3.12 (Holder space interpolation inequality). Let p > 0, e > and Qp C M'^+^. 
Suppose u E C^"^'^'"((5p). There exists a constant C = C{n, d, e, a, m) such that 

<ep'"^+-[d'"^'\]a;Q,+C\u\o;Q, 

holds. 

The interpolation inequality can be established by simple contradiction arguments or directly 
using the mean value theorem; see, for example, |0- 

3.1.3.1. Trudinger's method of mollification. A (parabolic) mollifier (of order 2m) is a fixed 
smooth function p G C^{R'^+^) with / f^^+i pdX = 1. For r > we define the scaled moUifier 

, , 1 fx t \ 

Let P G M"+i and u G Lj^ci^)- For < r < d{X, dP), the moUification of u is given by 

and satisifes spt Ur C Pt, where P^ = {X G P : d{X, dP) > r}. 
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Proposition 3.13. We have Ur G C^. 

Proposition 3.14. Let u G Lj^^.($7). The following estimates hold: 

(3.17) I'"t|0;P, < \u\o;P^ 

(3.18) \did{urix,t)\o-P^ < Ct-'-^'^^uIo-p^. 
Proof. To prove (I3.17t . we have 

1 f f fx — y ^ — 
^ \''\o;p.-:piT2^ J J p[—^^^) dyds 

= \uW,Pr- 

And for KWi : 

dldiur{x,t) = II dial p{^^,^-^^u{y,s) dyds 

<Cr-^-'™>|o;P.. 

□ 

Proposition 3.15. Let u G C^^^{P). The following estimates hold: 

(3.19) \Ur{x,t) -u{x,t)\o-P, <T°[li]„;P^ 

(3.20) \did{ur{x,t)\o-P^ < C7r"-*-2™^M„;P^. 
Proof. For estimate (I3.19l l we have 

Ur{x, t) - u{x, t) = II p (^—^^ ~2^^ (^(y^ ^) ~ ^(^' *)) dy d^ 

< OSCp^ U 

To prove the second estimate we have 

didlur{x,t) = -^ II dldip{^^,^-^^u{y,s)dyds 

= :^m^ I f ^xdip {^^, ^) {<y. s) - u{x, t)) dy ds 

+ dyds. 

The mollifier p is has compact support on Pr and so the last term vanishes by the Divergence 
Theorem. Continuing, we have 

dld{ur{x, t) = 11^ dldlp (^^-y^, iu{y, s) - u{x, t)) dy ds 

□ 

To motivate things a little in the pai^abolic settting, we first briefly show how Trudinger's 
method works in the elliptic setting by treating the Poisson equation. The crucial ingredient in 
Trudinger's method is the following norm equivalence: 
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Lemma 3.16. Let u € C7"(R'^), R> Q and a G (0, 1). There exists constant C = C{d, a) such 
that the norm equivalence 

y;[u]a;BB < sup T^'"'\dxUr\o;BR < C[u]a-BR- 
^ 0<T<R/2 

is valid. 

Proof. The inequality on the right follows directly from equation (13.201 ) (the elliptic version) 
by choosing the appropriate values for the indices i: choosing i = 1 (there is no j in the elliptic 
moUifier) gives 

The first inequality requires a little more work. Let x, y G M*^ and r € (0, R/2). For |x — y| < R, 
by the triangle inequality 

\u{x) - u{y)\ < \u{x) - Ur{x)\ + \ur{x) - Ur{y)\ + |n^(y) - u{y)\ 

< 2T°'[u]a;BR + \dxUr\o-BR\x - y\. 

Set r = e\x — y\, where e < 1/2. Factoring out and dividing by \x — we find 

Choosing e < (l/2)~" and taking the supremum over r S (0, i?/2) completes the proof. □ 

We now derive the Schauder estimate for Poisson's equation. For simplicity we consider 
solutions with compact support in M.^ (the techniques for treating the general case will be seen 
later on when we treat parabolic equations). Fix a G (0, 1) and suppose thatu G Cq'"{R'^) solves 

-a'^ {x)diju{x) = f{x), 

where we assume a*-', / G C"(E'^) and A|^p < a^^iiij < A|^p. We proceed by the method of 
freezing coefficients, and accordingly fix a point xq G M*^ a rewrite the above equation equation as 

-a'^{x^)d,ju{x) = {a'^{xo) - a'\x))d,jU + f{x) 
(3.21) := g{x) 

By a linear coordinate transformation we can assume a^^ {xq) = 5"^^ so that equation (13.211 ) be- 
comes the Poisson equation. We now mollify equation (13.21b to get 

-An^ = gr 

and then differentiate thrice with respect to x to obtain 

-Aa^u^ = dlgr. 

We choose a radius i? > and work in the ball Br. Using inequaUty (13.201 ) we can estimate 

\dlgr\Q-BR < C{n)T-^\g\Q.BR+^ 

<C{n)T-\R + Tr[g],,Bn^^ 

< C{n)T~\R + Tr{[a]^,BR+AdluW,Bn^^ + 

We now recall the solid mean value inequality for subharmonic functions: If v solves —Av{x) < 
on a ball Bji{x) C M'', then v satisfies 

. N C(n) f , , , 
v{x) < — — / v{y)dy. 
^ ' - R- Jb^ 

To apply this inequality to our situation, noting A|xp = 2n, we have 

-A (dlur + = -Mlur - \dlgrW,B, < 0. 
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Thus the function d^Ur + |9^5r|o;BflkP/(2n) is subhamionic and applying the mean value in- 
equality and estimating we obtain 



mur{xo)\ < C{n) ( R-" 
1 



d Ur{y)dy 

Br 



< C{n) {R''-\dlu]^,BR+T-''R\R + Tr{[a]^.,BR+^dluW,BR^^ + . 
Setting R = Nt and returning to the original coordinates we find 

T^-''\dlur{xo)\ < C{n,X,A,a) {N^-^d^^uj^-SR + N^^^{[a]^,BR+^\dlu\o-BR+^ + [f]a;BR+^)) 

Now taking the supremum over r > and using the norm equivalence we obtain 

Choosing N sufficiently large and using the Holder space interpolation inequality on the right 
gives the desired estimate, namely 

where C depends on n, A, A, and a. Having given a feel for Trudinger's method, we move on to 
use this method to derive the Schauder estimates for second order parabolic equations. The crucial 
equivalence of norms lemma in the parabolic setting is the following: 

Lemma 3.17. Let u G C"(M'^"'^^), R > and a G (0, 1). There exists constant C depending only 
on d and a such that the norm equivalence 

y;[u]a;QR < sup { T^"" | S^U^ |o; + T^"""" | StU^ |o;Q^ } < C[u]a,:Q^. 
0<T<B./2 

is valid. 

Proof. The second inequality follows directly from equation (13.201 ) by choosing the appro- 
priate values for the indices i and j. To prove the spatial part of the second inequality, choosing 
i = 1 and j = in estimate (13.201 ) gives 

The temporal estimate follows similarly. Let X,Y £ W'-^'^ and r G (0, R/2). For d{X, Y) < R, 
by the triangle inequahty 

\u{X) - u{Y)\ < \u{X) - UriX)\ + \ur{Y) - u{Y)\ + \ur{x, t) - Ur{y, t)\ + t) - Ur{y, 

< 2T"[u]a-Qji + y\\dxUr\0;Qji + I* " s| | ^tU^ 1 0; • 

Set T = ed{X, Y), where e < 1/2. Factoring out d{X, y)" we have 

HX) - U{Y)\ < d{X, Yr {2e^[u]a;QR + e^-W'-^\dxUr\o;QR + e'-^'^r^— | ^tU, | O; Q«) . 

The proposition follows by fixing e sufficiently small and taking the supremum over r G (0, i?/2). 

□ 



We now proceed similarly to Poisson's equation to derive the Schauder estimate for the non- 
homongeneous heat equation. Fix a G (0, 1) and suppose thatn G C7^'°(M'^+i) solves 

dtu{x,t) - a'^{x,t)diju{x,t) = f{x,t), 

where we assume a^^ , f G C"(M'^+^) and A|,^p < a^^iiij < A|,^p. Again we freeze coefficients 
at a point (xq, to) £ K'^^^, perform a coordinate transformation and mollify the equation to get 

(3.22) dtu-r — Au^(x,t) = Qr- 
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Given the form of the norm equivalence, the desired Schauder estimate will follow if we can 
establish the estimates (for the spatial component of the Schauder estimate) 

(3.23) \d^Ur{xo,to)\ < C (^^OSCQj^dlu + R^\dlgr\o;Qi 

(3.24) \dtdlurixo,to)\ < C (^-^ oscq^dlu + R^\dtdlgr\o;Qii^ 
and for the temporal part 



1 



(3.25) \d:rdtUr{xo,to)\ <C l^— oscQ^ dtu + R \dxdtgr\o;Q 

(3.26) \dtUr{xo,to)\ < C (^-^OSCQ^dtU + R'^\dfgr\o;Qi 

We show how to obtain the spatial estimates, as the time estimates follow in exactly the same way. 
We recall the mean value property for subsolutions of the heat equation: If is a subsolution to 
the heat equation on M'^"*"^, that is if v satisfies dy — Av < 0, then v satisfies 

1 /" /" / xko-yp 



"^r"- J JE{xo,to;r) [to-sy 

for each £'(xo, to; C M'^+^. Recall the heat ball E{x, t; r) is the set given by E{xq, to;r) = 
{{y, s) e W^+'^ : \xo - < A/-27rs log[r2/(-47rs)], s E (to - rV(47rs), to)}. We denote the 
radius of the heat ball by Rr{s) : y^—2TTs log[r2/ (— 47rs)]. For further information on the mean 
value property of the heat equation we refer the reader to [Ev l and [E|. Let us now show (13.231) : 
Differentiate (13.221 ) thrice in space. Since \d^gr\o; eIx]"^ / {2n) is independent of time we see 

dt (^dlur + \^x9t\0;E^-^^ " A (^d^Ur + l^xS'r |o; = dt{dlur) - A{d^Ur) - \d^gr\o; E 

= dlgr - \dlgr\o;E < 0, 

and hence the function d^Ur + {d^grlo; E{xo,to;r)\^\'^ / C^^) subsolution of the heat equation. 
From the mean value property of subsolutions we have 

(3.27) dlur{xo,to) < 7^ // {dlur{y,s) + \dlgr\o\yf) ~ dj/dg. 

4^ J J E{x,t;r) — S\ 

By translating coordinates we can assume that (xo, to) = (0, 0). All the desired estimates involve 
evaluation the integral 

1 [ Rr{s) 



An 



■ ds, 



where a and /3 are given integers. The constants can be computed explicitly, however we are only 
interested in the scahng behaviour with respect to the radius r (and that the integral is finite). We 
compute 



1 Rr{sT 1 (-2nslog[rV(-47rs)])"/^ 



47r 47r 



C(n, a, j3)r~ 



n+a-2p+2 ^ " ^a/2-/3 ( log(4vrt)) dt 



Att 

oo 



With further substitution this integral can be converted into the Gamma function, which is finite 
as long as a/2 > —1. Returning to (13. 271 ). we have 

(3.28) dlur{xoM<^r-^jjjlMy,s)^^dyds + A\dlgrW,Er-''jjJ^ 
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We estimate the first term on the right by 

|2 



4r- [ [ d^yuAy, sy-\ dyds < Cr-^ [ j S^n, dy ds 







47r 








47r 


S2 



< Cr"" / ^ ( / oscd^udy ds 



<Cr ^oscEd'tu I ' ds 



S2 



C{n) ^2 
r 



The second term on the right of (I3.28l l can be estimated more simply to give 

4\dlgrWEr-'^ jjj^dyds < C{ny\dlgrW,E. 

The estimates involving time derivatives can also be estimated in a similar manner. For example, 
by integrating by parts in time, we have 

,y|2 



4r-" j I dtd^Uriy, s)^ dyds < Cr-'' J j osc 



u- „ 

C{n) ^2 
< — — osce o^u. 

The derivation now continues in the exactly the same was as for the Poisson equation, using the 
estimates (13.231 ) - (13.251) . the equivalence of norms lemma and the Holder space interpolation 
inequality; we ultimately obtain the desired Schauder estimate: 

(3.29) [5'''«]a;R^ < C{[f]^.^, + |u|o;M.), 

where C depends on n, A, A, and a. The method extends to more general equations and domains 
by using cutoff functions and Simon's absorption lemma, as we shall soon see in the case of 
systems. 

3.1.3.2. Simon's method of scaling. As we have mentioned before, Trudinger's method ex- 
tends to systems of even order, where the application of the mean value inequalities are replaced 
by estimates and the Sobolev embedding theorem. For parabolic systems the method becomes 
a httle computationally cumbersome, and instead we shall use the Simon's method of scaling. 
For the derivation of the Schauder estimates for elliptic systems, in addition the Simon's original 
paper HSimll . we highly recommend his lecture notes on PDF IISim3i . Once one has defined the 
notion of a parabolic polynomial his method adapts immediately to parabolic systems. Here we 
simply quote the interior and near-bottom Schauder estimates for parabolic systems of even order 
on Fuclidean space, and refer the reader to |0 for complete proofs. Any errors or inconsistencies 
are due to us. 

Proposition 3.18 (interior Schauder estimate). Suppose u € C'^'^'^''^{Qr{Xq)) is a solution of a 
general linear 2m-order parabolic system 

(3.30) Lu" := dtu'' + {-1)"" ^ Af {x,t)diu^ = f . 

\I\<2m 



Suppose the following conditions are satisfied: 

\CLiljl---imjm 

1^ ^nii^jic^ iitc^ oyitiitvK.li y i..utiLiiiiuti .n.^ 



1) The leading coefficient A?^^''^ *mjm ^^^/^yj^^ ff^g symmetry condition ■ 



bjiii---j,„in 



2) The leading coefficient satisfies the Legendre-Hadamard condition with constant A 

3) There exists a uniform constant A < oo such that X]|/|<2ml^^U;Qi{(Xo) — ^- 



3.1. SHORT-TIME EXISTENCE FOR FULLY NONLINEAR PARABOLIC SYSTEMS 



31 



Then there exists a constant C = C{n, N, 9, A, A) such that the estimate 

[d'^'Ma; eQa < C{[fUQn + i^-^'""" |^z| 0; Q,) 

holds for each 6 £ (0, 1). 

Proposition 3.19 (near-bottom Schauder estimate). Suppose u e C^™'^'"((5^(^o)). withu{-,0) 
uq, is a solution of a general linear 2m-order parabolic system 

(3.31) Lu" := dtu"" + {-ir Yl Af{x,t)diu^ = r. 

|7|<2m 



Suppose the following conditions are satisfied: 

t a 



1) The leading coefficient ' satisfies the symmetry condition y^^ui'"*™-?™ 



bjiii—jmiv 



2) The leading coefficient satisfies the Legendre-Hadamard condition with constant A 

3) There exists a uniform constant A < oo such that X]|/|<2ml^^la;Qii(^o) — ^• 
Then there exists a constant C = C{n, N, 6, A, A) such that the estimate 



u\ 



holds for each 9 £ (0, 1). 

The above estimates are the localised counterparts to equation (I3.29I ). In order to localise the 
estimate, the following adsorption lemma is needed: 

Lemma 3.20 (Simon's adsorption lemma). Let S be a real-valued monotone sub-additive function- 
on the class of all convex subsets of Bji{xo){i.e. S{A) < J2iLi ^i^j) whenever A,Ai,..., An 
are convex subsets with A C ^jLi C Br^xq). Suppose that 9q G (0, 1), /i G (0, 1], 7 > 1 and 
I > are given constants. There exists an eo = eQ{l , 9 , n) > such that if 

p'S{Bep{y))<eop'S{Bp{y)) + j 

whenever Bp{y) C Bf({xQ) and p < pR, then 

R'SiBeRixo)) < C7, 

where C = C{n,6, p,l). 

The proof can be found in HSimll and IISim31 . In localising the Schauder estimate we need 
to apply the adsorption lemma in the case S{A) = [u\a-A- We confirm that the lemma holds in 
this case, that is S is monotone and sub-additive on convex subsets of Qr. Let i? > a given 
radius and A C Qr. Since the Holder constant is defined by taking the supremum over a set, 
monotonicity clearly holds. To show sub-additivity, suppose A C Ai U A2, where all sets are 
convex. Fix X,Y £ A. If either X,Y £ Aior X,Y £ A2, then 

\u(X) — u(Y)\ 

< max|[nJa;Ai, Nq:;A2} < Na;Ai + ma;A2- 



If on the other hand X £ Ai and Y £ A2, the choose Z £ Ai n A2 lying on the line segment 
between X and Y. Then 

\u{X) - u{Y)\ ^ \u{X) - u{Z)\ + \u{Z) - u{Y)\ 



d{X,YY ~ d{X,Z)'^ +d{Z,YY 

^ \u{X)-u{Z)\ \u{Z)-u{Y)\ 
- d{X,ZY d{Z,Y)^ 

< [u]a;A+ [u]a;A2- 

The general case follows by induction. 
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3.1.3.3. Global Schauder estimate. The above Schauder estimate holds on a small parabolic 
cylinder Qi? C M^. We now want to lift these local estimates to the vector bundle E x {0,u) to 
obtain Schauder estimate globally on E"^. Let ^p : V x I ^ M" x M_(_ be the coordinate map for 
a sufficiently small neighbourhood V x I c M x {0,uj). By definition of a vector bundle, there 
exists a bundle trivialisation 4' : E^\vxi V x I x M^. In fact, if i/j is the coordinate map for V, 
then = if) X id. Using the bundle trivialisation and the coordinate maps we can locally identify 
a section T{E x (0,a;)) as a subset of M" x M+ x M^. We will abuse notation slighly, and for 
U G T{E X (0, uj)), we shall write {il)~^)*U to mean the local trivialisation C/|vx/ pulled back to 

X M+ X via the coordinate map ijj. 

Next we want to control the norm of section measured with the bundle metric in terms of the 
Euclidean norm of the pulled-back section. 

Proposition 3.21. Let U G T{E x (0, w)) and {Vi, t/ji) be a covering ofM x (0, cj) by a finite num- 
ber of normal charts of sufficiently small radius Rq. Then there exists a constant C = C{n,RQ) 
such that in each neighbourhood Vi the equivalence of norms 

^l^°^~^l2m,l,a;V^,(y.x7,) ^ \U\2m,l,a;V,xh < ^1 ° V'"^ 1 2m,l,a; (V, x/,) 

is valid. 

Proof. For the parts of the Holder norm involving suprema this is easy to show, as one 
simple writes the covariant derivative in terms of ordinary derivatives and the Christoffel symbols 
and uses Lemma lXT] To deal with Holder semi-norm, we note that parallel translation is defined in 
terms of solving an ordinary differential equation. We then have control on the size of the Holder 
coefficient in terns of the initial condition for the ODE in a finite number of charts, thus it too is 
uniformly bounded. □ 

Using the above lemma, we can now patch together the local Euclidean Schauder estimates to 
give the desired global Schauder estimate. 

Proposition 3.22 (Global Schauder estimate). Let E x (0, uj) be a vector bundle over M x (0, uj), 
where M is a closed manifold. Let L : T{Ei^) — )■ T{Fi_j) be linear differential operator of order 
2m. In any local coordinate chart L is of the form 

(3.32) ^C/- + (_!)- ^ A'diU, 

I<2m 

or in full 

with U{-,0) = Uq. Suppose that in any coordinate chart the following conditions are satisfied: 

1) The leading coefficient satisfies the symmetry condition ■ = 

2) The leading coefficient satisfies the Legendre-Hadamard condition with constant A 

3) There exists a uniform constant A < oo such that Yli\i\<2m 

Then there exists a constant C = C{n, N, A, A, M, uj) such that the estimate 

2m, 1, a; Eui 

Proof. Because M is compact, we can cover M x (0, uj) by a finite number of coordinate 
patches (Vijipi) of sufficiently small radii R-i < Rq so that we can apply Proposition (321] Suppose 
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X, y G are any two points. If d{X, Y) < Rq, then we estimate 



<C^(|F|o,o,a;y, + |C/|c 



|0;K 



<C(|F|o,0,a;i^. + |C/|o;i?J. 

On the other hand, if d{X, Y) > i?o we estimate 

<0|V ' U\o;E^Ro 



(3.33) 



< C'|C/|2m,l,a;Si^ 
<C{\F\o,0,a;E^ + \U\o;E^). 

Note that we have again used the fact that we have a finite covering, as we have needed to take the 
supremum over the all Holder coefficients in each chart. □ 

3.1.4. Linear existence theory. The next step in our existence program is to prove existence 
and uniqueness for linear operators in Holder space. We begin with the 2mth order heat operator. 

Proposition 3.23. Consider the following initial value problem: 

idtU + (-A™)C/ = X G 

\C/(-,0) = [/o. 

Suppose that F € C°'°'"(£;<^) and Uq G C2™'1'"(£;^), where a G (0, 1). Then problem (1333]) 
has a unique solution U G C'^'^'^'°'{E^). 

Proof. As usual, we may assume without loss of generality that Uq = 0. By mollification 
we can construct a section G C°°{Ei^) such that 

\Fe\a;E:^ < 'i\F\a;E^- 

Now consider the approximate problem 

(dtU, + {-A^)U, = F,{X), XGM^ 
\C/(-,0) = 0. 

From the theory, there exists a unique smooth solution G C°°{E^) to the above approximate 
problem. A short contradiction argument (see [|Sim3J) shows we can estimate 

\Ue\o:E^ < (■\Ue\2m,l,a;E^ + c{e)\\Ue\\L^{E^), 

and then using the Hilbert space regularity estimates we may estimate 

\\Ue\\L^(E^) < C\\F^\\j^2(^E^) < C\F^\o,0,a;E^- 

We point out that in the case of second order equations, using the maximum principle it is a slightly 
simpler matter to estimate 

mo;E^<C\F,\o;E^<C\F,\^,E^. 

Combining this estimate with the global Schauder estimate, \Ue\2m,i,a;Eu, < C{\Ue\o;Eu,+\^e\a;Eu,)^ 
we get 

\Ue\2m,l,a; E^ < C\F\a-E^, 

where the constant C is independent of e. Given that F G C"{Ei^), the left hand side is uniformly 
bounded. The Arzela-Ascoli theorem now applies to give a subsequence such that C/^ — )• f7 
uniformly in C'^''^'^{E^) as e — )• 0, and moreover U G C^™'^'"(£'tj). Last of all, we show 



(3.34) 
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uniqueness by the energy method. Suppose that Ui and U2 are two solutions to (I3.33l l. and consider 
the problem for W := Ui — U2, where W now solves the homogeneous heat equation with zero 
initial condition. For < t < T we define the energy e(t) by 



Then 



e{t) = / W{x,tydVg. 

^At)=2 f WWtdVg 

at Jm 

= 2 [ W{-A"')WdVg 
< 0. 

Thus e{t) < e(0) for all t G [0, T], and consequently Ui = U2 and the solution is unique. □ 

With a solution to the heat operator in place we can now use the method of continuity to solve 
the general linear problem. 

Theorem 3.24 (method of continuity). Let B be a Banach space, V a normed linear space, and 
Lq and Li bounded linear operators from B to V. For t G [0, 1] define 

Lr := (1-t)Lo + tLi 

and suppose there exists a constant C such that the estimate 

\\u\\b < C\\Ltu\\v 

holds independent ofr. Then Li maps B onto V if and only if Lq maps B onto V. 

For a proof of the method of continuity we refer the reader to [GT, pg. 75]. 
Proposition 3.25. Consider the following initial value problem: 

(dtU- + {-irE\i\<2^m^fiX)ViU' = F-{X), X€M^ 

\ui;0) = Uo. 

Suppose that the following conditions are satisfied: 

1) The coefficients Af^^^'"^""^"^ satisfy the symmetry condition = 

2) The leading coefficient satisfies the Legendre-Hadamard condition with constant A 

3) There exists a uniform constant A < 00 such that '^\i\<2ni 

Then problem (13.331 ) has a unique solution U G C^™'^'"(£'a;). 

Proof. As always, we may assume without loss of generality that Uq = 0. Define the 
operators 

Lq = dtU" + {-A"')U'' 

Li = atC/'^ + (-1)™ Yl K\x)ViU\ 

\I\<2m 

Consider the family of equations 

L^ := (1 - t)LqU + tLiU = F, 

where r is a parameter with r G [0, 1]. The operator Lj- satisfies the assumption of the theorem 
with At- and A,- taken as A,- = min{l, A} and A,- = max{l, A}. Suppose that Ur is a solution to 
(13.351 ). Then in exactly the same way as for the heat equation, using regularity, the same short 
contradiction argument and the global Schauder estimate we obtain the estimate 

\Ue\2m,l,a; Eu: < C\F\a;Eu:^ 

where C is independent of r. We may now apply the method of continuity, and since Lq is solvable 
by Theorem 13. 231 Li is also solvable. □ 



(3.35) 
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3.1.5. Nonlinear existence theory. With all the linear existence theory now in place, we are 
ready to prove Main Theorem [T] We do this by an appliation of the inverse function theorem in 
Banach spaces. 

Theorem 3.26 (inverse function theorem). Let X and Y be Banach spaces, P : X ^ Y a map 
from X to Y, and Uq and element of X. Suppose that P satisfies the following: 

1) P is continuously differentiable at Uq 

2) The Frechet derivative of P at Uq is invertible. 

Then there exists an open neighbourhood 11 ofUo in X, and an open neighbourhood VofV := 
P[C/o] in Y such that P : U ^ V is an isomorphism. 



For a detailed proof of the inverse function theorem we recommend to the reader HAEI pg. 
215]. 

Proof of Main Theorem^ Let X = C'^'^^^^'^ {E^) and Y = C^iE^), where /3 < a. We 
consider the nonlinear operator P as a map P : X ^ Y.To begin, linearise the nonlinear operator 
P at the initial value Uq- The linearisation of P about Uq in the direction V is a linear system 
in the unknown V which uniquely solvable by the Schauder theory presented in the previous 
section. Call Ui the solution to this linear system. From the Schauder theory we also know 
Ui G C^'"'^'"(£'(^). Now linearise P about the solution to the linear problem Ui. Next we confirm 
that the conditions of the inverse function theorem hold for the nonhnear operator P about Ui. The 
(Gateaux) derivative of P at in the direction V is given by 

P'iUi)V = ^F{Ui + sV) 

OS s=0 

= dtV - F^^'-'^^-{x, t, Ui,VUi, V^Z.,2rr.UiWn ' ' ' ^^2^V 
(3.3d) . 

• • • - F'^ix, t, Ui,VkUi, V2™C/^)V,,y - F{x, t, Ui,VUi, V2™C/;)F. 

The regularity assumptions in the statement of the theorem ensure that P is continuously differen- 
tiable and Frechet differentiable. We have 



\P{Ui + V) - P{V) - P'{Ul)\o,0,a;E^ 

^ P'{Ui + sV)- P'{Ui)ds)v 



< \\P'{Ul + sV) - P'(f/;)||£(c2m,l,a(£;^)^c'0'0."(£;„))l^|2m,l,Q;£;„ 
= o{\V\2m,l,a;Ej- 

Because F is continuous in all its arguments, 

\\P'iUi + sV)-P'{Ui)\\^c^^,i.c(^E^)^co-°'-'{E^)^0 as s^O 

and the last line above follows. This shows that P is Frechet differentiable at [//. The hnearisation 
of P about Ui in the direction V is again a linear system in the unknown V that is uniquely solvable 
by the Schauder theory, and thus the Frechet derivative of P is invertible at Ui. 

The inverse function theorem applies and guarantees an open neighbourhood IX of Ui in X, 
and an open neighbourhood "V of P[Ui] in Y, such that P : U ^ V is an isomorphism. For 
convenience, set fi{t) := P[Ui]. Define the function /^(t) := x{'t)fi{'t), where x{t) is a smooth 
cutoff function with the properties xi't) = for t < and xi't) = 1 for t > te, and is small 
number to be fixed sufficiently small. We claim for t G [0, t^) where is sufficiently small, that 
is in V.Beginning with the supremum estimate, if t > t^, then \ fi — fx\o = 0- For t < we 
use the crucial fact that since Ui is the solution to the linear problem, // satisfies fi{0) = 0: 

\Mt) - fxit)\ = \fiit) - fx(.t)\ - l//(0) - /x(0)l 
< [fi-fxUP, 
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and so \ fi — fx\o ^ Ct"^^'^"^^ since Ui is Holder continuous. For the Holder estimate we consider 
two cases. We may assume without loss of generality that t > s. If |t — s| < t^, then we need 
to consider two further subcases: I) s < with < s < t < 2t^; and 2) t,s > t^. In the first 
subcase we begin estimating 

\Mt) - fxit) - (Ms) - /x(^))l < \Mt) - fi{s)\ + l/xW - /x(^)l 

< {[flU + \x\o[fl]a + [xUfl\o)\t-s\^. 

The second term on the right is easy to deal with, since |x|o < 1- To deal with the second, we note 
that in this case we can estimate 

Combining estimates we see 

\fi{t) - fxit) - (Ms) - fx{s))\ < C\t - s\^tP 

where (3 < a, and so [// — < Ct^" second subcase is easy, since if t, s > t^, 

then \fxit) — fx{s)\ = \fi{t) — fi{s)\- To treat the second main case, namely if |t — s| >te, then 

\Mt) - fxit) - (Ms) - fx{s))\ < \Mt) - fx{t)\ + \Ms) - fx{s)\ 

< m - fx\o 

< 2Ctf^ 

< C\t - Sl^^te^"" . 

Therefore \ fi — f^lfs can be made arbitrarily small on small time intervals, and so we can fix 
sufficiently small so that for all t G [0, t^), is in V. By the inverse function theorem there exists 
a unique element ^ X such that -P[C/x] ~ /x' ^^'^ moreover, for t < P[Ux\ = 0. Thus 
the element is the unique solution to the initial value problem (13.51 ) for some short-time te/2 
and the proof is complete. □ 



3.2. Short-time existence for the mean curvature flow 

Here we apply the nonlinear existence theory espoused in the previous section to give a proof 
of short time existence of the mean curvature flow. In this section we denote the mean curvature 
flow, considered as a differential operator, by M, and the mean-curvature-DeTurck flow by M D. 
We begin by showing that mean curvature flow is only a weakly parabolic quasilinear system, and 
as such we cannot immediately apply the 'standard' theory. With respect to the induced metric the 
Laplacian of F is just 

^ \dx^dxi '^dx^J 
= H. 

The mean curvature flow equation can therefore be written as 



The similarity is however deceptive: The induced metric is evolving in time, and this adds extra 
terms to the principal symbol that result in the presence of zeroes. The principal symbol can be 
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computed by 

^ ^ ydx^dx^ dx^ ) 



ij f a^F" _ 1 f ^ _d_ \ dF^\ 

^ ydx^dxj 2^ ydx"^^^^ dxi^^^ dx^^^'^ J dx'^ ) 

(7 — ^ — r g -^y — ■■ — 1 r : h • • • 

dx^dx^ 2 dx^ \dx^dx^ dx^ dx^ dx^dx^ 

a^F-^ .idF'^dF'' d'^F^ 
9^^rT^-9^9 



dx^dx^ dx^ dx^ dx^dxi 

^kl dF'^ dF'' 



dx^ dx^ J dx^dx^ 



Observe that the term 9^;k^r is the orthogonal projection onto the tangent space of the sub- 
manifold: for any ^ G 



kl ca 



dx^ dx^ dy^ 

To examine the principal symbol, without loss of generality we may assume at a point that gij = 
dij and also that |^| = 1, so we can choose d = 1 and = for i > 2. The principal symbol is 
thus 

mm = \c\''{id-7rT^{0) 

= lel'vr,vs(0, 

which is zero if ^ E TE. Another way to see that the mean curvature flow is only weakly parabolic 
is to observe from the start that the equation is degenerate in tangential directions. We have just 
computed that 



^dx^dx^ dx^ 
so the mean curvature flow can also be written as 



-F = TTATS 5^- 



For any ^ G 



pn+fc 



dt \ dx^dx^ 



d ^.i.^dF^dF^ d 



Qya Q^k Qyb 

so again we find the mean curvature flow is given by 

dF" iif.a kidF''dF''\ 92 F^ 



dt \ dx'^ dx^ J dx^dx^ 

The mean curvature flow is therefore not strongly parabolic and the almost standard parabolic 
theory cannot immediately be conjured to yield existence for a short time. To overcome this diffi- 
culty we are going to adapt a variant of the DeTurck trick first elaborated by Hamilton IIHa6l that 
combines the mean curvature-DeTurck flow and the harmonic map heat flow. As the next propo- 
sition shows, the mean curvature flow is invariant under a tangential parametrisation. This means 
that adding a tangential term to the mean curvature flow equation results in a solution that differs 
from the solution of the mean curvature flow itself only by a reparametrisation of the submanifold. 
The DeTurck trick involves adding a tangential term to the mean mean curvature flow to break the 
geometric invariance of the equation. The modified flow is then strongly parabolic and the almost 
standard parabolic theory can now be summoned to ensure short time existence. The solution to 
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the mean curvature flow is then recovered from the solution to the mean curvature-DeTurck flow. 
Hamilton's coupling of the modified flow with the harmonic map flow serves to provide a simple 
proof of uniqueness. 

Proposition 3.27. Let W be a time-dependent family of vector fields defined onTi x [0, T). Sup- 
pose that F is a solution to 

lw = ^aF + VwF 
\f(-,0) =Fo. 

Then there exists a solution F to the mean curvature flow with Fq = Fq. 

Proof. For the moment, assume that there exists a time-dependent family of diffeomor- 
phisms ipt : T, X [0, T) — )■ S. Computing in local coordinates {x*^} around ipt{p) we calculate 

dF, , OF, , , , ^ , , , dipAp)^ 
-g^{p,t) = -^{^t{p),t) + VuF{^t{p),t) ■ 

= AgF{^t{p),t) + {w\^t{p),t) + ^^^) VkF{^t{p),t) 

= KgF{p,t) + [w\^t{p),t) + ^^^) VkF{^t{p),t). 

Therefore, if we can show there exists a family of diffeomorphisms solving the initial value prob- 
lem 

['^ = -W{vt{plt) 
\vq{p) = ids, 

then F will be the desired solution to the mean curvature flow. In the case that S is compact, 
standard ODE theory (for example, see ULel ) guarantees that the above ODE problem has a unique 
solution for as long as W is defined. □ 

Let us now continue with Hamilton's argument. Fix a background connection V on S. For 
example, we could take the induced connection on S at t = 0. As the vector field W in the above 
proposition we take W := g^-^{T^j — F^^). Consider the mean curvature-DeTurck flow given by 



dt 



AgF" + VivF"" 



ydx^dx^ '^^ dx^ 
The principal symbol is now 

a[MDm = lei'id, 

so the mean curvature-DeTurck flow is strongly parabolic and Main Theorem [T] guarantees a 
unique solution to this modified flow for a least some short time. The conditions of Main Theorem 
[T]are easily confirmed for the mean curvature-DeTurck flow. For example, the leading term of the 
linearised operator in some direction V is given by 

dV d'^V 

— -g^ 



dt dx^dx^ 
d'^V^ 
dx^dx^ 

Hence A'^^-' = g^^S^, and A^*-' = a!'J^ ■ As the mean curvature-DeTurck flow possesses a unique 
solution for some short time, the family of vector fields W{t) also exist on this short time interval, 
and the above ODE problem has a unique solution on the same time interval. By Proposition 
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13.271 we recover a solution to the mean curvature flow F by puUing-back the solution of the mean 
curvature-DeTurck flow by the diffeomorphism (ft, that is F{p, t) = ipfF{p, t) = F{(pt{p),t). 

We now show uniqueness of the above solution to the mean curvature flow. Suppose that F is 
a solution the mean curvature flow and denote associated the induced metric by g. Let 999 : S — )■ S 
be a diffeomorphism. Fix a metric g and associated Levi-Civita connection on the target manifold 
S, and consider the harmonic map heat flow 

with respect to the domain metric g and the target metric g. The harmonic map heat flow is a 
strongly parabolic quasilinear system (see, for example, HHall or IIHa4ll ) and thus possesses a 
unique solution ipt for at least some short time. We now define F := ipuF = {ip~[^)*F and claim 
this is a solution to the mean curvature-DeTurck flow. Repeating the calculation in Proposition 
l3.27l shows 

where V{p) = — |<^-i(p)- Thus if we can show that V = W then this estabhshes the claim. 
This follows from following result: 

Proposition 3.28. Suppose {K^, k), {M^,g) and (A^™, h) are manifolds, il) : K ^ M a diffeo- 
morphism and (p : M ^ N a map. Then 

The geometric meaning of this proposition is that the harmonic map Laplacian from a domain 
manifold to a taiget manifold is unchanged if we reparametrise the domain manifold. For a proof 
of this proposition we refer the reader to IICLNI pg. 117] or IIHAI pg. 78]. We have V{p) = 
— (^*9f|^-i(p) = —Ag^g(p{ip^^{p)) = Ag gids, and then adapting the above proposition to our 
setting {il; = ip~^ , h = g) we see A^^^ids = g^H^ij ~ ^ij)' ^^^^ "^^e two vector fields V and 
W are in fact identical. 

We can now finish the uniqueness argument. Suppose that there exist two solutions Fi,i = 1,2 
to the mean curvature flow with initial condition Fi{-,0) = F2{-,0). For each domain metric gi 
we can solve uniquely the harmonic map heat flow problem 

= ^h,9^ 

¥^(-,0) =ids 

for the functions (pi, which then give two solutions Fi = 99* F to the mean curvature-DeTurck 
flow. Because these two solutions satisfy the same initial condition and solutions to the mean 
curvature-DeTurck flow are unique, Fi = F2. The two diffeomorphisms Piit) also solve the same 
ODE problem 

i^ = -W{cp,{p,t),t) 
\v>iiP,0) =P- 

and so they too are in fact equal on their common interval of existence. Therefore Fi = (p^Fi = 
'P2F2 = F2, which concludes the proof of uniqueness. 



CHAPTER 4 



Submanifolds of Euclidean space 



Our goal in this chapter is to prove Main Theorem 2: 

Main Theorem 6. Suppose Sq = Fq{Y7^) is a closed submanifold smoothly immersed in 
IfTiQ satisfies \H\min > and < c|-ffp, where 



c < 



if2<n<A 
ifn > 4, 



then MCF has a unique smooth solution F : S x [0, T) — )• W^^^ on a finite maximal time interval, 
and the submanifolds St converge uniformly to a point q £ M"+'^ as t ^ T. A suitably normalised 
flow exists for all time, and the normalised submanifolds converge smoothly as t ^ oo to a 
n-sphere in some {n + k)-subspace of^""^^ 



4.1. The evolution equations in high codimension 

We begin by deriving evolution equations for various geometric quantities; of particular im- 
portance are the evolution equations for and The mean curvature flow amounts to the 
prescription F^:dt = iH in the notation of the previous chapter. For the moment we allow the back- 
ground space to be an arbitrary Riemannian manifold. The timelike Codazzi identity ( 12.191 ) is 
precisely the evolution equation of the second fundamental form under the mean curvature flow: 

(4.1) Va,/i(n, v) = Vu^vH + h{v, W(u, H)) + vr [R{F^u, lH)F^v) , 
or with respect to arbitrary local frames for the tangent and normal bundles 

dt^ij = ViVjH + H ■ hiphpj + H"'Riaj^h'i3. 

Using Simons' identity (I2.24I ). this converts to a reaction-diffusion equation 

^ dt^ij — ^^ij ~1~ l^ij ' ^pqhpq ~t~ h^q * hqphpj ~t~ hjq • hqphpj ^hip • hjqhpq 
-\- IRipjqhpq Rkjkphpi Rkikphpj -\~ hijaRkcikjB^ j3 
— "^hjpaRipapl^p — 'ihipaRjpapT^li + ^k^kijli^P — '^iRjkkl3^l3- 

For the remainder of this chapter we are concerned only with the case N = M"+'^, in which case 
the equation becomes 

(4.2) ^ dthij — ^hij -\- hij ' hpqhpq -\- hiq * hqphpj -\- hjq ' hqphpj^ 2hip * hjqhpq. 

Taking the trace with respect to g we obtain an evolution equation for the mean curvature vector: 

(4.3) V9,H = AH + H -hpqhpq. 

To derive the evolution equation for first recall that Vtg = 0, and then at a point we compute 

dt\h\^ = dt{h,h) 

= 2{Vthij,hij) 

— 2{Ahij -\- h{j * hpqhpq -\- h{q * hqphpj -\- hjq * hqphpj 2h{p * hjqhpq^ hij^ . 

We now use A|/ip = 2(A/ijj , /ijj) + 2|V/ip, and then noting that three of the reaction terms factor 
into the normal curvature we obtain 

= A|/ip - 2|V/ip + 2^ +2 (^'^ hipahjpi3 - hjpahipi3^ . 

a, 13 i,j i,j,a,^ P 
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Similarly, using equation (14.31 ) the evolution for is given by 

d_ 
dt 



t,j a 

The last term in (14.11 ) is the length squared of the normal curvature, which we denote by For 
convenience we label the reaction terms of the above evolution equations as follows: 

Rl = ^ hijahijpj + 

i,j a 

The special connections we have been using are especially convenient for deriving the evo- 
lution equations in high codimension. This will become quite evident when we come to deriving 
the higher derivative estimates. Of course the special connections do not have to be used, and the 
methods used in the hypersurface theory can still be applied. Let us see how some of this works 
in high codimension. Since the ambient metric is fixed, the evolution of the induced metric can be 
computed by 

lg,^=dt{d,F,djF) 

= {di{H''ua),djF) + {i^j), 

then using the Weingarten relation: di^a = C^^vp — hipaQ^'^dqF and noting which terms are 
orthogonal to each other we have 

^ Q,j = -{H^{hipo.9'"'d,F),djF) + (i o j) 



dt 

= -2H ■ h 



V 

To easily derive further evolution equations in this way it becomes necessary to compute in a suit- 
ably chosen evolving local frame for the normal bundle. Since the normal bundle of a hypersurface 
in one-dimensional, any rotation of the normal bundle is necessarily tangential. In arbitrary codi- 
mension however, the normal vectors may 'twist' inside the normal bundle giving possibly both 
tangential and normal motion. Here we have 

^i^a = {dti^a,dpF)gP'idgF + {dtVc,v^)v^ 

= -{va, dtdpF)gP'idqF + {dtUa, 

= -{vo^, V^H)gP''dgF + {dtu^, v^)v^. 

Observe that mean curvature flow of the submanifold only imposes the tangential motion of 
the normal frame and so we are free to choose the normal motion. A convenient choice is of course 
that there is no normal motion. 

Lemma 4.1. Let i^a(O), n^\<a<n^-p,hea local orthonormal frame for the normal bundle 
and define the evolution of the frame by 

^^u^{t) = -{v^y^H)g^^dqF 
Then Va{t) remains a local orthonormal frame for the normal bundle as long MCF has a solution. 



Proof. We first note that the evolution of the frame is determined by an linear system of 
ODE's and hence has a unique solution as long as MCF has a solution. To show that the frame 
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remains nomial we compute 

j^{dkF, Mt)) = (dtdkF, i^a{t)) + {dkF, OtMt)) 

= {dkH, u^{t)) + {dkF, -{v^, V^H)gP'idgF) 

= {ViH-F4W{dk,H)),iy^{t)) - {ua,y^H)g^<igk, 

= -{F,{W{dk,H)),iyUt)) 

= -gP%H,h{dk,dp)){d,F,u^{t)). 

We set Yk = {dkF, Va{t)) and = —gP'^{H, h{dk, dp)), and the last line above now reads 
(4.4) |n = AIY,. 

Equation (14.41) is a homogenous linear system of ODE's with initial conditions Yfe(O) = and 
y^'(O) = 0, and thus its unique solution is given by Yk{t) = for all time as long as MCE has a 
solution. Erom this we conclude that if the frame is initially normal then it remains so. To show 
that the frame remains orthonormal we easily compute 



= 0. 



□ 



In the coming sections the reader will note that by using the special connections we avoid 
needing the evolution equation for the Christoffel symbols. In high codimension should one wish 
to commute the usual partial derivative in time with spatial covariant derivatives, it is also nec- 
essary to understand how the normal connection forms evolve. By differentiating the Weingarten 
relation in time and using the special evolving normal frame one finds the normal connection 
forms evolve by 

d 

Ol^ka = -(^a • ^^pH) hkp ■ + {1^13 ' ^pH) hkp ■ l^a- 

Note that the evolution equations for the Christoffel symbols and the normal connection forms are 
both of the form h*Vh. This information is contained in the temporal Gauss and Ricci equations: 
they too are of the form h * Vh (the usual spatial varieties look like h * h). 

Another evolution equation we shall need to use on occasion is that of the volume measure. 
This is derived in exactly the same manner as for a hypersurface: 

d , d n 

2^/detg~j dt 
= -y^detgijg'^H ■ hij 
= -\H\^dlJ.g(j,y 

The evolution equations in case where the background space is a sphere will be needed in the next 
chapter, and we delay their derivation until then. 



4.2. Preservation of curvature pinching 

In this section we show that a certain curvature pinching condition is preserved by the mean 
curvature flow. We will often refer to the next lemma as the Pinching Lemma. 

Lemma 4.2. If a solution F : S x [0, T) — t- M"^'^ of the mean curvature flow satisfies + a < 
c\H\^for some constants a < ^ + 4- and a > Oatt = 0, then this remains true for all < t < T. 



44 



4. SUBMANIFOLDS OF EUCLIDEAN SPACE 



Note that under the conditions of Main Theorem |6] (at least in the case where the inequalities 
hold strictly), there exist constants c < ^ and a > such that the conditions of Lemma 14.21 
hold. Thus the result implies both that H remains everywhere non-zero, and that the curvature 
pinching is preserved. Consider now the quantity Q = \h\'^ + a — c\H\'^ , where c and a are positive 
constants. Combining the evolution equations for and \H\^ we get 

(4.5) -^Q = AQ - 2(|V/iP - clV-ffP) + 2Ri - 2cR2. 

at 

By assumption this quantity is initially negative. If there is a first point and time where Q 
becomes zero, then at this point we necessarily have ^ > and AQ < 0. We will derive a 
contradiction by showing that the gradient tems on the right-hand side of equation (14.51 ) are non- 
positive, whilst the reaction terms are strictly negative. We begin by estimating the gradient terms: 

Proposition 4.3. We have the estimates 

(4.6a) |V/iP > -^—\VH\^ 

n + 2 

(4.6b) |V/i|2-i|Vi7p>^i!^|V/i|2. 

n 3n 

Proof. In exactly the same way as MHulll and IIIIa2[ . we decompose the tensor V/i into 
orthogonal components Vj/ijfc = ^ijfc + Fijk, where 

^ijk = — -nidij'^kH + gikVjH + gjkViH). 
n + 2 

Then |V/ip > \E\'^ = ^^|Vi/p. The second estimate follows easily from the first. □ 

Since c < under the assumption of Lemma 14.21 the gradient terms are non-positive. 
In order to estimate the reaction terms of (I4.5t it is convenient to work with the traceless part of 
second fundamental form /i = h—^Hg. The lengths of /i and /i are related by = — 
At a point where Q = 0, we certainly have \H\ ^ 0, so we can choose a local orthonormal frame 
{va '■ 1 < a < A;} for ?sf such that z/i = H/\H\. With this choice of frame the second fundamental 
form takes the form 




and 




At a point we may choose a basis for the tangent space such that hi is diagonal. We denote the 

o o 

diagonal entries of hi and /ii by Aj and Aj respectively. Additionally, we denote the norm of the 
(a 7^ l)-directions of the second fundamental form by that is, = + \h-\'^. We 

also adopt from the following piece of notation from HCdCKl : for a matrix A = (a-ij), we denote 

N{A)=tT {A.A')=Y,ia^,f■ 



o o o o ± 

In particular, we have J2a p ^i^ahp — hpha) = \R\ 
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To estimate the reaction terms we work with the bases described above and separate the (a 
l)-components from the others. The reaction terms of (I4.5t become 



n n 

oi,p i,j a>l i,j 

2 

+ 

a,l3>l i,j 

\R^f = 2Y,N{hiha-Lhi)+ N{hJip-hpK) 



. 2 



a 



Writing out all the reaction terms we now have 

2Ri - 2cR2 = 2Y,{Yhij^hij^y + 2\R\^ - 2c 



, 2 

a,l3 i,j i,j 

(4.7) = 2\h,\^ - 2{c - - -(c - -)\H\' 

n n n 

2 



E/ — ^ o o \ 2 ^ o o o o 

( 2^ hijahijfsj +2 2^ N{hah[s - hpha 



a,/3>l i,i a!,/3>l 



Now we use the fact that Q = to replace (c — \H\^hy |/ip + ain the first line of (I4.7I ). giving 



2\h^\^ - 2(c - -)\hi\''\H\'' - -(c - -)\H\^ 
n n n 



= 2\k\'' - 2\hi\^ [\hi\^ + + a) - + a) (\hi\^ + jLp + a) 

c — l/n n(c — 

where we use the fact that all terms involving a are non-positive, and we have a strictly negative 
term — ^{c-i/n) ■ need to control the last two lines of (14.71 ). In the second last line we proceed 

o 

by expanding the terms and using the fact that hi is diagonal: 

2 ^ — ^ / ^ — ^ o o \ 2 

Mia 

a>l i,j a>l 

\2 



5Z ( =Y{Y1 ^^^^ 

a>l i 



i 

a>l 



3 

a>l 

■ )2 
na I • 
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Also, 



N{hiK - Khi) = Y.^\i- \j?{hi, 

a>l 

Eo o „ o 

(A, - X,f{h, 

< J]2(A,2 + A/)(/;,,„)2 

a>l 

a>l 



SO 



«>1 i,j a>l 



OL>l 

<2\hr\^\h.\^. 

To estimate the last line we use an inequality first derived in HCdCKH for a similar purpose, 
and later improved ULLll to be independent of the codimension. In our notation we have 

E/ ^ — ^ o o \ 2 ^ — ^ o o o o 3 ° ^ 

yZ^hijahiji3j + 2^ N{hahi3 - hpha) < -\h-\ . 

Q,/3>1 i,j a,/3>l 



Proof of Theorem 14.21 Using the above inequalities we estimate the reaction terms by 



2i?i - 2ci?2 < fe - —A—-) \hi?\h-\'' + - — ^— 

\ n[c — \/n)J \ n[c—l/njj 



The terms are nonpositive for c < ;^ + 3^ and the terms are nonpositive for 

c < ^ + The gradient terms are nonpositive for c < so the right-hand side of (14.51 ) is 
negative for c < ^ + while the left-hand side is non-negative. This is a contradiction, so Q 
must remain negative. □ 

To apply the pinching estimate in the case where equality holds in the assumptions of Main 
Theorem [6l we need the following result: 

Proposition 4.4. Suppose T,q = i<o(S") is a submanifold satisfying the conditions of Main The- 
orem\6\and let F : S x [0, T) — )• M"+^ be the solution ofMCF with initial data Fq. Then for any 
sufficiently small f > there exists c < ^ + ^ and a > such that the conditions of Lemma \4.2\ 
hold for Sf. 

Proof. We assume that Sq is not a totally umbillic sphere, since in that case the conditions 
of Lemma 14. 21 certainly apply. Since the solution is smooth, H remains non-zero on a short time 
interval. On this interval we can carry out the proof of Lemma \42\ with a = 0, yielding 



I (|.|^ - S A (|.r^ - 4Hi') - 2 (1 - I V,f + (3 - iLf. 

The coefficients of the last two terms are negative under the assumptions of Main Theorem |6] 
By the strong maximum principle, if — c\H\'^ does not immediately become negative, then 

o 

V/i = and /i_ = 0. The latter implies that lies in a (n + l)-subspace of R"+ , and then 
V/i = imphes that is a product x M""^ C M"+'' (see Chapter 5), and since Sq is not a 
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sphere we have p < n. But this is impossible since T,t is compact. Therefore for any small t > 
there exists a > such that — c\H\'^ < —a on and Lemma l4!2l applies . □ 



4.3. Higher derivative estimates and long time existence 

Here we consider the long time behaviour of MCF and establish the existence of a solution on 
a finite maximal time interval determined by the blowup of the second fundamental form. 

Theorem 4.5. Under the assumptions of Main Theorem^ MCF has a unique solution on a finite 
maximal time interval < t < T < oo. Moreover, maxsj/ip — )• oo as t — )• T. 

As a first step we observe that the maximal time of existence is finite. This follows easily from 
the equation for the position vector F: = A|Fp — 2n. The maximum principle implies 

\F{p, < i?2 _ 2nt and thus T < where R = sup{|Fo(p)| : p G S}. 

Next want to prove interior-in-time higher derivative estimates for the second fundamental 
form. We use Hamilton's * notation: For tensors S and T (that is, sections of bundles constructed 
from "K and !N by taking duals and tensor products) the product S * T denotes any linear combi- 
nation of contractions of S with T. 

Proposition 4.6. The evolution of the m-th covariant derivative ofh is of the form 
VjV™/i = AV"/i + V'h*V^h*V^h. 

i+j+k—m 

Proof. We argue by induction on m. The case m = is given by the evolution equation 
for the second fundamental form. Now suppose that the result holds up to m — 1. Differentiating 
the m-th covariant derivative of h in time and using the timelike Gauss and Ricci equations to 
interchange derivatives we find 

VtV™/i = VVtV"'-^h + V"'-^h *h*Vh 

i+j+k—m—l 
i+j+k=m 

The formula for commuting the Laplacian and gradient of a normal- valued tensor is given by: 

llVkT = Vk^T + Vra{R{dk,dm)T) + {{R{dk,drn){VT)) 

Since T and VT are valued tensors acting on "K, equation (12.41 ) gives expressions for R{dk,dm)T 

as R*T + R*T, and similarly R{dk,dm)VT = R*VT + R* VT, where R and R are the 
curvature tensors on "K and which are both of the form h* h. The terms arising in commuting 
the gradient and Laplacian of V"^~^h are of the form J2i+j+k=m * V^h * V^h, so we obtain 

VtV™/i = AV™/i + Y ^'^ * "^^^ * ^^'^ 

i+j+k='m 

as required. □ 
Proposition 4.7. The evolution o/lV^/ip is of the form 

^IV"'/!^ = AlV^/il^ - 2|V"+i/i|2 + Y V^/i * V-'/i * V^'/i * V^/i. 

i+j+k=m 
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Proof. Denoting by angle brackets the inner product on (g)™+^?f * (g) 7^, which is compatible 
with the connection on the same bundle, we have 

= 2(V™/i,VtV^/i> 

i+j+k=m 

= A|V"^/i|2 - 2|V™+^/i|^ + Yl V*/i * V^7i * V^'/i * V™/i 

i+j+k=m 

as required. □ 

Proposition 4.8. Suppose that mean curvature flow of a given submanifold Sq has a solution on 
a time interval t G [0, r]. // < K for all t G [0, r], then |V™/ip < (1 + l/f^)for all 
t G (0, r], where Cm is a constant that depends on m, n and K. 

The strength of this estimate is that assuming only a bound on the second fundamental form 
(and no information about its derivatives) we can bound all higher derivatives. The fact that these 
estimates blow up as t approaches zero poses no difficulty, since the short time existence result 
bounds all derivatives of h for a short time. While not crucial here, the interior-in-time estimates 
are useful in singularity analysis. 

Proof. The proof is by induction on m. We first prove the Lemma for m = 1. We consider 
the quantity G = t] V/ip + \h\'^, which has a bound at t = depending only on curvature. The 
strategy is now to use the good term from the evolution of to control the bad term in the 
evolution of |V/ip: Differentiating G we get 

8G 

— = |v/i|2 + t{A\Vh\^ - 2|V2/i|2 + /i * /i * V/i * V/i) 

+ (A|/ip - 2|V/ip + h*h*h*h) 
<AG+ {cit\hf - l)|V/ip + C2\h\\ 
For t < l/{ciK) we can estimate 

at 

and the maximum principle implies max^^ ^ G < K + C2K'^t. Then | V/ip < G/t < K/t + C2K^ 
fort G (0,l/(ciK)]. Ift > l/{ciK) we apply the same argument on the interval [t~l/{ciK),t], 
yielding |V/ip(t) < (ci + 02)-?^^. This completes the proof for m = 1. Now suppose the estimate 
holds up to m - 1, and consider G = t"'\V'^h\'^ + mt'""^|V"'~^/ip. Differentiating G gives 

i+j+k=m 

+ m{(m- l)t™-2|V™-i/i|2 + r-i(A|V™-i/i|2 -2|V™/i|2 

i+j+k=m—l 

Noticing that in the quartic reaction terms there can only be one or two occurences of the highest 
order derivative, using Young's inequaUty we can estimate 

^ /~i ^ j.m— 1 ii-rm 7_ 1 2 , j.m f a i^?Tt7_i2 , _ i^mr_i2 , ''I 



-G < mt'^-'\V"'h\' + AlV^/il" + C3|V™/i| 
ot I t 



+ m{(m - l)i'"-2|V"*-^/ip + t'"-^(A|V'"-^/i|2 - 2|V"*/i|2 + c5\V"'-^hf + ^)}- 
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We split the gradient term of order m out of the second line, and then since m is at least two, all 
other terms are bounded by the induction hypothesis for t < 1, giving 

< AG + (est - m)i™- + C7. 

Thus -^G < AG + C8 if t < min{l, m/ca}, so by the maximum principle | V™/ip < C/i™ for 
t < min{ 1, m/cs}. The same argument on later time intervals gives the result for larger t. □ 

Proof of Theorem 14. 5 1 Fix a smooth metric ^ on S with Levi-Civita connection V. g 
extends to a time-independent metric on IK, and V extends to "K by taking Vq^u = whenever 
[dt,u\ = 0. The difference T = V - V restricts to a section of Ji* (g) Ji* (g) Ji. If S is a 
section of a bundle constructed from "K, J{ and F*TN, VS denotes the derivative of S with the 
connection on this bundle induced by the connections V on IK, V on 3sf, and V on F*TN, so 
that VS -VS = S *T. 

To prove Theorem 14. 5 1 we assume that \h\ remains bounded on the interval [0, T), and derive 
a contradiction. This suffices to prove the Theorem, since if |/i| is bounded on any subsequence 
of times approaching T, then Equation (14.11 ) implies that \h\ is bounded on S x [0, T). Under this 
assumption the boundedness of Vtg = —2H • h implies that the metric g remains comparable to 
g: We have for any non-zero vector v e TS 











9{v, 


v) 


dt 






g{v,v) 


9{v, 


v) 



< 2\H\ 



9{v,v) 
'g{v,v) 



so that the ratio of lengths is controlled above and below by exponential functions of time, and 
hence since the time interval is bounded, there exists a positive constant cg such that 



(4.8) 



1 

eg' 



-g < 9 < egg. 



Next we observe that covariant derivatives of all orders of F with respect to V can be expressed 
in terms of h and T and their derivatives: We prove by induction that 

(4.9) 



F^V'^-^T + F, 




V^F 



+ {i + F^ 
This is true for k = 2, since 



V 



^k-3 



,X;(m+l)Pm=i "1=0 



(4.10) 



F 



To deduce the result for higher k by induction, we note that equation (14.101 ) implies a formula for 
the derivative of F^, : 

(VF,)(y) = F*r(., V) + ih{., V) = F^T*V + Lh* V, 
while equation (12.211 ) gives 

(VO(e) = -F,W{.,C) = F,h*^. 

The result for A; + 1 now follows by differentiating the expression (I4.9l l. and writing V(V"/i) = 
V"+i/i + V"/i * T. It follows that if \V^F\g is bounded for j = 1, . . . , /c - 1, then 



(4.11) 
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The above observations allow us to prove C'^ convergence of F as t — )• T for every k: We have 
VtF = iH, so the boundedness of H implies that F remains bounded and converges uniformly 
as t — r. Differentiating as above, we find by induction that 
(4.12) 

fc-l / fc-2-j . \ / j+1 

ViV'^F = (F,+0* n (^"^)' )* ^ WiV^hf- 

i=0 \E("+l)in=A;-l-j n=0 / \£{m+l)pm=j+2 m=0 

Suppose we have established a bound on \ F\g for j < k — 1. Then using the estimate (14. Ill ), the 
bounds on |V"/i|g from Lemma l4~8l and the comparability of g and g from (I4.8l l we can estimate 

|VtV'=F|^ <c(l + |V*-'"2r|g) <c(l + |V^'F| 

so that |V''F| g remains bounded, and V'^F converges uniformly as t — T. This completes the 
induction, proving that F{.,t) converges in C°° to a limit T) which is an immersion. 

Finally, applying the short time existence result with initial data F{.,T), we deduce that the 
solution can be continued to a larger time interval, contradicting the maximality of T. This com- 
pletes the proof of Theorem 14.51 □ 

4.4. A pinching estimate for the traceless second fundamental form 

In this section we show that the pinching actually improves along the flow. This is the key 
estimate that will imply that the submanifold is evolving to a "round" point. 

Theorem 4.9. Under the assumptions of Main Theorem^there exist constants Cq < oo and 5 > 
both depending only on Sq such that for all time t G [0,T) we have the estimate 

(4.13) I^P < Co\H\^-^. 

We wish to bound the function /o- = - l/n|i7p)/|_ffp(^^'^) for sufficiently small a. As 
in the hypersurface case, a distinguishing feature of mean curvature flow when compared to Ricci 
flow is that this result cannot be proved by a maximum principle argument alone. Somewhat more 
technical integral estimates and a Stampacchia iteration procedure are required. We proceed by 
first deriving an evolution equation for fa-. 

Proposition 4.10. For any a G [0, 1/2] we have the evolution equation 

(4-14) < A/. + ^-^{V,\H\,V,fa) - |^p(T,,) |ViJ|^ + 2a\h\^fa. 

Proof. Differentiating in time and substituting in the evolutions equations for the squared 
lengths of the second fundamental form and mean curvature we get 

A|/i|2 - 2|V/i|2 + 2Ri 1 (Al/^l^ - 2|VF|2 + 2R2) 



dtfa 



(4.15) 



:^-'')<l"l'-'/"l«fl(A|if|^-2|V//P+2fl.). 



2\2-a 



The Laplacian of is given by 

(1 - a){\h\' - yn\H\') , ^ (2-a)(l-a)(|/^|2-i/n|g|2) 



Using this and the identity 

-fa\Hf\V\H 
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equation (14.151) can be manipulated into the form 

2(1 - a) , ^ ^ , 2 ^,^^,2 \h\^ ,^^,2^ , 2alhfa 



Ml^f lFf|2|Y7lFf||2 2a(|/.p-l/n|gp) 2 
2 / 

We discard the terms on the last two lines as these are non-positive under our pinching assumption. 
The gradient terms on the first line may be estimated as follows: 



(\H\^)i-<ry' ' |i/|2'' V- (|if|2)i-^ Vn + 2 

and also R2 < Importantly, observe that if c < 4/(3n), then ey := 3/(n + 2) — c is 

strictly positive. □ 

The small reaction term 2a\h\'^ f^r in this evolution equation is positive and hence we cannot 
apply the maximum principle. As in the hypersurface case, we exploit the negative term involving 
the gradient of the mean curvature by integrating a suitable form of Simons' identity: Contracting 
equation (12.241) with the second fundamental form we obtain 

(4.16) ^A|^|2 = hij ■ V^VjH + |V^|2 + Z, 

where 

2 



a,P i,j i,j,p 



Lemma 4.11, T/'S" is a submanifold ofW"'^^ that satisfies H and < c|//p, where 

f<^, n = 2,3 

then there exists ez > such that Z > ez\h\'^\H\'^. 

The example given in the Introduction shows the best value of c that can be expected is l/(n — 
1). In dimensions greater than four, 4/(3n) > l/(n — 1) and so somewhere in the analysis 
the condition c < l/(n — 1) had to manifest itself. For a submanifold of Euclidean space, the 
condition |/i|<l/(n — 1) implies that the submanifold has positive intrinsic curvature. Just as in 
the hypersurface case (where strict convexity implies positive intrinsic curvature), it is the positive 
intrinsic curvature that makes this lemma, and indeed the Main Theorem true. 

Proof of Lemma |4. Ill Working with the local orthonormal frames of Section I4.2l we ex- 
pand Z to get 

Z = -\hi\'^ + + - 2 JJ^J] kki,^'^ -2^N{hiK - Khi) 

E/ ^ — ^ o o \ 2 — ^ 00 00 

[y^hijghijp] - 2^ N{hahp -hpha) 

a,/3>l i,j a,/3>l 



a>l a>l 
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We estimate the first summation term on line one and the two terms on line two as before, namely 

2 



a>l i,j a>l 

i 

X — ^ /x — > ° ° \ ^ X — > o o o o 3 ° /I 

- Z^[Z^hijahiji3j - 2^ N{hahp-hpha)>--^\h-\ , 

a,/3>l i,j a,0>l 

however we need to work somewhat harder with the remaining summation terms. 
Proposition 4.12. For any ij > Owe have the following estimate 



- 2 J2 ^(hiL - Khi) + Y,\H\\ii}iijaf 



a>l a>l 



> -max {4, l}\k\\\h-? - Y.(huaf) - f - Y.^kr.f)- 



An 

a>l ' a>l 



Proof. Using the Peter-Paul inequality we estimate 

- 2 ^ N{hiha - Khi) + Y}H\\i{hi 



H]a) 

a>l a>l 



= -Y, {2(A, - Xjf - ^(A. + A,)}fe)2 

>-Y^ {2(A, - A,f + ^(A, + A,f - ^|//p(|L|2 - ^(La)') 

a>\ ' a>l 

= - E {(2 + |)(A^' + A/) + (| - 4)A.A,}(/;,,«)2 

a>l 



' a>l 



If ?7 > 8 we estimate 

^ - E {(2 + |)(A^' + A/) + - 2)(A,2 + A/)}fe)2 

a>l 

' a>l 

i 

a>l ' a>l 

i^j i 

a>l ' a>l 
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while for 7/ < 8 we can similarly estimate 

^ - E {(2 + + A/) + (2 - ^)(A.2 + A/)}(/;.,„)2 

Q>1 



Q>1 



a>l ' a>l 

i^j i 

> -Mk\'{\h-\' - Y^Chuaf) - ^\H\'{\L\' - Y^Chua?). 



□ 

To estimate the remaining two terms we use the following two inequalities from HAdCI and 
EH: 

^ ^yn{n - 1) 

Y.\H\xSuaf > --^=l=\H\\k\ J2ih 

i i 

and further estimate them using the Peter-Paul inequality to obtain 

^2 



2 

iia j 1 



^ 2' ' 2/in(n-l)' ' ' ' 

i i i 

Note that in dimension n = 2 above two terms are actually zero and there is no need to further 
estimate them in this way. For n = 2 the remaining quantities can now be estimated as we have 
done before to give the estimate for c < 3/4. For the higher dimensions, putting everything 
together we obtain 

Z>-\k\' + -\k\^\H\^ + - 2\hl\' Tihuaf - Ifh^f 

n n ^-^ I 

i 

- max{4, + max{4, ']^\h^\^ ^fe)' - ^\}if\hA^ 

i 

i i 
2 



4p nin — 1) 



2 

ua ) • 



We now need to choose the optimal values of the constants rj, p and p. First, choose p to be 
equal to n — 2 and p = {n — 2)/2. Next we want to choose tj to make the ^ {huaf' terms 

i 

non-negative, that is, we would like to choose 77 so that 



max|4, II - 2 - p > 0. 
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As /9 = (n — 2)/2, we want 

( , rj) n + 2 

max < 4, — > > , 

I ' 2i - 2 ' 

thus we are able to choose r/ = n + 2. In dimensions 3 to 5 this term is positive and we discard it, 
while forn > 8 it is identically zero. The only mildly troublesome term that remains is 

1 1 (n-2)2 



4r/ 4/9 n{n — 1) 



iia ) • 



With our choices of t] and p this term this term is negative for n > 3 and we estimate 

2nn-l 4n + 2W' ' Z^^ - V2n(n-1) 4(n + 2 W ' ' 

i 

After substituting in our choices for p, p and rj we have, in dimension three to five: 

z>-\h,\' + -l/^iplFp + l|Lp|F|2 - "^M' - 4|/.ip|Lp 

n n 2 

2' ' 2n(n-l)' ' ' ' 2n(n - 1) ' ' ' ' ' 
and in dimensions six and higher: 

,° .A 1,° ,o, ,o 1,° ,o, ,o ?^ — 2° ,^ 71 + 2 , o ,o,° ,o 
- I il -r ^1 ill I -r ^1 III 2 ' ' 2 ' ' ' ' 

2' ' 2n(n - 1)' ' ' ' 2n(n - 1) ' ' ' ' 

We now group like terms, estimate from below by |/ip/(c — 1/n) and calculate the maximum 
value of c permissable in each case such that the coefficients are all strictly positive. For n = 2 and 
n = 3 the most restrictive term is the cross-term, and the best value of c is given by c = 3/4 and 
c = 11/24 respectively. For the corresponding value of n note that both 3/4 and 11/24 > 4/(3n) 
and so we have simply used 4/(3n) in the statement of the lemma. For n > 4 the most restrictive 
term is the \hi\ term, which is identically zero when c = l/(n — 1). Thus, for the values c 
stated in the proposition, we have now shown there exist strictly positive constants C2, C3 and C4 
depending on Sq such that 

Z > C2\hi\^ + C3|^i|2|L|2 + Ci\h-f 
(4.17) > cS\\ 

where C5 = min{c2, 03/2,04}. To prove the desired estimate we note that by using Peter-Paul on 
various terms of Z we can estimate 

Z > C6\h\^\H\'^ -cjlhl"^. 
Combining this with (14.171 ) gives for any a G [0, 1] that 

Z > a{c6\h\^\H\^ - C7\h\^) + (1 - a)c5\h\^. 
Choosing a = 05/(05 + C7) gives 

~ C5 + C7 ' ' ' 

and the lemma is complete by setting ez = 0506/(05 + 07). □ 

Next we derive the integral estimates. 
Proposition 4.13. For any p >2 and rj > Owe have the estimate 
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Proof. Using the contracted fomi of Simons ' indentity and A|i7p = 2\H\A\H\+2\V\H\\^, 
the Laplacian of can be expressed as 

2 1,^„,2\ 2(l-a)(l-2a)^ ,^,^,,2 



9 



The combination of the last two terms is non-negative and we discard them. We muhiply the 
remaining terms by f§^^ and integrate over S. On the left, and in the last term on line one we use 
Green's first identity, and in integrating the first term on the right we use the Divergence Theorem 
and the Codazzi equation. The term arising from integrating on the left is non-negative and we 
discard it. Two other terms arising from the integration combine, ultimately giving 

-4(1--) j^^j^^^,{^m-k,v,H)d,, + ^-^^ /^^fr^iv^l'rf^, 

- 2(1 - a){p - 2) ^ !^(y,\HlV,f^) dfig + 2(1 - a) ^\V\Hf dfig. 

Note the terms with an inner product do not have a sign. Using the Cauchy-Schwarz and Young 
inequalities, the inequalities < c|iJp'^, |V|i?|p < |Vi7p, 1 - a < 1, c < 1, and = 
fa\H\'^^^~'^^ we estimate each term as follows: 

<^//r'|V/.|^d/., + (p-l)r? l^jJ^\^H\'d^,,; 
-2{l-a){p-2) I lL^{V,\H\,VU)df,g 

"'SI I 

<^/^/r'|V/.pd/., + (p-2)^ I^^Jl^lVHfd^g; 



fa f fS~^ 

-^|V|i/||2 dfig < 2 / |^p(i_,) I V - I ^^^g■ 



2(1-'^) J^T^MHW'df,g<2 I ^^Z,_j VH\'dfig. 
Putting all the estimates together we obtain 



Our use for this inequality will be to show that sufficiently high norms of fa are bounded. We 
are not interested in finding optimal values of p and consequently we are going to be a httle rough 
with the final estimates in order to put the lemma into a convenient form. Setting n = r], and using 
p — 2<p— l<p and Lemma l4. lll we get 

2ez I fl\H\^<{2pr^ + ^) [ J^i^l VFp d/., + / f^'l^ f.\' d^^g. 

Dividing through by 2ez completes the Lemma. □ 
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Proposition 4.14. For any p > max{2, 8/ (ey + 1)} we have the estimate 



Proof. Differentiating under the integral sign and substituting in the evolution equations for 
ffj and the measure d^g gives 

<-p{p-l) j fP-^\VU\^di^g + 4il-a)p f lL^\V\H\\\VU\di,g 
(4.18) ' ' 



We estimate the second integral by 



A{l-a)p jJ-^MH\\\VU\d^Xg 

<J |^fr'N^\'d^^g + 2pp |^^Jl^\VH\'d^,g, 



and then substituting this estimate back into (15.131 ) gives 
d 



, /? < ( - MP - 1) + ^) fS-'NUf d,, 

We now want to choose p so that 1 — 2/{p{p— 1)) > 1/2 and p so that 1 — p/ey > 1/2. Choosing 
p = 4/(p — 1) and p > max{2, 8/(ev + 1)} gives the result. □ 

Lemma 4.15. There exist constants cg and cg depending only on Sq such that if p > cg and 
<7 < cg/ y/p, then for all time t £ [0,T) we have the estimate 



1 



where Ci is a uniform constant. 

Proof. Combining Propositions 15.51 and 15 . 6 1 we get 



Pev / ^^^7TT^|Vif| dpg 



Suppose that 
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Set ri = i:a/ez, then 



( 2£z ^ 1 
' er; 2 



1 



For the second last inequahty to hold we must assume p > 64/ev. We conclude that 

d 



dt 



This implies the lemma with cs = max{2,8/(ev + l)j 64/ev}, cg = e^A/ev/S and Ci = 
(iSo + l|)max^g [0,1/2] (maxso /^). □ 

An important corollary of this lemma is the following, which states that for larger values of p 
and smaller values of a, powers of H can be absorbed into f^^. This property is key in the final 
iteration argument. 

Corollary 4.16. Forp > max{c7, 4n^Cg} and a < {c%/2) / the estimate 

HV^df^g< J fl,d^Jig. 

s s 

holds ont£ [0, T). 

Proof. We need a' = a + n/p < C'g,! ^ for sufficiently large p and small a . Suppose that 

p > max{c7, 4n^/c8} and a < {cs/2)/^). Then 

/ n C8 I n Cs 
a = a -\ — < 1 < — 

as required. □ 

Lemma 14.151 shows that sufficiently high LP norms of are bounded. We now proceed to 
derive the desired sup bound on /^r by a Stampacchia iteration argument. The argument rests on 
the following well-known iteration lemma. 

Lemma 4.17. Let ip{t), ko < t < oo, be a non-negative and non-increasing function which 
satisfies 

for h > k > kQ, where C, a, and (3 are positive constants with /3 > 1. Then 

ifiko + d)=0, 

where 

d° = C7|(^(A:o)|'^-^2°'^/('5-i). 

For a proof of this lemma we refer to reader to [KSl. Continuing with the iteration argu- 
ment, set ko := maXg.g[o,i/2] ™axso fa- For any k > k^, define the truncated function f^^k '■= 
max{(/o- — k),0} and the set A{k,t) := {x G Sj : f^^k > 0}. In exactly the same manner as 
Proposition 15. 5 1 we derive the following evolution equation for f^^k- 

d f rp ^ / P{P-^) f .P-2|v7^ |2 , f fl,k 



dt 



I fl,u di^, < [ /^^ I V/.,.P d^., - pesj [ jji^Wm' d^, 

JAik.t) ^ JA(k,t) JA{k.t) Li^^"^ 

+ 2pa [ m'^fl.d^ig. 

JAik.t) 



For p > 8 we estimate 
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then setting = /^^ and discarding the second temi on the right we get 

(4.19) / Vldflg+ [ \VVk\^dflg<2pa [ \H\yP^ ^d^lg. 

JA{k,t) JA{k,t) JA{k,t) 

We now make us of the Michael-Simon Sobolev inequality MMSB . which states that for any function 

u G Co'\s) we have 

n-l 



S 



u\'^-^ dfig) " <Cs / (\Vu\ + \H\u) d^g, 



where Cs (the Sobolev constant) is a constant that depends only on n. The Michael-Simon 
Sobolev inequality is a generalisation of the standard Sobolev inequality to functions on a sub- 
manifold. The form of the Michael-Simon Sobolev inequality stated above corresponds to the 
p = 1 case of the standard Gagliardo-Nirenberg-Sobolov inequality. To obtain the inequality in 

the case 1 < p < n we set t> := |n|''', where 7 = p{n — l)/{n — p) > 0, and after a use of Holder's 
inequality we find 

jWfdy.g)^<Csl[j^Vu\Pdfig)'^+Cs[jm''d^lg)''[jH''*d^^^^ 

We want to take advantage of the good gradient term on the left of ( 14.191) . and so we need the 
Sobolev inequality with p = 2. Squaring both sides, using (a + 6)^ < 2(a^ + 6^) and then setting 
g = n/(n — 2) if n > 2 or any number finite number if n = 2, we obtain 

(4.20) [j^vlUflgY <Cl0[jjVVk\^dflg^ + Cli(^ jjH\^ dlXg^^' j^vlUlXgY . 

Using that by definition is zero outside of A{k,t) and CoroUarv I4.16l we estimate 

fP X 2/n 



JAiKt) ' " ^JA{Kt) ' 

^JA{k,t) 



2/n 



(4.21) 

^JA{k,t) ' 

Therefore we can fix a /ci > /cq sufficiently large such that for all k > ki the second term on the 
right of (14.201 ) can be absorbed into the left giving 

(4.22) [J^vlUfigY <ci2[ljVvk\^dfXgy 
Combining equations (14.191 ) and (14.221 ) we obtain 

(4.23) ^ / vldfig+( [ vlU^igX'" < C13 / m^f^dflg. 

"I JA{k,t) ^ JA{k,t) ' JA{k,t) 

Integrating this equation from t = until some time r G [0, T] we get 

[ vldfig- [ vld^g+ [ ( [ vl'^dfig) ^ dt < Ci^i [ [ \H\'^ fP dfig dt. 

JT,r J^o JO ^JA{k,t) ' Jo J A(k,t) 

By definition of v^, the integral evaluated at the initial time is zero. Denote by the time when 
the first integral on the left achieves its supremum, that is = sup^gjo.r] /s^ ^fc d^J-g- We integrate 
(14.231) until and until T and add the two inequalites, then discarding two unwanted terms on the 
left and estimating by T on the right we obtain 

(4.24) / vldfig+r([ vl''df,gy^'dt<cn r [ \H\^fPdf,gdt. 
JT.t^ Jo ^ J A{k,t) ' Jo JA{k,t) 
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We now need to estimate the remaining two integrals on the left. Recalling the standard interpola- 
tion inequality for W spaces: 

II II ^ II ll^ll l|i-^ 
Irllgo — Irllilrllij i 

where 1 < < q and < < 1, we interpolate with 6* = 1/go to get 

Using the above interpolation inequality, and the Holder and Young inequalities we see 



7 

J A(k,t) ' 

( [ ( [ vldiigX" ( / vl^d^xg) '''dt) 

^JO ^JA{k,t) ' ^JA{k,t) ' ' 



< 



< sup ( [ vl d^g) ■ ( [ if Wfc" dng\ '"dt" 
t&[o,T] ^ JA{k,t) ' ^Jq ^JA{k,t) ' 

cT 



< sup / vld^^g + qoU-^) V (/ vl'^ dfig) ^''dt 
te[o,T]JA{k,t) ^qo-i-'' Jo ^JA{k,t) ' 

< sup / vld^g+ f [ I d^g) ^ dt, 

te\0,T] JAik.t) Jo ^JA(kd) ' 



and using this to estimate (14.241 ) from below we obtain 

( r / dlXgdt)'"" < C13 r [ \H\^fPdlIgdt. 

^Jo J A(k,t) ' Jo J A{k,t) 

Set \\A{k^ II = X4(fc d^Xg dt. Estimating the left from below by the Holder inequality gives 



l/r 



vl!'' dfigdty^'' >( [ [ vld^igdt) . P(A;,t)|ho \ 
JA{k,t) ^ ^Jo J A(k,t) ^ 

and the right from above with Holder's inequality 

r [ Vldflgdt < Ci3P(A:,t)||2-l/.0-l/r/' r f \H\^rfpr^^^^^ 
Jo J A{k,t) ^Jo J A{k,t) 

Choose r sufficiently large so that 7 := 2 — \/qo — l/r > 1, and then by the same argument as 
(I4.21I ). the second factor on the right can be bounded by a constant. For h > k > ki we estimate 
the left from below as follows: 

[ vldfigdt>\h-k\P\\A{h,t)\\, 

J A{k,t) 

finally obtaining 

\h-knA{hM<ciA\\A{k,tW 
which again holds for all /i > A; > fci. From Lemma 14.171 and the defintion of A{k, t) it follows 
that fa <ki + d, where dP = C52P^/^^+'^^\\A{ki,t)p-'^. Since f^(^k^,t) ^ 1^*1 ^ l^ol and 
the maximal time of existence T is finite, we conclude that fa- < Co, where Co is positive uniform 
constant, and the theorem follows. 



4.5. A gradient estimate for the mean curvature 



In this section we derive a gradient estimate for the mean curvature. This will be used in the 
following section to compare the mean curvature of the submanifold at different points. 
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Theorem 4.18. Under the assumptions of Main Theorem^ for each r] > there exists a constant 
Cn depending only on r] and Sg such that the estimate 

(4.25) IVH]"^ <Tj\H\^ + Cr, 
holds onT, X [0, T). 

We begin by deriving a number of evolution equations. 

Proposition 4.19. There exists a constant A depending only on Sq such that 
8 

— iViJp < AlViJP - 2|v2i?|2 + ^|i?|2|v/iP. 
dt 

Proof. Differentiating jV-fTp in time gives 

||Vi^P = |(VF,Vi^> 

= 2 {VtVH, VH) 

= 2g{VkVtH + R{dk, dt)H, VkH) 

(4.26) =2g{Vk{AH + H ■ hpghpg),VkH) + 2g{R{dk,dt)H,VkH). 
To manipulate the last Line into the desired form we need the following two formulae: 

A|i7p = 2g{AVkH, VkH) + 2\V'^H\'^ 

AVkH = VkAH + Vp{R{dk, dp)H) + R{dk, dp)VpH + RcpkVpH. 

Substituting these into (14.261 ) and observing that the Gauss equation (I2.14al ) and the Ricci equation 
(12.161) are of the form R = h* h and R = h * h, and that the timelike Ricci equation (12.171 ) is of 
the form R{-, dt) = h* Vh, we find 

-^iViJP = AlViJP - 2|V^i?P + h*h*Vh*Vh. 
dt 

The proposition now follows from the Cauchy-Schwarz inequality and the Pinching Lemma. □ 
Proposition 4.20. For any Ni, N2 > Owe have the estimates 

(4.27) ^\H\^ > Ali^P - 12|iJp|ViJp + -\Hf 
at n 

^((iVi + iV2|i^|')|/i|') < A((iVi + N2\H\^)\h\'') - ^^"l" ^\ n2 - l)\H\^\Vh\^ 

(4.28) . _ . '''' 

- 3^ ^ iVi - ci(iV2))|V/i|2 + C2{N^,N2)\h\WH\'' + 1), 

where ci and C2 depend only on Sq, Ni and N2. 

Proof. The evolution equation for |//|^ is easily derived from that of 

d_ 

dt' 

Equation (l5?T6l ) follows from the use of |V|iJ|p < \VH\'^ and R2 > l/n\H\^. To prove (IST/l ). 
from the evolution equations for and \H\'^ we derive 



\H\'^ = A\H\'^ - 2|V|i?|2|2 - 4|iJp|Vi7p + 4i?2|^f|^- 



^^{{N^ + N2\H\')\h\') 

= A{{Ni + iV2|F|2)|^|2) - 2iV2(Vi|i7|2, Vi\h\^) - 2N2\h\'^\Vh\'^ + 2iV2i?2|^|^ 

- 2{Ni + 7V2|^|^)(|V/i|2 - -\VH\'^) + 2{Ni + N2\H\'^){Ri - -i?2). 

n n 
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We estimate the second teiTti on the right as follows: 

-2N2{Vi\H\^,Vi\h\^) < 8iV2|/i||^||VF||V/i| 

< 8N2\H\^/^\Vh\'^y^\H\^-^/^ 

in 

o 

Using Young's inequality, R2 < \h\'^\H\'^, and i?i — l/ni?2 < 2|/ip|/ip we estimate 

2N2R2\h\'^ + 2{Ni + N2\H\'^)iRi - -R2) < C2iNi, N2)\h\WH\^ + 1), 

n 

and equation (15.171) now follows. □ 

Proof of Theorem |4J8J Consider/ := \VH\'^ + {Ni + N2\H\'^)\h\'^. From the evolution 
equations derived above we see / satisfies 

|/ < A/ + A\H\'\Vh\' - ^^^iN2 - m\'m' - ^^^(^1 - ciiN2))\Vh\' 

+ C2{Ni,N2)\h\^{\H\^ + l). 

Choose large enough to consume the positive tern arising from the evolution equation for 
|Vi/p. This leaves 

|/ < A/ - ^^^{N2 - mfm' - ^^^{Ni - c,{N2mh\' 

+ C2{Ni,N2)\h\WH\'' + l). 
Now consider g := f — r]\H\'^. From the above evolution equations we have 

< - ^^^{N2 - m\'m' - ^^^{N^ - c,{N2))m' 

+ C2{NuN2)\h\WH\^ + 1)) + 12r/|/7|2|VF|2 - ^\Hf. 

n 

By choosing sufficiently large the gradient term on the last line can be absorbed, and then we 
choose A'^i larger again to make the |V/ip term negative. We finally discard the negative gradient 
terms to get 

^5 < + C2{Nr,N2)\h\WH\'' + 1) - ^\H\\ 
ot n 

Using Theorem I4.9l and Young's inequality we further estimate 

g^g < A5 + C3, 

from which we conclude g < c^. The gradient estimate now follows from the definition of g. □ 

4.6. Contraction to a point 

In Section l43] we established that MCF has a unique solution on a finite maximal time interval 
< t < T determined by the blowup of the second fundamental form. With the results of the 
previous two sections in place, we can now show that the diameter of the submanifold approaches 
zero as t — )• T, or put another away, the submanifold is shrinking to a point. This combined with 
Theorem I4.5l then completes the first part of the Main Theorem |6l 

Theorem 4.21. Under the conditions of Main Theorem^ diam — )• a^' t — )• T. 

The proof is an adaption of Hamilton's use of Myer's Theorem in Section 15 of IIHa2ll . however 
here our pinching condition gives a strictly positive lower bound on the sectional curvature of St 
and we can use Bonnet's Theorem instead. A proof of Bonnet's Theorem can be found in many 
places, for example [P|. 
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Theorem 4.22. Let M be a complete Riemannian manifold and suppose that x £ M such that the 
sectional curvature satisfies K > K^m > along all geodesies of length Tr/^/K^from x. Then 
M is compact and diam M < tt/^/K^. 

We will also need the following result due to Bang- Yen Chen: 

Proposition 4.23. For n > 2, if S" is a submanifold of M"'''^, then at each point p ^ TP the 
smallest sectional curvature K^m satisfies 

K^Up)>l{:^\H\{pf-\h\{pf). 

The proof is a consequence of careful estimation of terms appearing in the Gauss equation and 
can be found in [Cj Lemma 3.2] Combining this with our pinching assumption we see 

(4.29) K,^M > - c) \H\{pf = e'l^l(p)' > 0. 

2 Vn — 1 / 

Lemma 4.24. The ratio \H\,„ax/\H\min 1 as t ^ T. 

Proof. From Theorem l4.18l we know that for each rj > there exists a constant Cn such that 
\VH\ < rj\H\'^ + Cri on < t < T. Since |i^|max — )• oo as t — )• T, there exists a T{r]) such that 
< l/2r?|i7|^^ for all r < t < T, and so \ VH\ < ■n\H\'^^^ for all t > r. For any a G (0, 1) 

choose r] = ""(^~°")^ _ Let t e [T{r]),T), and let x be a point with \H\{x) = |^^|max- Then along 
any geodesic of length from x we have \H\ > |if|max - i^^^^??l-f^lmax = o-|^^|max, and 

consequently the sectional curvatures satisfy K > 

^Imax- The Bonnet Theorem applies to 
prove that diamM < — , so that |i/|min > c|-f^|max on the entire submanifold St for all 
[T(r?),r). " □ 

Since |i/|max — )• oo as t — )• T, the last lemma show that the same is also true for |-ff|mm- 
Bonnet's Theorem now implies that diam — > as t — )• T, which completes the proof of the 
first part of Main Theorem |6] 



4.7. The normalised flow and convergence to the sphere 

The second part of the Main Theorem |6] deals with the asymptotic shape of the evolving 
submanifold as t — )• T. Here we shall show that a suitably normalised flow exists for all time and 
that the (normalised) submanifold converges to a sphere as time approaches infinity. This clarifies 
the sense in which un-normalised submanifold shrinks to a 'round' point. 

We denote quantities pertaining to the normalised flow by a tilde. We are going to define the 
normalised flow in such a way so that the size of the area of the evolving submanifold remains 
constant. We do this by multiplying the solution of MCF at each time by a positive constant 'ip{t) 
so that the measure of the the normalised submanifold is equal to the measure of the initial 
submanifold Sq: 

Fi;t) = mn,t) 

such that 

(4.30) / d^lg^t) = |So|. 

The above rescaling is so far only a rescahng in space. 
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Proposition 4.25. Suppose we rescale an immersion F by F := ip ■ F, where ip is a positive 
constant. Then various geometric quantities rescale as follows: 

9 = i^^g 

h = tph 
H = tP~^H 
\h\' = ^p-'\h\' 

dp-git) = 

V = V 
A = V"^A. 

Proof. Beginning with the metric, we have 

gij = {F^di,F^dj) 

= i?{F,di,F,dj) 
= ip'^gij- 

We also have = ijj~'^g~^. For the second fundamental form, using Gauss' formula and noting 
i^a = we have 



which is just hij = tphij. The mean curvature follows from the inverse metric and the second 
fundamental form: 

H = = ^-''H. 

Similarly for the length squared of the second fundamental form we get 

For the measure we have 



dfig{t) = V^et gij = ^"V'det^ij = Vd^J.g(t). 
The Christoffel symbols are scale-invariant: 

^% = \g^\digjk + djgik - dkgij) 

= \'4^~'^g^\tp'^digjk + ip'^djgik - ip'^dkgij) 

and thus so too is the coimection. Finally, the Laplacian is given by 



□ 
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We now derive the evolution equation for the normalised flow with respect to the time variable 
t. Differentiating (14.301 ) with respect to t we have 



which implies 



dt n J diig{t) 

where the last line follows because the rescaling factor is a function of time and the integration 
is over spatial variables. Define the average of the squared length of the mean curvature over the 
submanifold by 

^.^ !Wd^igi^t) 

fdf^git) 

The evolution equation for the normalised flow with respect to the time variable t is now given by 

dF di^^ , dF 

We now rescale in time to divide out the factor of tp'^ in the above equation. Note that from now 
on a tilde represents a rescaling in both space and time, and not only a rescaling in space as was 
previously the case. We define the rescaled time variable by 

i{t) := l%\T)dT, 

and so di/dt = ijj'^. We now have 







dF , _2 dF 



dt " dt 

= H + -hF, 

n 

where this normalised flow is now defined on the time interval < i < T. Next we want to 
show how various estimates and evolution equations for the normalised flow can be obtained from 
their un-normalised counterparts. The scahng-invariant estimates are the easiest, since the the 
normalising factor Tp simply cancels from both sides and the same estimates hold: 

Proposition 4.26. The following estimates hold for the normalised flow: 

(4.31a) < c\H\^ 

(4.31b) LJmin _^iasi-^f 

\H\^ 

(4.31c) K^in > e^\H\^ 



The following lemma shows how evolution equations for the normalised flow can be easily 
obtained from their un-normalised counterparts: 
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Lemma 4.27. Suppose that P and Q depend on g and h, and that P satisfies the ( un-normalised) 
evolution equation dP/dt = AP + Q. If P has "degree" a, thatis, P = ip^P, then Q has degree 
{a — 2) and P satisfies the normalised evolution equation 

dP ~ ~ ~ a~ ~ 
^ = AP + Q + -hP 
ot n 

For a proof of this lemma see Lemma 17.1 of IIHa2l and Lemma 9.1 of BHulL The evolution 
equation for the metric does not follow from this lemma, but is easily derived in the same way as 
the un-normalised equation. 

Proposition 4.28. The evolution equation for metric under the the normalised fiow is given by 



2H ■ hij + -hgij. 



dgij ^-^ I , 2~. 
c 

Proof. We compute 



dt ' n 



i^-'^{^t^p'^gij+^lJ'^^tg^j) 
^-2 (2i}^-hgij -tl;^2H ■ hij 
2 

-2H ■ hij + -hgij 
n 

~ ~ 2~ 
-2H ■ hij + -hgij. 

n 



□ 



Next we want to show that the mean curvature of the evolving normalised submanifold is 
bounded below by a constant greater than zero, and bounded above by a finite constant. As we 
know of no suitable isoperimetric inequality in high codimesion, we adapt Hamilton's intrinsic 
arguments in IIHa2ll to our setting. We will need to use the following fundamental results in 
comparison geometry to prove these estimates: 

Theorem 4.29 (Bishop-Gromov, Giinther volume comparison theorem). Let M be a complete 
Riemannian manifold and B{r) a ball of radius r in M. Denote by V^{r) the volume of a ball of 
radius r in the complete Riemannian manifold of constant curvature k. 

1) If the Ricci curvature of M is bounded below by Re > {n — l)kg, then 

vol{B{r)) < V''{r). 

2) If the sectional curvature of M is bounded above by some constant K > 0, then 

vol{B{r)) > V^{r). 

For a proof of these theorems we refer the reader to [IGHL1 . 

Lemma 4.30 (Klingenberg's Lemma). Suppose that M is a compact manifold and denote the 
length of the shortest closed closed geodesic in M by Ishon- If the sectional curvature of M is 
bounded above by some constant K > then the injectivity radius of M is bounded below by 

injg{M) > min I^L, ^Z,/,orf}- 



For a proof of Klingenberg's Lemma we refer the reader to [0. Since the second fundamental 
form of the evolving normalised submanifolds is bounded above, the length of the smallest closed 
geodesic must be bounded below, and therefore so too the injectivity radius: In order for a small 
loop to be forming, the second fundamental form must be blowing-up, and this is not the case. 
Heintze and Karcher derive an explicit lower bound for the length of the shortest closed geodesic 
in BHKH . although it suffices for our purposes to note that this is greater than zero. 
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Proposition 4.31. We have 
for all time t G [0, T). 

Proof. From equation (14.3 Icl ) the intrinsic sectional curvature of satisifies Kmin > 0. 
The Bishop-Gromov volume comparison theorem now implies vol(I]) < Cd", where d the diam- 
eter. From Bonnet's Theorem we also have d < C/^f\F\~, and thus V < C\H\J'J'^. In the 
normalised setting vol(S) = |So|, so |^|min < C and then (I4.31bl ) implies that |max < C □ 
Proposition 4.32. There exists as constant Cmm depending only on Sq such that 

\H\min ^ > 

holds for all time i £ [0, T). 

Proof. We work with the universal cover S of the normalised submanifold S. By the 
Giinther volume comparison theorem, the volume of S is some multiple its injectivity radius: 
vol(I]) > Cinj(S)". From the Gauss equation and the Pinching Lemma, the intrinsic sectional 
curvature of S is bounded above by some multiple of |i?|max> which is uniformly bounded above 
by the previous proposition. Moreover, since the second fundamental form of the normalised sub- 
manifolds is also bounded above, from Klingemberg's Lemma we obtain a lower bound for the 
injectivity radius. We may now estimate 

(4.32) vol(E) > Cinj(l;)" > > C\H\~l. 

The evolving submanifold is not undergoing any topological change before the singularity time, 
so by Bonnet's Theorem the first fundamental group of S is finite and constant in time. We have 

vol(S) = |7ri(S)| vols, 
and since both the first fundamental group and volume of the normalised submanifold are constant 
in time, vol(S) is also constant. This combined with equation (14.321 ) gives a lower bound on 
|^|max> and then equation (I4.31bl ) gives the desired lower bound on |min. D 

Proposition 4.33. We have 

/ \H\lax{t) dt = OO. 

Jo 

Proof. Follow the proof of Theorem 15.3 in IIHa21 with i?max replaced by |^ax ^^'^ 

^^\H\^ < A\H\^ + 2c\H\ljH\^ 

□ 

Proposition 4.34. The normalised flow exists for all time, that is, T = oo. 
Proof. We have di/dt = ip"^ and \H\'^ = V'^^I-H^P, so 

/ fi{i) di= [ h{t) dt = oo; 
Jo Jo 

however h < i/^ax ^ C'^ax ^^'^ therefore T = oo. □ 

The key step in the convergence argument is to show that the length of the traceless second 
fundamental form decays exponentially in time. Similar to the un-normalised setting, one consid- 
ers the scale-invariant quantity The reaction terms of the evolution equation for 
are again not quite favourable enough to use the maximum principle directly, and one proceeds in 
a similar manner to the un-normalised setting via intergral estimates (see [Hul] and below). The 
Stampacchia iteration is not needed, but only the Poincare inequality obtained from integrating 
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Simons' identity (see Propostion l4.36l below). We shall present a new argument based on the max- 
imum principle, which simplifies the existing argument by avoiding the use of integral estimates. 
Since this argument in new even for the case of hypersurfaces, we first treat the codimension one 
case as considered by Huisken in LHulj . 

Proposition 4.35. Suppose is an initially strictly convex hypersurface smoothly immersed in 
W"^^ moving by the normalised mean curvature flow. For all time t G [0, oo) we have the estimate 



|V/i|2 + |/i|2 <Ce 



Proof. The idea is to consider / := e\Vh\'^ + N\h\'^ /\H\'^ , where e > will be chosen small 
and sufficiently large. For the moment we work in the un-normalised setting. The evolution 
equation for |V/ip is of the form 

d 

— |V/ip = A|V/i|2 - 2\V^h\^ + h*h*Vh*Vh, 



so we obtain the estimate 



^|V/i|2 < A|V/i|2 - 2|v2/i|2 + ci\Hf\Vh\^. 



The evolution equation for is given by 

" 2 , . I ^,|2 , n , lll2 



(see Lemma 5.2 of IIHull and set a = 0). The importance of including the gradient term |V/ip 
in / is the following: the antisymmetric part of |V/ip contains curvature terms which we can use 
to obtain exponential convergence. We split V^/i into symmetric and anti-symmetric components, 
and upon discarding the the symmetric part we obtain 

where the last line follows from Simons' identity. Some computation shows 

\Rikjphpl + Riklpf^jpl — 4 ^ ^ (^j — K^K'j), 

then using that K^m > we estimate 



(4.33) > n4iJ/i|2 := ei\h\' 



The next important step is to estimate the term \H • Vihf^i — ViH • /ifcjp from below in terms 
of |V/ip. It is a relatively simple matter to estimate this term from below in terms of \VH\'^, 
however we want to use this good negative term to control the bad reaction term ci|i/p|V/ip of 
the evolution equation for |V/ip, so we need an estimate in terms of V/i. To do this, as always 
let h denote the second fundamental form and B a totally symmetric three tensor (we have Vh 
in mind). Consider the space A := {h,B : = l,|-Bp = 1}, and we also assume strict 
convexity of h. The conditions on h and B imply this space is compact. Now consider the 
function G{B) = \hpp ■ Bijk — hijBkpp\'^. We claim G{B) > 6 for some 6 > 0. Since A is 
compact, by the extreme value theorem G assumes its minimum value at some element of A. 
We show by contradiction that G / which proves the claim. The anti-symmetric part of G is 
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\Bipp ■ hjk — Bjpp • /life p. We compute at a point where G obtains its minimum, and rotating 
coordinates so that ei = VH/\VH\ we have 



n 

\Bipp • hj}^ — Bjpp • /ijfel = ^ ^ ^ifc^ ) 

k=l 

so \Bipp\'^ = or = X]fc=i ^ife- The latter implies that = h\i, which contradicts the strict 
convexity of the hypersurface. Therefore, if G{B) = 0, then I-Bjppp = 0. From the definition of G 
it now follows that the full tensor = 0. This contradicts = 1 and the claim follows. The 
term /ipp-i3jjfe — /iy i?fepp is a quadratic form, so for arbitrary /i and i? we obtain > (5|/ip|i?p 

by scaling. Applying this to our situation, we have \B\'^ = |V/ip, then estimating > nK^j^^ 
we obtain 

(4.34) \H ■ V,hu - V^H • hkif > <5n4i„|V/i|2 := eajV/ip. 

Returning now to the evolution equation for /, converting to the normalised setting and using the 
estimates (14.331 ) and (14.341 ) we get 

< Af - ei\h\^ + ci\H\^\vM^ + -1-{V,\^^^ 

at \H\^ \Hy 

\H\^ n 

In the normalised setting the second fundamental form, and therefore all higher derivatives, are 
bounded above. We can therefore estimate 

(V.I-H-p, Vi(e|V/ip)> < A\H\\VH\\Vh\\V^h\ < C|V/ip. 

Using < Cmin < I I min < |-^|max < C'max, wc make N Sufficiently large to consume the bad 
|V/ip terms and then we discard these terms. Using again Cmm < l-ff |min we estimate 

-ei\h\^ - —h\Vh\^ < -6f 
n 



for some small 5. We ultimately obtain 
This implies 



^f<Af + U''Vj-5f. 
ot 



ot 

and from the maximum principle we conclude e^^f<C and the theorem follows since |-H^|max < 

Cmax • I— I 

Note that we obtain exponential decay of both | V/ip and at the same time. Since we have 

pointwise control on the decay of exponential decay of the higher derivatives can be proved 
by the maximum principle in a similar manner as the un-normalised estimates. The important 

modification needed is that one adds in which is exponentially decaying, rather than 
which is not, to generate the favourable gradient terms. In our high codimension setting it is 
a simpler matter to estimate the good gradient term in the corresponding evolution equation for 

o 

in terms of |V/ip, and the same proof goes through provided we can estimate a lower 

bound for 

|V2/i|2 > ^\x7iVjhki-VkVihij\^ 



-:\RikaphjiaVp + Rikjphpl + Riklphj 



4 



'JP\ 



in terms of Such an estimate could hold for c < — 1), although at this stage we 

can no longer muster the patience to attempt the index gymnastics involved. In the absence of 
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this calculation, we give a sketch of the original proof contained in BHulll . with the necessary 
adjustments made for the high codimension. We remark that it would be nice to use the same idea 
in the un-normalised setting, and avoid the integral estimates. Unfortunately, at the moment we 
can only make such an argument work if the submanifold is already extremely pinched. 

Proposition 4.36. There exist positive constants C and 5 both depending only on Sq such that the 
estimate 

holds for all time t £ [to, oo), where to is some sufficiently long time. 
Proof. Consider the function 

|^|2 n' 

which is scale-invariant. The evolution equation for / is easily obtained from equation (14. 141) by 
taking a = 0: 

^f = Af- S{V,H,Vj) - ^\VH\^. 
dt \H\ \H\^ 



In the same manner as Proposition 15 . 5 1 we obtain the differential inequality 

where 6 is some small positive constant and the second integral on the right arises from differenti- 
ating the normalised measure. Using estimate (14.3 Ibl ) and l^j^ > Cmm we see there exists some 
time to such that for all t > to we have 

d_ 
dt 

for some smaller 5. This implies 



(4.35) ^ / jPdflg < -6 I hm^djlg 



e 



min 



from which the proposition follows easily. □ 
Proposition 4.37. We have the estimate 

for all time t > ti, where ti is some sufficiently long time. 

Proof. Consider / := \VH\'^ + N\H\'^\h\'^. This function is of degree -4 and from the 
relevant un-normalised evolution equations and Lemma l4.27l we derive 

The second term on the right can be absorbed by choosing N sufficiently large. The term 2 ( Vj | p , Vj | /i | ^ 

can be estimated by C{N)\h\ |V/ip, which after some time ti will be absorbed by the negative term 
V/ip. We use |i/|min > C'min to estimate the last term on the right, obtaining the differential 
inequality 

p =Kf + Ce-'' -5f 
which holds for all time t > ti. We then have 
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and from the maximum principle conclude 

/ < C{1 + i)e-^^ 
< Ce"^* 

for some 6 smaller again. The proposition now follows by integrating this estimate along geodesies 
and using that the diameter is bounded above. □ 

Proposition 4.38. For every m > 1 and p > 2, the estimate 

holds for all time t G [t^, oo). 

For a proof we refer the reader to Lemma 10.4 in HHull . The usual Sobolev inequality on a 
compact manifold now implies that max^^J V™/i| < Cm- With these higher derivative estimates 

o 

in place we can prove the crucial pointwise bound on \h\: 

Lemma 4.39. There exist positive constants C and 6 both depending only on Eq such that the 
estimate 

holds for all time t G [t4, oo). 

For a proof we refer the reader to Theorem 10.5 in [Hull . See also IIHa2i for the above two 
results. In particular, the reason why the Sobolev constant is uniformly bounded, and thus why 
can in fact use the Sobolev inequality is explained in IIHa2l . 

Proposition 4.40. The normalised submanifold converges uniformly to a smooth limit subman- 
ifold Soo as t ^ CO. 

Proof. The first step is to show Sqo is continuous. As we have done in the un-normalised 
setting, using Lemma 14.2 of liHa2J it suffices to show 

dt < C <oo. 



poo 


dg 


Jo 


di 



We estimate 

r-OO Q~ 







dt 



dt 



H -h 



n 



< C 

< C 
< 



Ce 



n 

hij Hgij 
n 



dt 

m 
1 

+ - 
1 

+ - 



(l^^ma) 



K)gv 



Hi) 



di 



I -^min I ) 9ij 



dt 



which is finite as desired. In going to the last line we have used Proposition |437]and Lemma |4.39l 
The proof that Scxd is smooth mimics that of the un-normalised setting, where here the exponential 
decay of the normalised estimates guarantees that the indefinite integrals in time which arise are 
finite. □ 

Proposition 4.41. The limit submanifold Sqq is a n-sphere lying in some {n + l)-dimensional 
subspace of^""^^ 



Proof. Lemma |4 . 39 1 implies that Sqo is totally umbilic. By the Codazzi Theorem (see |Sp 
n. 26] for a proof), the only closed, totally umbilic n-dimensional submanifold immersed ir 
^'^ is a n-sphere lying in some (n + 1) -dimensional subspace of M"+'^'. II 

The last proposition completes the proof of the second part of the Main Theorem |6] 



CHAPTER 5 



Submanifolds of the sphere 

In the previous chapter we studied the evolution of submanifolds of Euclidean space by the 
mean curvature flow. We now want to consider the situation where the background space is a 
sphere of constant curvature K. Our main result is the following: 

Main Theorem 7. Let = Fq{TP) be a closed submanifold smoothly immersed in S"+^. Tjf Eq 
satisfies 

|/i|2 < ^|^|2 + 2(^^^ n = 2,3 

< ;^|i?P + 2^, n>4, 
then either 

1) MCF has a unique, smooth solution on a finite, maximal time interval < t < T < oo 
and the submanifold Tit contracts to a point as t ^ T; or 

2) MCF has a unique, smooth solution for all time < t < oo and the submanifold 
converges to a totally geodesic submanifold Sqo- 

We highlight again that no assumption on the size H is required. The pinching condition 
< l/(n — + 2K implies that the submanifold has positive intrinsic curvature. A 

natural question to ask is whether some other geometric flow will deform all submanifolds of 
positive intrinsic curvature to either round points or totally geodesic submanifolds. In the case 
of hypersurfaces this problem has a very nice resolution due to Andrews. Beginning with the 
assumption that positive intrinsic curvature is preserved by some flow, in |EI the desired speed 
of the flow is found as an explicit solution of an ordinary differential equation. He then goes on 
to show that this flow does indeed deform an initial hypersurface of positive intrinsic curvature to 
either a point or a totally geodesic hypersurface. We point out that in the high codimension case 
such a theorem cannot be true (in dimension two) because of the Veronese surface. 

This proof of this theorem proceeds similarly to [Hu3l using the high codimension techniques 
developed in the previous chapter. After the relevant evolution equations are derived, we prove 
a version of the Pinching Lemma that holds in a sphere. The Pinching Lemma states that if the 
initial submanifold satisfies a certain curvature pinching, then the mean curvature flow preserves 
this pinching. A stronger pinching estimate is then deduced by a Stampacchia iteration argument. 
The essential content of this estimate is that in regions of large mean curvature, or after sufficiently 
long time, the submanifold is nearly totally umbilic. This estimate allows us to characterise the 
long time shape of the evolving submanifolds, which is completed in the last sections. 

5.1. The evolution equations in a sphere 

In the previous chapter we derived the evolution equation for the second fundamental form of 
submanifolds of arbitrary codimension in an arbitrary background space: 

^ dt^ij — ^l^ij ~\~ l^ij ' l^pql^pq ~\~ hiq ' hgphpj -\- hjq • hgphpi 2hip • hjghpg 
-\- 2Ripjghpg Rkjkphpi Rkikphpj -\- hijoiRkak/S^fi 
— 2hjpaRipoipl^l3 — 2hipaRjpal3'^(3 + ^kRkijji^P — ^ iRjkkpT^ji- 

In the case where the background space is a sphere the above evolution equation can be sim- 
plified significantly. If Ca, 1 < a < n + fc, is an arbitrary local frame for the background sphere, 
then in such a frame the Riemann curvature tensor takes the form 

(5.1) ^(ea,eb,ec,ed) = K((ea, ec)(eb, e^) - (ca, Cd) (e;,, Cc)). 
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The derivation of the evolution equation for follows that of the Euclidean case, however extra 
terms are now present due to the background curvature. We will show how to deal with these extra 
terms. First of all, as a sphere is a symmetric space, the first derivatives of the the background 
curvature are zero. The extra ambient curvature terms that remain are 

4:hijQ.hpqQiRipjq ^hijoi^hipoiRi^ji^p + '^hijahijjjRi^Q^i^p ^hijjjhipf^Rjpo^p 

Now, using the form of the Riemann curvature tensor of the sphere given by equation (15. It . for 
example, Rkak/s = K{6kk^ai3 — ^kp^ak)^ one finds various terms are zero or cancel, ultimately 
leaving only 

AK\H\^ - 2nK\h\^. 
The evolution equation for is therefore given by 

d 

(5.2) — = A|/i|2 - 2|V/ip + 2Ri + AK\H\'^ - 2nK\h\'^, 
or equivalently 

(5.3) — = A|/iP - 2|V/iP + 2Ri + 2K\H\'^ - 2nK\h\^ , 
at 

The ambient curvature terms appearing in the derivation of the evolution equation for can be 
dealt with similarly, and we obtain 

d 

(5.4) ttIHI'^ = A|i?P - 2|VFP + 2Ri + 2nK|ifP. 
ot 

The contracted form of Simons' identity takes the form 

(5.5) ^A|^|2 = hij ■ ViVjH + |V^|2 + Z + nK\h\\ 
where again 

Z = — ^ hijahijp^ — + Hahipahijphpjp. 

And finally, the basic gradient estimate 

(5.6) |V/i|2 > -^\VH\^ 

n + 2 

carries over unchanged. 

5.2. Curvature pinching is preserved 

We now prove the version of the Pinching Lemma that holds in sphere. Whenever we make 
reference to the Pinching Lemma in this chapter we obviously mean the following lemma. 

Lemma 5.1. If a solution F : T, x [0, T) — ?• S"+'^ of the mean curvature flow satisfies 
(5.7) 

at t = 0, then this remains true as long as the solution exists. 



\h\^ <i^\H\^ + ^K, n = 2,3 
\h\'^ <^\H\^ + 2K, n>4: 



Proof. The proof closely follows the Euclidean case. Here we consider Q = — aH^ — 
I3K, where a and /? are constants. Because we are allowing the initial submanifold to have H = 0, 
in order to compute in a local frame for the normal bundle where ui = H/\H\ we need to consider 
two cases: I) H = and 2) H ^ 0. For the first case, from the evolution equations for and 
we derive 

(5.8) —Q = AQ-2\Vh\'^ + 2Ri-2nK\h\'^. 
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In this case there is no need to split up directions of the second fundamental form as we did in the 
Euclidean case, and using the estimate of BLLIU on all the normal directions of Rl we get 

a,l3 i,j a, 13 



The reaction terms of (I5.81 i may therefore be estimated by 

2Ri - 2nK\h\^ < 3\h\'^ - 2nK\h\^ . 
If Q doesn't stay (strictly) negative, then = and 

31^1'' - 2nK\h\^ < -/3(2n - 3f3)K^ 

which is (strictly) negative as long as /3 < (2/3)n. This is a contradiction and the lemma follows 
in this case. Now consider the case H ^ 0. We may now work in the special local frames of the 
previous chapter, and the evolution equation becomes 

= AQ - 2(|V/iP - alVi/p) 
(5.9) dt ' 

+ 2Ri - 2aR2 - 2nK\h\'^ - 2n{a - l/n)K\H\'^. 
Arguing as in Euclidean case, if Q doesn't remain (strictly) negative, we may replace \H\'^ with 

o _ 

{\h\ — f3K)/{a — 1/n), and estimating as before we get 
2Ri - 2aR2 - 2nK\h\'^ - 2n{a - l/n)K\H\^ 

< 2\h^\^ - 2{a - -)\hi\^\H\^ + -\hiW " -(« " -)l^l' + 8|/liP|Lp + 3\L\^ 
n n n n 

- 2nK{\hi\^ + |L|2) - 2n(a - l/n)K\H\^ 
V n(a--)/ V n(a--)/ 



+ 2/3-4n+ ^ " . |/ii|"ir + 4 " ^ -n]\h.\'K 
V n(a--)/ \n{a--) > 

Vn(a--) / 

6 - (l/iiplLl^ + |L|^) + (2/3 - 4n + 

n{a--)) V n{a--)) 



-3\h-\'' + A[—^-—--n\\h_YK -2p[—^-—--n]K\ 
\n[a--) ) \n[a--) > 

We have rewritten the last line as such to highlight that after choosing the coefficient of the 
term as large as we can (namely 4/(3n)), we still have the good term — 3|/i_p left 
over. The last line above is a quadratic form, so by requiring that its discriminant be negative we 
will have a strictly negative term. The discriminant is 

n 



which is negative for our values of a and /3 in dimensions two to four. For dimensions n > 4 the 
best value of a we can expect is — 1), and so with this restriction, the amount of the good 
terms is increases to —2{n — 4) — 3. The discriminant is now 

^ = «(^£t)-"){2(^£t))-(2(»-4)+3)^}. 

^ n' ^ n' 

and which is strictly negative for /3 = 2 for all n > 4. The most restrictive condition on the size of 
/3 comes from the coefficient of the term, which gives the values of (3 in the statement of 
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the lemma. With the chosen values of a and /3 the right hand side of equation is strictly negative, 
which is contradiction, and so Q must stay strictly negative. □ 

We now want to formulate a slightly different statement of the Pinching Lemma that will be 
useful in later setions. For e > 0, set 

ra> ^ = 2,3 f/3,:=f(l-6), n = 2,3 

;i=T+7, ri>4 := 2(1 - e), n > 4. 

If the strict inequality < a|i/p + fiK holds everywhere on the initial submanifold, then there 
exists an e > such that < Qe|//p + fi^K on Sq. On the other hand, if equality of the 
pinching condition holds somewhere on the initial submanifold, that is < a|//p + fiK, and 
the pinching does not immediately improve, then the same strong maximum principle argument 
as in Proposition 14.41 of the previous chapter shows that Sq = x S"~p, where < p < n. 
If p = 0, then So is a totally umbilic sphere, in which case there exists an e > such that 
< Oel-ffp + ^^K holds. If p 7^ 0, the above-mentioned product of spheres all he outside 
of the pinching cone being considered. Therefore, if the equality of the pinching condition holds 
initially, after some short time the submanifold satisfies < a^lffp + ji^K for some e > 0. 

5.3. Pinching improves along the flow 

In this section we prove an important estimate that allows us to characterise the asymptotic 
behaviour of the solution. As mentioned in the introduction to this chapter, the essential content 
of this theorem is that in regions where the mean curvature is large, or after long enough time, the 
submanifold is increasingly becoming totally umbilic. This can be interpreted by saying that the 
pinching improves along the flow. 

Theorem 5.2. There exist constants Co < oo, o"o > 0, and (5o > all depending only on Eo such 
that for all time < t < T < oo, the estimate 

\h\^ < Co{\H\^ + K)^~''o^-^ot 

holds. 

For technical reasons it is more convenient to work initially with the auxiliary function := 

\h\^/{a\H\^ + peKf'", where a := l/{n{n - 1 + e)). 

Proof. We begin by deriving the evolution equation for fa-. 
Proposition 5.3. For any a £ [0, e/2] we have the evolution equation 

(5.10) a\H\+l3,K^ " {a\H\^ + I3,KY'''^ ' -"^ 

+ 2ct|/i|2/<x- 

Proof. From the evolutions equations for and \H\'^ we get 
A|/i|2 - 2|V/i|2 + 2Ri + AK\H\'^ - 2nK\h\'^ 



dtfa 

(5.11) 



(a|i/|2 + /3,i^)i-'^ 
1 {A\H\^ - 2|Vi?p + 2R2 + 2nK\H\'^) 



The Laplacian of is given by 

_ a{l-a){\h\'-l/n\H\-) , a'{2 - a){l - a){\h\' - l/n\H\') , , 
(a|if|2 + /3,i^)2— I I ^ (a|i/|2 + /3,i^)3— MM- 
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Using this expression for the Laplacian as well as the identity 



alHP + l3,K^ ' '' ' ' (a|Jf|2+AK)2 
equation 15.1 It can be manipulated into the form 
(5.12) 



+ 



2aaR2fa 
a\H\'^ + P,K' 



The gradient terms on the third line are non-positive under our pinching assumption and we discard 
them. Using equation (15.61 ) and the Pinching Lemma we estimate the useful gradient terms on the 
second line as follows: 

|V/i|2 - i|VF|2 - , ,> |Vg|^ 



a\H\'^ + p,Ky ' n' ' a\H\'^ + P,K' 

< -2 [ 3 _ 1 _ a((«,-l/n)|gp + /3,^) ) 

< , /^-l-4ivg|^ 

~ a|i7|2+/3,K\n + 2 n J ' ' 
_^ -2ev|Vgp 
■~ a|//|2 + /3,i^' 

where importantly, ey is positive for all n>2. Next we estimate the reactions terms on the second 
last line of (15.121 ). Expanding these reaction terms in the special local frames and estimating 

W>{\h\'^- PeK)/{a, - l/n) we obtain 



(a|i7|2 + /3,K)i-^ [ ' n ' ' a\H\^ + P^K a\H\'^ + jS^K 

2 { a / ^ 1 ^ ' ? '4i ? i2 



Of — - V 2 niaf — -]^ Of — -V2 n — - k 

3 a ^ 2a ^3 an[2 — (cr + £2)] ^ ^ ^^|° ,4 



2a,-- n(a, --)2 2 /3,(a, _ _) 
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-Pe[ , ' 1, +n(l-ei) ^ 1 ]\h^YK 



n 



- ^ { + n(l - e,) - -J^^JliZp^) 

\n{a,--y a,-- / 

- nl3,eiK^\h\'^ - ane2K\h\'^\Hf 

Provided ei, 62 and a are all chosen sufficiently small all terms in the above expression can be 
made negative, and we discard them with the exception of the two terms on the last line. The 
above terms contain two quadratic forms, which are estimated in a similar manner to the Pinching 
Lemma. We have 



{a\H\^ + fSeKy-"" { n ' ' a\H\'^ + p^K a\H\'^ + 

< -2nmin{ei,e2}i^/CT 
:= -2ne'Kf,. 

Finally, we estimate the last term on the right of equation (15.101) by i?2 < and the 

proposition is complete. 

□ 

As in the prevous chapter, the small positive 2a\h\'^ prevents us from using the maximum 
principle and we proceed by deriving integral estimates and an iteration procedure. The thrust of 
this iteration procedure is to exploit the good negative jV/fp term in (15.101 ) using the contracted 
Simons' identity and the Divergence theorem. In order to do this we need a lower bound on the 
Laplacian of f^j, and as Huisken points out in fHu3], this can be achieved because the pinching 
condition (compare the pinching condition in [Hu3 |) implies that the submanifold has positive in- 
trinsic curvature. The next estimate is the part of the argument that relies on the intrinsic curvature 
of the submanifold being positive. 

Lemma 5.4. Let Sq be a n-dimensional submanifold immersed in a spherical background of 
constant curvature K. If^Q satisifies < a|-ffp + 13 K, where 

J^A. " = 2,3 _|^, „ = 2,3 

then there exists a positive constant e depending only on Sq such that the estimate 

Z + nif > e|^p(a|i?p + f3,K) 

holds for all time. 

The proof of this lemma is similar to Lemma |4. 11 l of the previous chapter. For the same rea- 
sons as in the Pinching Lemma, the cases H = and H need to be examined seperately, and 
again the case H = Ois treated easily. Let us briefly examine the case H ^ 0. The computations 
are the same as those in Lemma l4.11l and one finds, in dimensions two to five 

1,°,0, ,0 1,° ,0, ,0 n 2 ° .A ,°,0,° ,0 

Z + n h^K > - hi ^ + - hi f H f + - /i_ pLff 2 /ii ^ -4 /ii 2 P 

n n 2 

- - h\^ r hif H ^ 7 +n( hi ^ + h\^)K, 

2' ' 2n(n-l)' ' ' ' 2n(n-l)' 1 1 1 ^ vi n ti \ j , 
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and in dimensions six and higher 

The size of a is computed in the same way as the Euclidean case, and for /3, in all dimensions 
n > 2 we require 



/3 /I n-2 



Q, _ i Vn 2n(n — 1) 



+ n > 0. 



In dimension two and three this gives /3 < 2(n — l)/3, which is more restrictive than that required 
by the Pinching Lemma. For n > 4 we require /? < 2. We now commence with the integral 
estimates. 

Proposition 5.5. For any r] > Owe have the estimate 
[ fP{a\Hf +^,K)dVg 

(2pr/ + 5) f fr' \^rn\iV i " ^) / fP'^W f \' HV 



< 



Proof. The proof this lemma follows Lemma l4.11l of the previous chapter. Using the con- 
tracted Simons' identity and A I i7 1 2 = 2\H\A\H\ + 2|V|i7| the Laplacian of can be written 
as 

2 ° 2 2 ° 



2^ 



a|iJ|2 + (3,K- 



We want to estimate Af^ from below. The first term on the third line is non-negative and we 
discard it. Working with the last two terms of line three, using the Kato-type inequality |V|ff| p < 
I V/fp, equation (15.61 ) and the Pinching Lemma we estimate 

2 i 3 I a\h\' 1| ,2 

- {a\H\^ + (3,Ky~^\n + 2 n a|i?|2 + J ' 

which is non-negative and we discard this term. Having discarded these terms we are left with 

2 2 2 
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We now multiply this equation by fa~^ and integrate it over the submanifold. The terms integrate 
as follows: 

/ fp-'Au = -{p-i) [ /r^iv/.p dv,- 

In performing the integration we have made use of Green's First Identity, the Codazzi equation 
and the Divergence Theorem. We discard two terms that have the appropriate sign, noting that 
the terms with and inner product do not have a sign. After some factoring and rearranging, we 
estimate the terms with an inner product using the Cauchy-Schwarz and Kato inequahties to obtain 

^^^X. HHpfft.).-^ '^"'''''-' 

o 

Using the Peter-Paul inequality, as well as the inequalities < /^(ali/p + P^KY~'^, < 
{a\H\'^ + jie^Y and (1 — a) < 1, we estimate each term on the right as follows: 

^ ^ X. + c - X. (.|g/;lAT-- '^^i'^'^'^ 
^^'"-^'"X. (»i^/f'ftg)'-' '"^- 
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We use Lemma l4.1 ll to estimate the two terms on the left: 

Putting everything together with a little rough estimation of the coefficients to coax them into a 
more convenient form we obtain 

2e / fP{a\H\'' + PeK)dVg 

<{3prj + 10)[ fP(a\H\^ + f3,K)dVg = ^^^^ [ /r'| V/.p dF^. 

JSi V JT,t 

Dividing through by 2e completes the proposition. □ 

The next step is to show that sufficiently high norms of fa are bounded, and in fact decay 
exponentially in time. 

Proposition 5.6. For any p > 8/(ev + 1) we have the estimate 

-2neKp f f;dV, + 2ap [ /J(o|ff p + ftK) rfV,. 

Proof. We differentiate under the integral sign and substitute in the evolution equations for 
/cr and the measure dVg to get 



d 
dt 



s 



-id fa 



^-^^dV, 

(5.13) 

rp-l 



< _ 1) ^ /r'lV/.l' dVg + 4(1 - a)p l^^^j^^^\H\\V\H\\\Vfa\ dVg 



We estimate the second integral by 
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then substituting this estimate back into (15.131 ) gives 

2p 



dt 



<[-pip-l) + j) /r'l V/.P dVg 
+ 2neKp [ fPdVg + 2ap [ \h\^fPdVg 

-2neKp [ fPdVg + 2ap [ \hff^dVg. 

JT. JT. 

We want to choose /x so that 1 — 2/{^{p — 1)) > 1/2 and p so that 1 — /x/ey > 1/2. We 
therefore choose /x = 4/(p — 1) and p > max{2,8/(ev + 1)}. In the last term we estimate 
< + f^eK), and the proposition is complete. □ 

Lemma 5.7. There exist constants C2 and C3 both depending only on Sq such that ifp > C2 and 
c < Cs/^yp, then for all time t G [0, 00) we have the estimate 

Proof. Combining Propositions 15.51 and 15 . 6 1 we get 



_ . M2p^x / ^ , \VH\UVg 

-2neKp [ f^dYg. 

JT 

Recall we are already assuming that p > max{2, 8c/ (ey + 1)}. Now suppose that 

cr < - 

Set 7] = Aca/e, then 

{4£ _ 1 

We require p > 25/ev for the last inequality to hold. With these assumptions on p and a we have 

I ^ dVg < -2neKp £ dVg, 



^-2neKpt 



and thus 

d_ 

dt Jy." " " Jt " "4=0 
This implies the lemma with Ci = (|So| + 1) ™aXo.g[o^i/2] Sq/o-), 61 < 2neKp, 
C2 := max{8/(ev + 1), 25/ev} and C3 := ' □ 



fldVg<- f fldVg 
Jt 



5.4. A GRADIENT ESTIMATE FOR THE MEAN CURVATURE 81 

Lemma 1577] shows that for a sufficiently small, sufficiently high norms of are bounded 
and exponentially decaying in time. We can now proceed as in IIHu3ll via a Stampacchia iteration 
procedure to uniformly bound g^r := /cre^'^^/^)*, from which the theorem easily follows. We point 
out that during the course of this argument, a is fixed sufficiently small once and for all. □ 

5.4. A gradient estimate for tlie mean curvature 

Here we establish a gradient estimate for the mean curvature. This estimate is required in the 
following section to compare the mean curvature at different points of the submanifold. 

Theorem 5.8. For each rj > there exists a constant depending only on rj such that for all 
time the estimate 

holds. 

We begin by deriving an evolution equation for \VH\'^. 

Proposition 5.9. We have the evolution equation 
(5.14) 

^1 Vi^p = A| V/?|2 - 2|v2Fp + 2^(i?(4, d^)H - Vp{R{dk, dp)H) + R{dk, dp)VpH 
+ RcpkVpH + H ■ hpqhpq + nKH, VkH) . 

Proof. We compute 

d d 

— \VH\'^ = — {VkH,VkH) 

= 2g{VkVtH + R{dk, dt)H, VkH) 

(5.15) = 2^{Vk{AH + H ■ hpqhpq + nKH) + R{dk,dt)H, VkH) . 

The proposition now follows after the use of following two identities: 

A|Vi?|^ = 2g{AVkH,VkH) + 2|V^i7|^ 
AVkH = VkAH + Vp{R{dk, dp)H) + R{dk, dp)VpH + RcpuVpH. 

□ 

Corollary 5.10. There exist constants A and B, depending only on Sq, such that we have the 
estimate 

d 

Proof. We estimate the reaction terms of (I5.14l i. Using the spacelike Gauss and Ricci equa- 
tions, all the reaction terms except the first one look like /i*^ * VH*"^ and VH*"^. We need to use 
the timelike Ricci equation to estimate the first reaction term, however simple estimation of the 
ambient curvature term in the timelike Ricci equation gives rise to a term that looks like H*^*VH. 
A closer inspection of this term shows that in fact it is zero: 

g{R{FA, F,dt)H, VkH) = g{R{F,dk, H)H, VkH) 

= {F,dk,H){H,VkH) - {F,dk,VkH){H,H) 
= 0. 

All the reaction terms now look like h*'^*VH*'^ and VH*'^, which we can estimate by yl|//p|V/ip 
and i?| V/ip using the Pinching Lemma the Cauchy-Schwarz inequality. □ 

We need two more estimate to complete the proof. 
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Proposition 5.11. We have the estimates 

(5.16) ^\H\^ > AliJp - 12|iJP|ViJp + -\H\^ 
at n 

^((iVi + N2\H\^)\h\'') < A((iVi + N2\Hf)\hf) - ^\ n2 - l)\H\^\Vhf 



(5.17) 

-(iVi - ci(iV2))|V/i|2 - C2{N,,N2)\h\'{\H\' + !)• 



3n 

Proof. The evolution equation for \H\'^ is easily derived from that of \H\'^: 
8 

— = AliJp - 2|V|i/P|2 - 4|ifp|Vii'P + iR2\H\'^ + inKlHl'^. 
ot 

We discard the last term and the proposition follows from the use of |V|ff|p < |Vffp and 
R2 > l/n\H\'^. To prove (15.171 ). from the evolution equations for and \H\'^ we derive 

^^{[N, + N2\H\%\') 

= A((iVi + N2\H\'^)\h\'^) - 2N2{Vi\H\'^,Vi\h\'^) - 2N2\h\'^\Vh\'^ + 2N2R2\h\'^ 

- 2{Ni + N2\H\'^)i\Vh\'^ - -|VF|2) + 2(iVi + N2\H\'^)iRi - -R2) - 2nKNi\h\'^. 

n n 

We estimate the second term on the right as follows: 

-2iV2(Vi|F|2,Vi|^|2) < 8Ar2|/i||^||Vi/||V/i| 

< SN2\H\./^\Vh\^Co{\H\^ + 

<^!^\H\'\Vh\'+c,{N2)\Vh\\ 

o 

Using Young's inequality, R2 < \h\'^\H\'^, and Ri — l/ni?2 < 2|/ip|/ip we estimate 

2N2R2\h\'^ + 2{Ni + N2\H\'^)iRi - -R2) < C2iNi, N2)\h\WH\^ + 1). 

n 

The constants depend on more that just A'^i and A''2, however we only highlight the dependence on 
N as this is relevant in the following proof. We discard the last term on the right, and equation 
(15.171) now follows. □ 

Proof of Theorem [5T81 Consider / := |ViJp + (A^i + N2\H\'^)\h\'^ . From the evolution 
equations derived above we see / satisfies 

y<M + A\H\^\Vh\'' + B|V/i|2 + ^^^^(iV2 - 1)1^1 

+ ^^^^(^1 - ci(iV2))|V/i|2 + C2{N^,N2)\h\\\H\'' + 1). 

We choose A^^i and N2 large enough to consume the positive terms arising from the evolution 
equation for |Vi/p. This leaves 

|/ < A/ + ^^^(iV2 - m\'m' - ^^^(ivi - C3(iv2))iv/.|2 

+ C4(iVl,iV2)|^|'(|^|^ + l). 



5.5. ASYMPTOTIC BEHAVIOUR OF THE SOLUTION 
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Now consider g := e^''"/^)*/ — From the above evolution equations we have 

|(e(^o/2)*/-r/|i?r) 

< |e(5o/2)i (|v^|2 + (iVi + N2\H\^)\h\^) 

+ e(^o/2)t (a; + ^^^{N, - l)\H\'\Vh\' - ^^^{Ni - csiN,))\Vh\' 



+ c^iNi,N2)\h\WH\^ + 1)) - 7](A\H\^ - 12\H\^\VH\^ + -\Hf). 

The terms on the first line can be absorbed into those on the second line by suitable estimation. 
By choosing sufficiently large the gradient term on the last line can be absorbed, and then we 
choose A^i larger again to make the |V/ip term negative. We finally discard the negative gradient 
terms to get 

^g<Ag + C5e('^«/2)*|/;|2|/7|4 _ ^\h\\ 
at n 

Using Theorem 15.21 then Young's inequality we obtain 

ot 

from which we conclude g < cj. The gradient estimate now follows from the definition of g. □ 

5.5. Asymptotic behaviour of the solution 

In this final section we study the long time behaviour of the solution. Two limit profiles 
are possible, determined by whether or not the mean curvature blows up. We first examine the 
case where the mean curvature becomes unbounded. We do this by using the gradient estimate 
and Bonnet's Theorem to compare the submanifold at different points. In the case of a spherical 
background, the Chen's estimate combined with our pinching condition gives 

(5.18) K^nix) > li-^ - a^\H\'{x) + ^l^iMK. 

Theorem 5.12. If\H\ 

ma.v — > oo as t — > T, then T must be finite and diam — )> t — )■ T. 

Proof. From Theorem 15.81 we know that for any r/ > there exists a constant such that 
\VH\ < r?|i7p + on < t < T. We highhght that at this stage, T could be infinite. Since 
by assumption |i/|max — )• oo as t — )• T, there exists a r(r/) such that C,^j2 < l/2ri\H\'^^^ for all 

T <t <T. Thus \VH\ < r]\H\l^^ for all t > r. Fix some a G (0, 1) and set r] = ^^i^. 
Let t G [T(r/), T), and x be a point with \H\ {x) = |//|max- Along any geodesic of length 
from X, we have \H\ > |i^|max — — ^l-f^lmax = (^\H\maJ^■, and consequently the sectional 
curvatures satisfy K > e'^a'^\H\'^^^. From Bonnet's Theorem it follows that diamS < — , 
from which we conclude that |min > cl^^lmax on the whole of for t G [r(r/), T). 

The previous line shows that by choosing r sufficiently large, |if|min can be made arbitrar- 
ily large. It follows from Theorem 15.21 that after some sufficiently large time the submanifold 
is as pinched as we like (and in particular can be made to satisfy < l/(n — l)|ifpin 
dimensions n > 4 and < 4/{3n)\H\'^ in dimensions 2 < n < 4). We now show that 
once the submanifolds are pinched as such, the maximal time of existence must be finite. De- 
fine Q = — a|/ip — b{t), where a = ^ and b is some time -dependent function. Because 
l-f^lmin > and the submanifolds are as pinched as we like, for some sufficiently large time r we 
can choose a 6(r) = 6,- > such that Q > for t = r. The evolution equation for Q is 

—Q = AQ- 2(\VH\'^ - alV/iP) + 2R2 - 2aRi + 2(n - a)K\h\^ + 2anK\H\^ - b'it) 
ot 

>AQ- 2(|Vi?p - a|V/ip) + 2R2 - 2aRi - b'{t). 
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Estimating the reaction ternis as before we obtain 

2R2 - 2aRi - b'{t) 

= Y.{Y1 ^"^ii") ^ - 2a ( hij^h,jp^ ^ - 2a\R\^ - h'{t) 

ij ot a, (3 ij 

000 2 ° 00 

> 2|^ip(a|^ip + a|^__p + b) + — --(a|^_p + b){a\hi\^ + a\h_\^ + b) 

n(l — a/n) 

- 2a\hi\^ - 8a|^i|^|Lp - 3a|L|^ - b'{t). 

Equating coefficients, we find Q > is preserved if ^ < We can therefore take 

nbo 



b{t) 



n - 8bo{t - t)' 



This is unbounded as t — )• r + so we must have T < r + g^. 



□ 



Let us now consider the case where the mean curvature stays bounded for all time. 
Theorem 5.13. If\H\ 

max remains bounded, then the flow exists for all time and converges to a 
totally geodesic submanifold Sqo- 

Proof. Since |-fr|max is bounded, from Theorem I5.2l and Theorem 15.8 1 we have the estimates 

(5.19) 1^1^ < Coe-''°* 
and 

(5.20) |Vi?p < Ce-('^«/2)*. 

From (I5.I8I 1 we know that smallest sectional curvature is positive, so from Bonnet's Theorem it 
follows that the diameter of is bounded. Using this fact, integrating the second estimate of 
(15.201 ) along geodesies gives 

(5.21) \H\^,, - l^/Uin < Ce-(^°/2)t_ 

Now observe that if the time of existence is infinite, then |-?/|min must remain zero: From (I5.19t 
it follows that after sufficiently a long time the submanifolds are again as pinched as we like, 
and if |i?|min > 0, then the same argument just given in the previous case would show that T 
must be finite. Therefore we must have |-?/|min = 0. From equation (I5.21l i it now follows that 
l^^lmax < Ce^^*^"/^)*. Thus \H\^ decays exponentially and consequently 

We now have all the necessary estimates in place to repeat the convergence arguments of the pre- 
vious chapter to obtain smooth exponential convergence of the submanifolds to a totally geodesic 
submanifold. □ 



CHAPTER 6 



A partial classification of type I singularities 

In Chapter |4] we show that if a submanifold satisfies a suitable pinching condition, then the 
mean curvature flow evolves the submanifold to round point in finite time. In this chapter we 
relax the pinching of the initial submanifold and seek to understand the asymptotic shape of the 
evolving submanifold as we approach the maximal time of existence. We still assume that \H\ 
is everywhere positive initially, however having relaxed the pinching assumption, we no longer 
necessarily expect the entire submanifold to disappear at the maximal time. In the case of mean- 
convex hypersurfaces, a classification of type I singularities was achieved by Huisken in IIHu4ll 
and fHu51. A key ingredient in this analysis was Huisken's monontoncity formula, introduced 
in |Hu4 |, which also holds in arbitrary codimension. The singularities classified by Husiken in 
IIHu4i are a special kind of type I singularity called a 'special' type I singularity. The more general 
kind of singularity is naturally called a 'general' type I singularity and in order to have a complete 
understanding of type I singularity formation, it is desirable to be able to treat general singularities 
(definitions of the various kinds of singularities follow). In the case of embedded hypersurfaces, 
the classification of general type I singularities is due to Stone [St]. 

Here we follow |Hu4| and |Hu5 1 to give a partial classification of special type I singularities of 
the mean curvature flow in high codimension. Instead of using the continuous rescaling argument 
used in IIHu4L we proceed slightly differently by considering a sequence of parabolically rescaled 
flows. Huisken's original argument is recast in terms of rescaled flows in [E| and also IIW5I . We 
cannot make Stone's argument to classify general type I singularities work in high codimension, 
essentially because a pointwise curvature condition does not seem enough to conclude that em- 
beddedness is preserved. We point out that even in the codimension one case, the classification of 
general type I singularities of immersed hypersurfaces is an outstanding problem. For an excellent 
account of singularity analysis in the mean curvature flow of hypersurfaces (as well as a wonder- 
ful introduction to the mean curvature flow) we recommend to the reader the recent book [M | by 
Mantegazza. 

A similar classification of type I singularities of the mean curvature flow in high codimension 
has previously been obtained by Smoczyk [Sm3|. In [Sm3| Smoczyk classifies blow-up limits of 
the the mean curvature flow that have flat normal normal bundle. Although this curvature condition 
is much more restrictive than the pinching condition we have been working with, Smoczyk's 
classification includes additional submanifolds that do not feature in our classification, namely, 
products of Euclidean space with an Abresch-Langer curve, which also appear in the hypersurface 
classification. These spaces do not appear in our classification as they do not satisfy the pinching 
assumption. It's worthwhile to point out that the condition of having flat normal bundle is not 
preserved by the mean curvature flow. 

With regards to this last chapter, we wish to express our gratitude to Patrick Breuning for 
sending us a draft of his PhD thesis [Br], in which he improves upon and extends Langer's com- 
pactness theorem [La] to submanifolds of arbitrary codimension. We would also like to thank 
Andrew Stone for friendly correspondence and for sending us a copy of his PhD thesis. 



6.1. The blow-up argument 

We shall need some basic concepts from measure theory in this chapter. We remind the reader 
that S is a fixed manifold, and that St := -Ft(S) refers to the immersed submanifold. For a 
function / G C°° x M) defined on the ambient space, we follow standard abuse of notation 
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and denote 



f{F{p))dflg{p)= / f{p)dflg{p). 



Integration over the manifold S with respect to dfig and integration over the image in M"+'^ 
are linked by the area formula. We denote the pushforward measure by ^ = F{fig), where for 
U G M"+'^ an open set, /x(C/) := fig{F^^{U)). In order for the pushforward measure /i to be a 
Radon measure, the immersion F must be a proper immersion. Recall an immersion is proper if 
the inverse image of a compact set is also compact. The area formula relates the induced measure 
on S to the Hausdorff measure on R^^*^ restricted to the image F{T,) of the immersion. We 
denote Hausdorff measure on the ambient space by d'X'"'^^ or simply by dJi. For a /ig -measurable 
function / : S — )■ M, by the area formula we have 



[ fip)d^g{p)= [ I fip)] dJi'^-'H^)- 

Choosing / = X[f-i(f(S))] gives 

[ df,gip)= [ IE) ^^""''(^)- 

Thus, denoting by 6 the multiplicity function, we have fi = "KlB. In particular, if F : — )• ]R"+'^ 
is a properly embedded submanifold, then 6 = 1 and 

dugip) = [ d5{"+^(a;). 

For more details on Hausdorff measure and the area formula we refer the reader to IEG1 . One 
final piece of notation before getting underway, for a point p G S, we put linit^+T F{p) '■= P ^ 

In order to study the asymptotic shape of the evolving submanifold St around a singular point 
as the first singular time is approached, we progressively 'magnify' the solution around this point 
by considering a sequence of rescaled flows. The limit of such rescaled flows is called a blow-up 
limit. Our first task is to show how to obtain such a limit. In order to obtain a smooth blow-up, we 
assume that the submanifold is developing a so-called type I singularity. This imposes a natural 
maximum rate at which the singularity can develop, which then enables us to rescale at a rate that 
keeps the maximum curvature of the rescaled solution bounded. 

A submanifold is said to be developing a type I singularity at T if there exists a constant 
Co > 1 such that 

2/. ^\ ^ C'o 



max|/i| [p, t) < 



pes ' I ^ - 2(r - t) ■ 
The blow-up rate of any singularity also satisfies the lower bound 

max|/ip(p, t) > 



(see HHuli or ||Ml)^ and so in the case of a type I singularity 

1 I T 1 2 / \ ^0 

< max|/i| [p, t) < 



2{T-t) - pes' ' ' - 2{T-t)' 

Let g G E be a fixed point and assume that the type I condition holds. We want to rescale 
the solution around the point q G M"+'^ by remaining time. Let (tfc)fcgN be any sequence of times 
such that tfe — T as A; — oo. For example, we could take tf^ := T — 1/k. To rescale by remaining 
time we set the scale = l/y^2(T — tk). We then define a sequence of parabolically rescaled 
flows 

(6.1) Fk{p,s) ■.= \k{F{p,T + s/Xl)-q). 



6.1. THE BLOW-UP ARGUMENT 
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Then for each k, Fk : x [— A^.T, 0) is a solution to the mean curvature flow (in the time vari- 
able s) that exists on the time interval s G [— A^T, 0). Under our parabolic rescaling the second 
fundamental form rescales like = |/ip/A|, so using the type I hypothesis 

-^^^"^^•^■2(r-r-s/A2) 

^ -Co 
2s 

which holds on s £ [— A^T, 0). Consequently, on the time intervals := (— A|T, the 
rescaled flows have bounded second fundamental form. We would now like to apply a com- 
pactness theorem for immersed submanifolds in order to obtain a limit flow. The compactness 
theorem usually quoted in this context is llLal . The result presented in llLal is for a sequence of 
two-surfaces of Euclidean three-space with -bounded second fundamental form and a global 
area bound, whereas we need to apply the result to a sequence of n-dimensional submanifolds of 
codimension k in the presence of bounds on all higher derivatives of the second fundamental form 
and only a local area bound. Very recently we learnt that in his PhD thesis Patrick Breuning has 
extended Langer's result to submanifolds of arbitrary codimension in the presence of a local area 
bound I Br I . We record Breuning 's compactness theorem as follows: 

Theorem 6.1 (Breuning-Langer compactness theorem for immersed submanifolds). Let : 

— )• be a sequence of proper immersions, where is a n-manifold without boundary 
and G Fk{Mk). Assume the following conditions are satisfied: 

1) Uniform curvature derivative bounds: 

For each /c G N, for every m G N there exists a constant Cm (R) depending on m and R 
such that I V^/ifcliT'j. < Cm- 

2) Local area bound: 

For every R > there exists a constant Cr depending on R such that ^^{Br) < Cr. 

Then there exists a proper immersion F^ : M^o — )• W^^^, where M^o is again a n-manifold 
without boundary, such that after passing to a subsequence there exists a sequence of diffeomor- 
phisms (j)]. : Uj. {FjS)~'^[Bij) C M^, where C M^o are open sets with CC C/fe+i and 
Moo = Uj^i ^^'^^ f^'^t (t^*kPk\uj converges in C°°{Uj,M.'^) to Foo\uy 

This is the essentially the statement of the Breuning's theorem in his thesis; we have simply 
changed some notation to conform with our own. Note Breuning states the local area bound in 
terms of the pushforward measure. Before we learnt of Breuning's compactness theorem we did 
not know whether Langer's theorem did in fact hold in arbitrary codimension and we produced 
the following compactness theorem for immersed submanifolds in arbitrary codimension using 
the well-known compactness theorem of Cheeger and Gromov for abstract manifolds. We refer 
the reader to IIHAII for an introduction to Cheeger-Gromov convergence, and to HHaTl lPl for proofs 
of the Cheeger-Gromov compactness theorem. We have been influenced by the treatment in [HA|. 

We consider the following notion of convergence of a sequence of immersed submanifolds: 
For each G N, let be a complete smooth manifold, F^ : — )■ a smooth immersion 
and pk G Mfc a basepoint. We say that (Mk,Fk) converges to (Mqc-^oo) on compact sets of 

X M if there exists an exhaustion {^7fc}feGN of -^oo and a sequence of smooth diffeomorphisms 
(t)k:Uk^Vk C Mk satisfying: 

1) For every compact K C Moo, (j^l^klK converges in C°°{K,M^) to Foq\k 

2) For any compact ^ C there is some ko £ N such that {4>lFk){Uk)nA = Fk{Mk)nA 
for all k > ko 

We remark that the Langer-Breuning compactness theorem in the form we have stated it is not quite 
satisfactory for our purposes, as it does not address the second criterion of the above definition of 
convergence. 



88 



6. A PARTIAL CLASSIFICATION OF TYPE I SINGULARITIES 



Theorem 6.2. Suppose {F^, Mk)k^f^ is a sequence of proper immersions : — >■ M of 
smooth complete n-dimensional manifolds that satisfy the following conditions: 

1) Uniform curvature derivative bounds: 

For each k G N,/or every m € N there exists a constant Cm{R) depending on m and R 
(and independent ofk) such that |V^/ife|i?j. < Cm 

2) The sequence (-Ffc)fceN does not disappear at infinity: 

There exists a radius R> such that Bji{0) Ci F^^M^) =^ %for all k E N. 

3) Local area bound: 

For every R > there exists a constant Cr depending on R (and independent ofk), such 
that 



Then there exists a subsequence of {F^, Mk)keN which converges on compact sets o/R^ x R to a 
complete proper immersion (Mqo, F^) that also satisfies the the same local area bound. 

Before we prove this theorem, we mention two issues that need to be dealt with in the proof. 
First of all, the limit produced by applying the Cheeger-Gromov compactness theorem is an ab- 
stract limit that a priori loses all knowledge of the background space. Second, the limit may be 
disconnected (e.g. a lengthening cylinder), and so the metric produced by the Cheeger-Gromov 
compactness theorem only sees the connected component of which is is part. We therefore need 
to take care to capture all coimected components of the Umit. 

Proof. Let gk denote the metric induced by F). . We want to use the Cheeger-Gromov com- 
pactness theorem to extract a convergent sequence of manifolds {Mk, gk)kefi- The second as- 
sumption of the theorem guarantees that there is at least one sequence of points {pk)keN £ 
whose image lies some ball of finite radius. As we have already mentioned in Section 7 of Chapter 
4, bounds on all higher derivatives of the second fundamental form imply a lower injectivity radius 
bound. We may apply the Cheeger-Gromov compactness theorem, and upon passing to a subse- 
quence, we obtain a complete pointed limit manifold {M^o, goo,Poo), an exhaustion {Uk}kGn of 
Moo, and a sequence of diffeomorphisms {4>k ■ Uk ^ Vk C Mk)km such that 4>l.gk converges 
smoothly to g^o on each compact set K C M. By induction, similar to the proof of the higher de- 
rivative estimates in Chapter 4, it follows that all higher derivatives of 4'k^k are uniformly bounded 
with respect to goo on compact sets of Moo. Passing to a futher subsequence, we obtain smooth 
convergence of (p^Fk to a limit immersion Foo on compact sets of Moo. At this stage we have 
shown the first condition in our definition of convergence on compact sets is satisfied. 

We now need to show the second condition of our definition is also satisified. The fact that 
the area bound holds on the limit is a simple consequence of the C^-convergence of the metrics. 
The remaining argument is accomplished by induction and a Cantor diagonal sequence argument. 
We begin by looking inside a ball -Bi(O) in the ambient space. Suppose that there exists no 
ko eN such that (()lFk{Uk) n BiijS) = Fk{Mk) n i?i(0) for all k > ko, then we can pass to a 
subsequence such that there exists pk G M^ such that for all k, pk ^ Vk, whilst Fk{pk) € -Bi(O). 
Passing to a further subsequence, we can assume the sequence of pointed manifolds {Mk, gk,Pk) 
converges to a limit {Moo, goo, Poo), so that there is an exhaustion {Uk}keN and diffeomorphisms 
^k '■ Uk ^ Vk <Z Mk with ^k{p) = {Pk) such that c^^^fc converges to g smoothly on compact sets 
in M. As before, by the Arzela-Ascoli theorem, passing to another subsequence, we can assume 
that i^^-Ffc converges smoothly on compact subsets to a limit immersion F^o- Now we replace Moo 
with Moo U Moo and repeat the process again. All of these components intersect with -Bi(O), and 
have area inside ^2(0) bounded below, so by the local area bound this process must stop after 
finitely may steps, and we have produced a manifold M with finitely many coimected components 
with both parts of the compactness theorem holding on -Bi(O). 

We complete the proof by induction on the size of the balls in the ambient space: If we have 
subsequence for which both parts of the theorem hold on 5^(0), then we add in more components 
if there are points in B„+i(0) thatareinFfe(Mfe) butnotin ^j^Ffe(?7fe). By the same argument, after 
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adding in finitely many components we produce a subsequence and a limit immersion satisfying 
both parts of the compactness theorem on i?„+i(0). □ 

We can use the above compactness theorem for immersed submanifolds to obtain a compact- 
ness theorem for mean curvature flows. The proof follows Hamilton's compactness theorem for 
Ricci flows. One applies the compactness theorem for immersed submanifolds at the initial time, 
and then using the higher derivative bounds and the Arzela-Ascoli theorem combined with a di- 
agonal sequence argument, one obtains a properly immersed limit solution to the mean curvature 
flow that satisfies the first condition of our definition of convergence on compact sets of x M. 
That the second convergence criterion is also satisfied again follows quickly from the curvature 
bounds. In particular, we can apply the compactness theorem to our sequence of rescaled flows 
(16.11 ) (where we have assumed the type I hypothesis). By analogy with Hamilton's compactness 
theorem for Ricci flows ([HaV]), we only require a bound on the second fundamental form itself 
(and not any higher derivatives), as the type 1 assumption ensures that all higher derivatives are 
indeed bounded above. The missing essential ingredient is the local area bound, which we shall 
address in the next section, it being a consequence of Huisken's monotonicity formula. 

Theorem 6.3 (Compactness theorem for mean curvature flows). Suppose that {Ff^, Mj^)^^^^ is 

a sequence of proper time-dependent immersions of smooth complete n-dimensional manifolds 
Mk that satisfy the mean curvature flow on the time interval I = \tQ,T). Assume the following 
conditions are satisfied: 

1) Uniform curvature derivative bounds: 

For each /c G N, there exists a uniform constant Cq such that \ hk\Fk ^ Co on x / 

2) The sequence doesn't (initially) disappear at infinity: 

There exists a time to <^nd radius R > such that Bff {0) n Fk{Mk, to) ^ 0/or all A; G N. 

3) Initial local area bound: 

For every R > there exists a constant Cr depending on R (and independent ofk), such 
that 

F-\;to){Bii{0)) 

Then there exists a subsequence (F^, Mfc)^^^ which converges on compact sets o/M^ x M to a 
complete proper time-dependent immersion (Mqo, -Foo) that is also a solution to the mean curva- 
ture flow on the time interval I. 

6.1.1. Huisken's monotonicity formula. For a fixed point {xQ,to) G W"'~^^ x M we define 
the backwards heat kernel centred at (xq, to) by 

P^oM^^t) - (4^(ig_t))n/2^^P(, 4(to-t) 

which is well-defined on X 

(— oo, to)- The centre of our backward heat kernel will most 
often be {p, T) G M"+'^ x M. Note the backwards heat kernel is defined on the ambient space and 
so we are adhering to the abuse of notation mentioned at the beginning of this chapter. Huisken's 
montonicity formula, which holds in arbitrary codimension, is the following: 

Theorem 6.4 (Huisken's monotonicity formula). Let F : S x [0, T) — ;> W^^^ be a solution of the 
mean curvature flow. For any flxed point p G the formula 

Pp,T dpgt = - [ Pp^T " ^ 



d_ 
di 



H 



2{T - t) 



dpgt < 



holds for all time < t < T. 

For each pair of times < ti < t2 < T, the monotonicity formula implies that 

/ Pp,Tdm2 < / Pp,Tdmi 
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and being the limit of a monotone sequence of decreasing functions, the Umit 
certainly exists and is finite. We shall also use the notation 



and 



Q(p) := lim e{p,t). 



Since 6 is the limit of a monotone sequence of continuous functions, it follows that Q is upper- 
semicontinuous. We refer to as the heat density and 6 as the limit heat density. An important 
property of the monotonicity formula is that it is invariant under paraboUc rescaUngs. By the 
definition of our parabolic rescaling, for each k we have 



1 f i^-pp 



(-47rs)"/2 



RecalUng that t = T + s/\\, for each fixed s G [-A^T, 0) and all k we have 

J^Pp,Tdpg, = J^pdpf;^^'^\ 

and consequently 

(6.2) lim / PriTdpat = lim / pdp^'^. 

When it is (reasonably) clear which point we are rescaling around, we will often omit the notation 
{p, T) above the measure as we have just done to reduce clutter. An important application of the 
monotonicity formula is that it provides the local area bound (independent of k) necessary to apply 
the compactness theorem for mean curvature flows. It suffices to obtain the area bound on bounded 
subintervals Ii := [—XfTjl/l] C [—XfT, 0), as the final argument will be completed by a diagonal 
sequence argument sending I to infinity. Let us fix a point p G S and some fco >> sufficiently 
large. With these choices of p and ko, then for all s G Iko and every k > ko monotonicity formula 
gives the estimate 



(47rr): 



We then compute 



JF-\Ba) JF-\Bn) 



JFk {Br) 

< e 4 / 6-4" du„'' 

]F-\Bn) 

<e^A^„/x,,^(S), 
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and thus 



We can now apply Theorem 16.3 1 to our sequence of rescaled flows : S x [— A^T, l/l] 
^n+k defined by 

Fk{p,s) = \k{F{p,T + s/\l) - q). 

We highhght that here S is fixed, and by assumption closed, however Sqo is complete and not 
necessarily compact. The existence of the limit flow (Foo,5]oo) on the time interval (— oo,0) 
follows by diagonal sequence argument letting / — oo. 

Another consequence of the monotonicty formula is the following important result, which 
enables us to pass the limit through the integral in the rescaled heat densities. The result is due 
independently to Ilmanen 11] and Stone [St], who proved it slightly different contexts. Ilmanen 
proved it in the setting of Brakke flows, while Stone proved it in the context of Huisken's original 
continuous rescaling argument. We recast their proof in our setting. 

Proposition 6.5. Let Fk ■ x [— A^T, 0) — )• be a sequence of proper mean curvature flows 
of a closed manifold S that subconverges on compact sets ofW^ xM^to a proper mean curvature 
flow Foo : Soo X (oo, 0) ^ M^, where is <^ complete manifold. Assume that for all > the 
initial submanifold satisfles the area bound 



L 



Fo (Br) 



Then for any given e > and any flxed point p G S, there exists a sufficiently large radius radius 
R = R{e, So, s) such that for each fixed s G [— A^^^T, 0) and all k > k^we have 

Proof. By localising the mononicity formula (see Ifl or ||E]) and using the initial area bound 
we see for each fixed s G [— A|r, 0) and all k > k^ that 



for every R > Rq. For every R > Rq and each fixed s £ [— A|T, 0) we estimate 



< 



„ oo 



For each fixed s G [— A^T, 0), the term on the right can be made as small as we like by choosing 
R sufficiently large, so for any given e we can fix R sufficiently large so that the desired estimate 
holds for afl i? > □ 

The proposition is, by definition, the statement that the family of weighted measures p dpg^ 
is tight for each fixed s. By Prohorov's Theorem we immediately obtain the following important 
corollary: 

lin.^l^pdp'^^=l^jdp'-<oo. 

Let us dwell for a second on why this result is important: The limit manifold Sqo we obtain from 
the compactness theorem is complete, and not necessarily compact. Certainly if Sqo contains a 



92 



6. A PARTIAL CLASSIFICATION OF TYPE I SINGULARITIES 



compact component, then this component is diffeomorphic to S by definition of the convergence. 
However, if Sqo is only complete, as it often will be, then the integral 

could very well be infinite. The fact that the weighted family of measure is tight ensures that the 
measure 'does not escape to infinity' in the limit. We remark that the -convergence of Fj^ and 
Q]^ obtained from the compactness theorem implies that ji^ 1^, that is the pushforward measures 
converge weak-* in M""'''^. 



6.2. A partial classification of special type I singularities 

In order to probe the shape of the evolving submanifold as the first singular time is approached, 
we want to rescale the monotonicity formula around the singular point p. A point p G S is called 
a general singular point if there exists a sequence of points pk ^ p and times tk ^ T such that 
for some constant 5 > 0, 

\h\\pk,tk) > 

A point p G S is called a special singular point if there exists a sequence times 
for some constant 6 > 0, 

\h\\p,tk)> 



T such that 



T-tk 

This distinction between singular points is not made in |'Hu41, and the points studied in |Hu41 are 
actually special singular points (see Defintion 2.1 of [Hu4|). The analysis to cope with moving 
points was subsequently contributed by Stone in ||Stl. We now give a partial classification of 
special type I singularities in high codimension. 

Proposition 6.6. Let S : x [0, T) — )• M"+'^ be a solution of the mean curvature flow. If the 
evolving submanifold exhibits a special type 1 singularity as t ^ T, then there exists a sequence 
of rescaled flows Ffc(S) that subconverges to a limit flow Fod^oo) on compact set of] 



pn+k 



as k - 
plane. 



oo. Moreover, Fqq : X (—00,0) 



pn+fc 



satisfles H = —l/{2s)F-^ and is not a 



Proof. The existence of the limit flow, which exists for s G (oo, 0), was shown in preceeding 
section. It remains to show the last two assertions of the proposition. Suppose the special type I 
singularly is forming at some point (p, T) G x M, so by definition there exists a sequence 

times tfc — )• T such that for some constant 5 > 0, we have \h\'^{p, tk) > jSj-- Rescaling Huisken's 
monontonicity formula at each scale = 1/ \/2(T — tk) about the single fixed point p gives 



A. 

ds 



2s 



which holds for all k and s G [— A^T, 0). For any fixed sq G 
this from sq — a to sq and rearrange a little to get 



-A|T, 0) and u > we integrate 



2s 

We take the limit as /c — )• oo, and by equation (16.21 ) and Proposition 16 .5 1 we have 



pdjjb 



9SQ 



lim 

t-*T 



P{P,T) dfJ-gt 



We then conclude, using Proposition 



s 

again, that 



pdfi. 



< oo. 



lim 

k—^oo 



so 



SO' 



2s 



so 



so- 



Hx + —F- 
2s ' 



dfi 
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and therefore Hx^ = — l/(2s)F^ on s E [sq — a, sq]. Finally, for every scale A^, at the fixed 
point p at time Sk = X\{T — tk) = —1/2 the rescaled second fundamental form satifies the lower 
bound 

I,|2 ( N \h?{P,tk) 
\h\x^ip,sk) = 72 

> 2{T - tk) ^ 



T-tk 
= 26. 

Thus the the limit flow also satisifies (p, — 1/2) > 26 and consequently it is not flat. □ 

We have just shown that the blow-up limit of a type I singularity is self-similar. In order to give 
a partial classification of these solutions, in addition to assuming that Sq satisifes |i?|min > 0, we 
also assume it satisfies the pinching condition < 4/ {3n)\H\'^. The pinching condition allows 
us to eventually reduce the problem to that of classifying hypersurfaces of a M"+^. This classifi- 
cation result was also used in the application of the strong maximum principle in Chapter |4] and 
for completeness we give a proof, adopting the proof in HCdCKl to the case of a flat background. 
We mention that this classification first appeared in [ILawl . where different techniques were used. 

Proposition 6.7. Let F : M" — )• M"+^ be an immersion of a closed manifold. If F{M) satisfies 
V/i = 0, then F{M) is of the form §P x where < p < n. 

Proof. The proof is a very nice application of the method of moving frames and Frobenius' 
Theorem. Recall from Chapter |2]that the structure equations of M"+^ restricted to the hypersurface 

F{M) are 

(6.3) duj' = -ufj A 

(6.4) uj] = -ioi 

(6.5) duj] = -ujI a - a a;!^+i , 
and that the first covariant derivative of h is 

(6.6) /ijjfc'^'' = dhij - hiiujj - hijUj\. 

Choose a local frame {ei, . . . , e„, i/} for M that diagonalises the second fundamental form. So 
hij = for i / j. If the principal curvatures are all zero, in which case M = M, the lemma is 
clearly true, and so from now on we assume that at least one of the principle curvatures is non- 
zero. If hijk = 0, then setting i = j in the above equation and using that hij = for i / j we 
get 

= dhii - 2hiiu}l, 

and since hu is symmetric in i and /, and is antisymmentric (from (I6.4I )). we have dhu = 0, so 
hii is constant. Since hijk = 0, dhij must be zero, and with these two conditions equation (16.61 1 
becomes 

= hiiujj - hijul = {hi - hj)uj, 
which shows that cuj = whenever ha ^ hjj. Thus if ha ^ hjj, equation (16.51 ) reads 

= dJj = -Lol A ujj - w^+i A ujj-^^. 

The term —ul. A Uj must be zero because / and Ujj^O would imply ha = hjj = hkk, 
which contradicts our assuption that ha ^ hjj . Therefore, 

= A 
= hikhjiuj^ A J 
= hiihjjOj'^ A oj^ . 
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We conclude that if ha / hjj, then either ha or hjj is zero, but not both. By reordering the indices 
of the frame if necessary, for each < p < n and a constant k / we have now shown 



(6.7) 



( _ _ 

. . . K/p — 

= for 1 < i < p and p + 1 < j < n. 



Now define two distributions hy = . . . = oj^ = and cj^^^ = . . . = = 0. Frobenius's 
Theorem states that a distibution cj'^ = 0, 1 < A; < n, is integrable if and only if doo^ = for every 
k. From the structure equation (I6.3l l we have dw* = — A , and so by the third equation of 
(16.7b both the distributions just defined are integrable. We therefore obtain a local decomposition 
at every point of M given by x M"^^, where Q < p <n. □ 

Let us now commence with classification in the compact case. 

Theorem 6.8. Suppose Fao '■ x (— oo, 0) — )• M"^"^ arises as the blow-up limit of the mean 
curvature flow F : T,^ x [0, T) ]R"+'^ about a special singular point. Additionally, suppose that 
Eq satisfies \H\,„in > and < 4/(3n)|/fp. If Foo{T,oo) is compact, then at time s = —1/2, 
Foo (Soo) must be a sphere S'^{m) or one of the cylinders S™'(m) x M"""^, where 1 < m < n — 1. 

Proof. Choose a local frame {cj : 1 < i < n} for S. We advise the reader that in the 
following we are making the indentification ej = F^,ej. Take inner product of H = —F-^ with H 
and differentiate in the ambient space in the chosen local frame to get 

2{VjH,H) = -{VjF,H) - {F,VjH). 

Using the Gauss relation and VjF = ej we continue to compute 

{VjH,H) = {H, {F,ep)hip) 

± 

and therefore VjH = {F, ep)hjp. A further differentiation gives 

jH — (ej, Gp)hjp + {F, hip^hjp + {F, Cp)^ phij 
— (cj) (^p}hjp {H, hip^hjp -\- (i*^, 6p)^ ph-ij 

(6.8) = hij - H ■ hiphjp + {F, ep)Vphij. 

Contracting (16.81 ) with gij gives 

AH = H-H- hiphip + {F, ep)VpH, 
and after taking the inner product with H we obtain 

(6.9) A|/f|2 = 2\H\'^ - 2^{H • hij)^ + {F, ep)VpH • H + 2| Vi^l^. 
On the other hand, contracting (I6.8I 1 with g^j we get 

(6.10) hij ■ ViVjH = - H ■ hiphij ■ hjp + {F, ep)Vphij ■ hij. 

Now recall Simons' indentity: A|/ip = 2hij ■ ViVjH + 2|V/ip + 2Z. Combining Simons' 
indentity and (16.101) gives 

(6.11) A|/i|2 = 2\h\^ + 2{F, ep)Vph,j ■ hij + 2\Vh\^ " ( 2^ h^jaKj^) - \R\^. 

Note that the term H • hiphij • hjp cancels. The idea now is to examine the scaling-invariant 
quantitiy and to do so, we first establish \H\ 7^ in order to perform the division. The 

strong elliptic minimum principle applied to equation (16.91 ) shows that either |//| = or > 
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everywhere. Since Fqo (Scxd) is assumed to be compact, it must be that \H\ > everywhere. Using 
equations (16.91) and (16.11b we compute A(|/ip/|ffp) and obtain 



Since Foo(Soo) is assumed to be compact, the function attains a maximum somewhere 

in S. At a maximum Vi(|/ip/|ffp) = and A(|/ip/|ff|^) < 0, and so at a maximum we have 

Moreover, from the basic gradient estimate (I4.6bt and the Pinching Lemma we can estimate 
(6.13) < A (}^) - ci(n)|V/i|2 - C2{n)\hi\^\h_\^ - C3(n)|L|^ 



where ci, C2 and C3 are positive constants that depend only on n. We conclude from the strong 
elliptic maximum principle that must be equal to a constant and |V/ip = = 0. 

This implies that Foo(Soo) is a hypersurface of some (n + l)-subspace of W^~^^ with covariant 
constant second fundamental form, and since was assumed to be compact, from Proposition I6.7l it 
must be a n-sphere. □ 

If Foo(Soo) is no longer compact then we cannot apply the maximum principle as we have 
just done. In this more general case, following [Hu5 |, we multiply equation (16.121 ) by VjC^I^I 
and integrate by parts. The following theorem includes the previous one as a special case. 

Main Theorem 8. Suppose F^o : x (—00, 0) — >• R"'*'^ arises as the blow-up limit of the mean 
curvature flow F : S*^ x [0, T) — )• M"^^ about a special singular point. //'So satisfles \H\,„in > 
and < 4/{3n)\H\'^, then at time s = —1/2, Foo(Soo) must be a sphere S^{m) or one of the 
cylinders S™(m) x M"^™, where 1 < m < n — 1. 

Proof. We multiply equation (16.121 ) by Vje"'^'^/^ and integrate the term involving the Lapla- 
cian by parts to achieve 

-/Soc VI I / -/ Soo II II 

2 ^2 



+ 2 / n^(^i-^^2)e 2 d^xg 



The above equation again implies that must be equal to aconstant and |V/ip = = 

and the theorem follows. 

6.3. General type I singularities 

As we mentioned in the introduction to this chapter, because the mean curvature flow in high 
codimension does not preserve embeddedness we are not able to extend Stone's hypersurface 
argument to high codimension. Let us explore a little why this is the case. Stone's result for 
hypersurfaces is the following: 

Proposition 6.9. Let F : x [0, T) — )• M"+^ be a solution of the mean curvature flow. Suppose 
that So is embedded and satisfies |-fr|min > 0. If the evolving submanifold develops a type I 
singularity at some point p £ Ti as t ^ T, then p is a special singular point. 

Stone's analysis shows that it is in fact enough to understand special singular points. We 
follow closely fSF], adapting his proof from the continuous rescaling setting to that of rescaled 
flows. Stone's argument requires the classification of special type I singularities for hypersurfaces 
obtained by Huisken in IIHu4l and IIHuSII : 
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Theorem 6.10. Let F^{T:^) C 



pn+l 



be a hypersurface that arises as a blow-up limit of the 



mean curvature flow. If Tjq is embedded and satisfles H > 0, then must be a hyperplane, the 
sphere S'^{m) or one of the cylinders S"^(m) x M"^*", where 1 < m < n — 1. 

Our equivalent theorem for submanifolds is Main Theorem |4] We also need to know that 
embeddedness of hypersurfaces is preserved by the mean curvature flow, and that the blow-up limit 
is also embedded. A proof of the former follows the next proposition, whilst for embeddedness of 
the limit we refer the reader to llMl . 



Proof of Proposition [6T91 Suppose that is developing a general type I singularity at 
some point {p, T). By definition, there exists a sequence of points Pk ^ P and times ^ T such 
that for some constant 5 > 0, 

\h\\pk,tk) > ——. 

-tk 

As before, we want rescale the monotonicity formula, but now we need to rescale about the moving 
point pk- Rescaling the monotonicity formula about the moving points pk gives 



d_ 

ds 



P 



9a 



which holds for each k and s e [— A^T, 0). For any fixed sq G [—X^T, 0) and a > we integrate 
this from sq — a to sq and rearrange a little to get 



(6.14) 



so 



so- 



F 



The difficulty now is that in general, limfc_j.oo 0{pk,tk) / Q{p). The proof is now by contradiction. 
If p is a general singular point but is not a special singular point, then by definition there exists 
some function e{t) with e{t) — as f — )■ T such that 

\h\\p,t) < 

for all time t £ [0, T). This implies that any blow-up about the single fixed point p would satisfy 
= 0. From Theorem OOl we know that a blow-up around a special singular point is one of n+ 
1 different hypersurfaces. Furthermore, the heat density function evaluated on these hypersurfaces 
takes on n + 1 distinct values, of which 1 is the smallest, which corresponds to a unit multiplicity 
plane. Full details of these calculations can be found in the Appendix of [Sfl. Crucially, since Sqo 
is also embedded, it can only be a unit multiplicity plane, and not a plane of higher mulitplicity. 
Since © is upper-semicontinuous, it is actually continuous at p, and therefore = 1 in a whole 
neighbourhood of p. Dini's Theorem on the monotone convergence of functions now implies for k 
sufficiently large, that 6{pk,tk) &{p) uniformly. This is the point at which the argument breaks 
down in high codimenion: since embeddedness of the initial submanifold is not preserved, the 
blow-up limit may be a plane of higher multiplicity, and thus G [p) could be any integer. Therefore, 
we cannot conclude that is continuous at p, and Dini's Theorem is no longer applicable. 

We complete Stone's argument: Returning now to equation (16. 14I ). for every fixed sq and every 
fixed point pk the monontonicty formula implies 



- / pdpif^'^^ < 

for aWl > k. Estimating as such, for alH > A; we have 

(•so r F- 



P 



2s 



Sending / oo and using Proposition 16.51 we obtain 



so 



0—cr JT, 



F- 



2s 



< l^pdpf^^ 



(Pk,T),Xk 



QiPk). 
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By Dini's Theorem, given any e > 0, there exists a ko such that for all k > kowe have 



P 

Sq—ct J S 



2s 



and thus 



lim 

fc— >-oo 



2s 



Using the blow-up procedure of the previous section we obtain a limit flow on (— oo, 0), and by 
Proposition 16. 5 1 the limit solution satisfies Hx^ = — l/(2s)Fx^ and is again not flat. This is a 
contradicton, since by Proposition [ 



lim 9{pk,tk) = @{p) = 1, 

fc— >oo 

which implies the limit solution is a plane and hence flat. □ 

In order to extend Stone's argument to submanifolds, we must conclude that the blow-up 
limit is a unit multiplicity plane. As in the case of hypersurfaces, it would be enough to show 
that the limit is embedded. In high codimension embeddedness is not in general preserved by the 
mean curvature flow, and unfortunately for us, a pointwise pinching condition alone does not seem 
enough to guarantee the preservation of embeddedness. We give a proof of that the mean curvature 
flow preserves the embeddedness of hypersurfaces to highlight the problem the high codimension 
introduces. 

Proposition 6.11. Let be a closed manifold, and F : TP x [0, T) —5- a solution of 

the mean curvature flow. IfT,Q is embedded, then it remains embedded for as long as the flow is 
defined. 

Proof. We follow iMl pg. 25] initially, but give an alternate argument to show that the 
distance squared d? between two points is non-decreasing in time. Let F : x [0, T] be a 
closed hypersurface, initially embedded, moving by the mean curvature flow, and suppose for a 
contradiction that T is the first time at which the hypersurface fails to be embedded. The set S of 
pairs of points {x, y), x ^ y, such that F{x, T) = F{y, T) is a nonempty closed set disjoint from 
the diagonal in S x S, otherwise fails to be an immersion at some point of S. We may therefore 
remove a small open neighbourhood -Be (A) from around the diagonal such that i?e(A) n 5" = 0. 
We consider the quantity 

5= inf inf \F{y,T)-F{x,T)\, 

ie[0,t] (p,(7)g9B,(A) 



and note that 5 is positive, since -Be(A) n 5" = and dB^{IS.) is compact. Next we claim that the 
square of the minimum of the distance function 

d^{t)= min |F(y,T) -F(x,r)P 

is bounded below by min{(i^(0), 5} > on [0,T]. This contradicts the fact that S is nonempty 
and contained in S x S \ B^{/S). To this end, if at some time d^{t) < 5, then this must occur at 
points not belonging to dB^{A), that is at points (x, y) € S x E \ i?e(A). We now want to show 
that is non-decreasing on S x S \ B^{A) x [0, T], which proves the claim and the theorem. We 
compute the first, second and time derivates of d"^ in some choice of local coordinates {x*} near x 
and {y*} near y. The first derivatives are 



dyjd^ = 2{y - x,d. 



d^jd^ = —2{y — X, d. 



1 I 1 
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the second derivatives 

dy.dyjcf = 2gf- - 2{y - x, h\.Uy) 
d^.d^icf = 2gfj + 2{y - x, hfji^^) 
^^^^yJ(f = -2{^^^,^yJ)■, 

and last of all the time derivative is 

dtcf = 2{y - X, -HyVy + H^j^v^). 

Importantly, observe that at a minimum of the distance function the tangent planes at x and y are 
parallel to each other. We may therefore choose local coordinates such that {x*} and {y*} are 
parallel for each i. We now compute 

dt V^" dx'dxi ^ dy'dyi ^ ^ ' ^' ' dx'dy^ J 

= 2{y - x, -HyVy + H^Vx) - 2n - 2{y - x, H^v^) -2n + 2{y - x, HyVy) + An 
= 0. 

We conclude by the maximum principle that is non-decreasing in time. □ 

In the above proof it was crucial that we were able to choose parallel orthonormal frames at 
the points x and y : without the good 4n contribution from the cross-derivative terms the proof 
does not work. In high codimension this is not possible to do in general since the tangent planes 
could easily be orthogonal to each other at a point of minimum distance, which results in zero 
contribution from the cross-terms. 

6.4. Hamilton's blow-up procedure 

In the above singularity analysis the type 1 assumption was essential to obtain a smooth blow- 
up limit. It is tempting to think that the above blow-up analysis could be simplified by using a 
Hamilton blow-up argument, in which case even for a type II singularity we could obtain a smooth 
blow-up limit. If one performs a type II Hamilton blow-up and uses this in combination with the 
monotonicity formula in a similar fashion to what we have done above, then one again obtains a 
smooth limit solution to the mean curvature flow that satisfies \H\ = —F. The subtle problem 
with doing this is that the point of maximum curvature may not actually lie in this limit. As an 
example, blowing-up the grim reaper in such a fashion would in fact result in a cylindrical limit. 

Next we want to give a more successful application of a Hamilton blow-up to give a short proof 
of the limiting spherical shape of the evolving submanifolds considered in Main Theorem 2. The 
interested reader may like compare the following with the corresponding argument in the Ricci 
flow, which can found, for example, in [T|. Here the Codazzi equation performs the same role as 
the contracted second Bianchi indentity, and the Codazzi Theorem that of Schur's Theorem. For a 
proof of the Codazzi Theorem we refer the reader to [Sp. Thm. 26]. To begin, pick any sequence 
of times {tk)k&n such that ^ T as k ^ oo. The Pinching Lemma implies that and \H\'^ 
have equivalent blow-up rates, so we can in fact rescale by \H\'^. Then, since S is assumed to be 
closed, we can pick a sequence of points {pk)keN defined by 

\H\{pk,tk) = max\H\{p,tk). 

For notational convenience, set := \H\(pk,tk). We now define a sequence of rescaled and 
translated flows by 

Fk{q, s) = Afe {F(q, tk + s/Xl) - F{pk,tk)) , 

where for each k, Fk : T, x [A|r, 0] is a solution of the mean curvature flow (in the 

time variable s). The second fundamental form of the rescaled flows is uniformly bounded above 
independent of k and we can apply the compactness theorem for mean curvature flows to obtain a 
smooth limit solution of the mean curvature flow F^c : Eqo x {—oo, 0] — )• R"+'^'. Futhermore, at 
s = the limit solution satisfies = 1 by construction, so the limit is not flat. By definition 



6.4. HAMILTON'S BLOW-UP PROCEDURE 



99 



of the rescaling, the second fundamental fomi rescales as = |/ip/A^, and so estimate (14.13b 
of Chapter IHrescales as 

\h\\<C,X-,'\H\\. 

The limit therefore satisfies 

(6.15) |/i|L=0' 

and thus Foo{^oo) is a totally umbilic submanifold. By the Codazzi Theorem, Foo(Soo) must 
be plane or a n-sphere lying in a (n + 1) -dimensional affine subspace of R"+'^. We know that 
l-^Uoo = 1' ^rid so Foo(Soo) is not a plane. □ 
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